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SUMMARY 


A user's manual for the computer program developed for the prediction of 
propeller-nacelle performance reported in Volume I, "An Analysis for High 
Speed Propeller-Nacelle Aerodynamic Performance Prediction - Theory and 
Initial Application" is presented. The manual describes the computer program 
mode of operation requirements, input structure, input data requirements and 
the program output. In addition, it provides the user with documentation of 
the internal program structure and software used in the computer program as it 
relates to the theory presented in Volume I. Sample input data setups are 
provided along with selected printout of the program output for one of the 
sample setups. 
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INTRODUCTION 


The purpose of this manual is to provide the user with sufficient 
documentation to run the computer code for the propeller-nacelle performance 
analysis (PANPER) developed by the United Technologies Research Center (UTRC) 
under contract to the NASA Lewis Research Center. The computer analysis is 
capable of predicting the performance for high speed propeller-nacelle config- 
urations for either single or coaxial counter-rotating propellers for either 
an internal flow condition (wind tunnel) or external flow conditions (free 
flight). The analysis was developed by combining and modifying existing 
propeller performance prediction capabilities applicable to the high speed 
flight problem with an existing axisymmetric through flow analysis modified to 
calculate external flow problems. The resulting combined analysis couples the 
separate solution procedures by including in each solution portion in a con- 
sistent manner the appropriate effects due to the respective portions of the 
two solution procedures. The basic program structure is shown in the flow 
diagram of figure 1. The propeller performance solution is obtained including 
the influence of the nacelle body, and the flow field and nacelle performance 
are calculated including the influence of the work done by the propeller on 
the fluid. Either of the performance solutions (propeller or nacelle) can be 
obtained without the other if so desired. The technical aspects of the solu- 
tion procedure are detailed in reference (1) and will not be explained here. 
This manual has been written under the assumption that the leader has read the 
technical report (reference 1) and is familiar with the theoretical features 
of the analysis. 

This manual consists of two major portions: 1) a description of the 

setup and input data required to run the computer code; and 2) documentation 
of the internal program software as it relates to the technical aspects of the 
analysis • 

The input portion is broken into three sections, describing the basic 
program setup requirements and program mode operation, the propeller related 
setup and input and the nacelle related setup and input, in this order. The 
documentation portion consists of two sections which describe the subroutines 
and labeled common blocks used in the computer program for the propeller and 
nacelle portions of the analysis, respectively. 
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DESCRIPTION OF THE PROGRAM OPERATION 


This section is intended to describe the general features, program setup 
and input data of the PANPER computer program in sufficient detail so that the 
program can be operated successfully by the user. Input to this p'rogram 
consists of three parts; the program mode control data, propeller data and 
nacelle data, in this order. The input data for each of these parts will be 
described in the following subsections. 

The first subsection describes the various modes that the program will 
operate in and the input control switches. Special attention should be paid 
to these mode switches since this program may in effect solve one of three 
problems : 

(1) Propeller Lifting Line Analysis only 

(2) Nacelle Analysis only 

(3) Combined Propeller-Nacelle Analysis 

The second and third subsections present a detailed description of the input, 
output and diagnostics of the Propeller Lifting Line Analysis and Nacelle 
Analysis portions of the program, respectively. 

The computer program was written and developed in FORTRAN V Computer 
Language for use on a UNIVAC 1110 computer. Before execution of the PANPER 
Program, fourteen files must be assigned in the JCL RUNSTREAM. Information 
about these files is given in Table (l). Sample Runstreams for three cases 
are presented in Appendix A, and selected output for the second case is 
presented in Appendix B. 
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Input Mode Control 


This subsection describes the mode control input data card required to 
operate the PANPER Program. Information from this card will determine which 
mode of operation the PANPER Program will perform in. On this card the Mode 
of operation control information is read in as described below: 

Name Column Format Comments 

N0PPF 1-2 12 N0PPF = -1, The propeller analysis is performed 

without including the nacelle effects calculated 
directly from the nacelle portion of the program. 
N0PPF = 0, The nacelle analysis is performed 
uncoupled from the propeller lifting line code. 

The blade forces may be considered through input. 
N0PPF = 1, The propeller and nacelle analysis is 
performed through coupling of the nacelle portion 
and the propeller lifting line portion of the code. 

N0PPC 2-4 12 Indicates number of passes through the viscous flow 

algorithm of the nacelle portion and the propeller 
portion. See reference 1, Section entitled: 
Description of the Combined Analysis Solution 
Procedure. 

It should be observed that for N0PPF = 0, the propeller data will not be 
read in and for N0PPF = 1, the nacelle data will not be read in. Descriptions 
of the propeller and nacelle data are given in the following subsections of 
this section. Samples of this card can be seen in Appendix A. Finally, it 
must be noted that N0PPC is a cycle counter on the number of passes through 
the viscous propeller-nacelle flow solution. It may be desirable to cycle 
through the viscous flow solution and propeller solution until the propeller 
blade forces do not significantly change. However, experience to date has not 
demonstrated a necessity to perform this cycle for propeller performance 
applications . 
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Propeller Portion 


Major Input Features 

The major input features consist of four basic groups of input data along 
with the appropriate program control data. These four groups of data are 
blade geometry, airfoil characteristics, inflow properties and the wake 
geometry. In the following subsections, these groups of input are described 
in some detail so that the user will understand their importance. 

Bljguie t_in£ L.in_e_G£ome_tr_^ 

The geometric description of the blade lifting line representation of the 
blade is of primary importance for obtaining the most accurate solutions given 
the assumptions inherent in the analysis. The hub-pitch axis centered 
cartesian coordinates for each lifting line segment boundary must be input 
(XSB, YSB, ZSB) consistent with the input blade twist distribution (THET) so 
that the program can correctly rotate this geometry about the pitch axis for 
the required blade angle. The coordinates for the definition of the blade tip 
for the tip Mach cone calculations (XMC, YMC, ZMC) must also be input 
consistent with the twist distribution. For counter rotating coaxial 
propellers, each set of coordinates is input referenced to its respective hub 
and pitch axis centers. The technical description of this coordinate system 
is detailed in reference 1, see figure 2. The selection of the blade 
segmentation is of primary importance in regions of severe loading gradients. 
For these regions (generally the tip of the blade) finer segmentation is 
required as compared with regions of weak gradients. 

Propeller^ Blade Ai_rfoij_ Characteristics 


To calculate the blade airloading, it is necessary to specify the 
distribution of airfoil type along the blade radius. There are three 
different sets of airfoil data available in this analysis: two sets of NACA 

16 series isolated airfoil data, whose characteristics are described in 
reference 1, and one cascade airfoil data set for NACA 65 series, also 
described in the above noted reference. The user specifies the use of these 
airfoil data sets by specifying the radial location which denotes the outer 
boundary of the region (RADCAS) for which it is desired to use the cascade 
data set. Outboard of this region, the distribution of the identification 
number (23 or 24) for the isolated airfoil data sets is input through the 
airfoil type designation number distribution input (AIRN). If desired, it is 
possible to model the cascade effects on the isolated airfoil data by 
application of an analytical cascade correction (CASCAD). This model is also 
described in reference 1. The use of this model may be desirable for two 
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reasons: first, if the inboard section of the propeller blades is not 

adequately modeled with NACA 65 series airfoil sections, and second, if the 
cascade influence extends beyond the region where the NACA 65 series airfoil 
types apply. 

Once the distribution of airfoil type and cascade regions are* determined 
for the design under consideration, the particular airfoil characteristics are 
defined by additional input. These characteristics are: the design lift 
coefficient (DECL), the thickness to chord ratio (T0VC) , and the chord 
(C0RD). 

_In^f l_ow_P_ro^£er^tJ : es_ at__th_e_B^ad_e_R^)w 

This analysis allows the user the ability to describe the noninduced 
inflow properties at the propeller blade rows if run independent from the 
nacelle portion of the analysis. It is therefore possible to prescribe the 
nacelle's influence (or any desired influence) on the inflow conditions at the 
propeller blades without running the nacelle portion of the analysis. This 
may be desirable if the variation of the nacelle's influence is small for 
slight changes in the propeller designs. The inflow properties are the axial 
(V0VO) and radial (URVO) noninduced inflow velocity ratio distributions, and 
the density (DENS) and speed of sound (S0UN) ratio distributions along the 
blade radius. These distributions scale the respective freestream values to 
define the local inflow conditions at the blade rows. 

Wake_M£d£l_D£S£r_ip£ic» n 

The description of an accurate wake geometry for the flight condition 
under investigation is of primary importance for accurate predictions of the 
induced inflow solution and the resulting propeller blade air loading. The 
wake models available have been described in detail in the technical section 
of reference 1; however, a brief review of the applicable wake models for the 
different flight regions follows. 

For static thrust conditions, the generalized wake model should be used 
(figure 3). It has been clearly demonstrated to be the most accurate model 
available and is necessary for accurate performance solutions. Wake rollup 
modeling must also be used for this flight condition (figure 4). In low speed 
flight conditions, the classical (figure 3) or modified classical wake model 
(standard model for high speed flights) is probably sufficient for reasonable 
performance predictions; however, it is clear that the wake model must have 
some of the features of the generalized wake model (radial contraction, in 
particular). Because of this, it is possible to use the generalized wake 
model for nonzero flight speed conditions. In this case the generalized wake 
model will have the inflow velocity distribution superimposed on it to 
describe the wake geometry. Thus, with careful selection of the input gener- 
alized wake coefficients, it is possible to model a low speed wake geometry if 
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the required characteristics are known. Wake rollup modeling should probably 
be used for these flight conditions. For high speed flight conditions, the 
wake is carried away from the propeller so rapidly that it is doubtful that 
any model other than the classical or modified classical wake will be 
required. Generally no wake rollup modeling is required at these flight 
speeds. 

Similar considerations must be given to the influence of the nacelle on 
the wake geometry. For static thrust conditions, the nacelle influence cannot 
be modeled by using the nacelle portion of the analysis. However, if the 
displacement of the wake due to the presence of the nacelle is known, it can 
be modeled through the wake geometry input option. For all other flight con- 
ditions, the nacelle's influence can be included directly in the analysis. 

If it is desired for any reason to use a wake model which is not geomet- 
rically compatible with the basic wake models available, the wake geometry can 
be input in cylindrical coordinate form. This allows for a wide range of 
possible wake modeling capabilities in this analysis. 
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Detailed Description of Propeller Data Input and Setup 

^t£n_da£d^I_npjjt^_D_at_a and_ Se_tu£ 

The propeller input data is grouped into 3 distinct data sets, the first 
set consists of input data which describes the propeller analysis modeling 
options, freestream flight conditions and primary propeller characteristics. 
The second set of data defines the physical location of the blade lifting line 
segment boundaries and the coordination used to define the location of the tip 
Mach cone. These items are referenced to their respective centers of rota- 
tion. The third set of data is used to describe the local blade element 
flight inflow properties (based on the freestream flight condition) and 
secondary blade characteristics. This data set consists of interpolation 
tables for the required items. All of these input data sets are described in 
the following subsections. For coaxial propellers, the second and third data 
sets are repeated for the second propeller following all of the data sets for 
the first propeller. Two sample propeller input data decks are listed in 
Appendix A (case 1 and case 3) for an isolated propeller mode and a combined 
coaxial propeller-nacelle mode. Each data set is initiated by a header card 
with an alphanumeric label in card columns 1 through 6 (left justified) and 
terminated by a card with the alphanumeric label END in card columns 1 through 
3. Within a given data set, there is no ordering dependency for the input 
items and if duplicity of the item occurs, the last value read will be used. 
All numbers are input in FORTRAN floating point or exponential format. 

Da_ta_Set_I 

The header card for this data set consists of the characters INPUT in 
card columns 1 through 5. The input data required for this data set is input 
one card at a time following the header card. Each card has a label and value 
punched on it. The label designates the item and the value for that item 
follows the label on the card. The label is an alphanumeric, three to six 
character name, left justified in card columns 1 through 6, with the input 
value following in a FORTRAN E20.8 format (columns 7 through 26). The 
description of the labels and corresponding input data items are listed below. 
For items with the numeric characters 1 or 2 on the end of the label, the 1 
and 2 designate the first and second propeller quantities respectively for a 
coaxial condition. If the input item is omitted, a value of 0.0 is used 
internally. 

Input Label Description 

BLADE 1 ,BLADE2 Blade number per propeller 

CASCAD Option switch to use an analytical cascade correc- 

tion on isolated airfoil data, a value of 1.0 
requests the model based on flat plate theory, a 
value of 2.0 uses a model based on empirical 
correlations of reference 44 in Volume I. 
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Input Label 


Description 


CBWAKE 


CNSECT 


C0FL0W 


C0MPRS 


CPI 


CTI 


DCPDT 


Option control for including the effects of compres- 
sibility on the induced velocity calculation from 
the bound lifting line vortex (this model is of 
questionable validity). A value of 1.0 sets this 
option, a value of 2.0 sets this option but the 
bound influence on the blade generating the effect 
is neglected. Should use 0.0. 

Fraction of the chord measured from the leading 
edge, used to determine the tip Mach cone inter- 
section location on the blade for the Evvard tip 
correction, generally the trailing edge (1.0) is 
used. A zero input sets this value to 1.0. 

Option control for overriding the limitation that 
the wake and bound vortex compressibility effects be 
applied only when the section Mach numbers are 
greater than 1.0. (This model is of questionable 
validity.) An input value of 1.0 engages this 
option. Should use 0.0, see section entitled: 
Compressibility Considerations for Induced Velocity, 
of reference 1 . 

Option control for including the effects of compres- 
sibility on the induced velocity calculation from 
the trailing wake geometry. A value of 1.0 sets 
this option. This option should be used. 

Requested power coefficient for performance iter- 
ation. A zero value assumes no iteration. This 
iteration option will not work for coaxial pro- 
pellers. 

Requested thrust coefficient for performance itera- 
tion. A zero value assumes no iteration. This 
iteration option will not work for coaxial propel- 
lers. The power iteration will override this option 
if both are requested. 

The derivative of power coefficient with respect to 
blade angle. If the power coefficient iteration is 
requested and this value is nonzero, the input value 
is used to determine the second iteration blade 
angle value, otherwise a change of 1.5 degrees in 
blade angle is used for the second iteration. 
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Input Label 


Description 


DCTDT 


DEBUG 


DENSTY 

DFRNAC 


DPRNAC 


DPSI 


EVAARD 


G400 


The derivative of thrust coefficient with respect to 
blade angle. If the thrust coefficient iteration is 
requested and this value is nonzero, the input value 
is used to determine the second iteration blade 
angle value, otherwise a change of 1.5 degrees in 
blade angle is used for the second iteration. 

Intermediate print option control, generally not 
used (0.0). A value of 1.0 requests printout of 
many quantities associated with the geometry trans- 
formations, geometric influence coefficients, circu- 
lation matrix and intermediate aerodynamic 
quantities. A value of 2.0 requests a full debug 
printout and should not be used. 

Freestream air density (slugs/ft 3 ) 

Input skin friction drag (lb^) due to the nacelle. 
Included in the performance calculation if input A 
positive value is opposite the direction of positive 
thrust. 

Input pressure drag (lbf) due to the nacelle. 
Included in the performance calculation if input. A 
positive value is opposite the direction of positive 
thrust . 

Maximum size of the azimuth increment (degrees) 
allowed to define the wake geometry and azimuthal 
interval in either single or coaxial mode. The 
program internally calculates the actual value. 

Option switch to request tip relief model. An input 
value of 1.0 requests tabled values for the Evvard 
model be used, a value of 2.0 requests that a func- 
tional form for the Evvard model be used which is 
slightly different than the tabled values, a value 
of 3.0 requests that conical flow theory be used 
with a variable Mach number distribution, while a 
value of 4.0 uses a fixed Mach number. 

Option switch to couple via an external file to an 
aeroelastic response analysis. If nonzero, the 
value identifies the device unit number to be used. 
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Input Label 
HUBQ1, HUBQ2 

PRINTI 

PRMAT 

PRNT0P 

PR0PT 

PR0PNM 
RADI, RAD2 

RDCASl , RDCAS2 

RDTRN1 , RDTRN2 

REV 


Description 

Input hub torque (ft-lb^). This input value will 
be included in the performance calculations and 
performance iteration loops. A positive value 
represents a power loss which the engine must over- 
come. However, the fluid does not sense this loss. 

Option to delete vector input listing, nonzero to 
perform this function. 

Geometric Influence Coefficient print option, 
generally not used (0.0). A value of 1.0 requests 
the printout of the geometric influence coefficients 
used to compute the induced velocity in both the 
cylindrical coordinate system and the blade element 
coordinate system. 

Option to delete performance printout of spanwise 
distribution quantities, nonzero to perform this 
function. 

Wake geometry print option, generally not used 
(0.0). A value of 1.0 requests that the wake 
coordinates be printed. 

Number of propeller blade rows (1 or 2). 

Input blade radius (along the pitch axis), this 
value may not be the true radius if the blade is 
swept off of the pitch axis (ft). 

Outermost fraction of the blade radius for which the 
cascade airfoil data will be applied. If zero no 
cascade airfoil data is used. 

Maximum radius to which the airfoil transition 
interpolation model can be applied. This value is 
also the flag which requests this option. 

Number of revolutions of wake geometry used to model 
the actual wake. The value chosen should be 
sufficient in length to approximately model an 
infinite wake's influence. Low flight speeds 
require a larger number of revolutions of wake 
geometry than high speed conditions. For high speed 
conditions use 2.0. 
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Input Label 
R0LUP1 , R0LUP2 

RPMRF1, RPMRF2 

RPM1, RPM2 
i SKIN0P 

S0UND 

I 

i STACK 

i 

STN 

TAUEXP 

THETA 1, THETA2 


Description 


Option switch to model trailing wake rollup, A 
nonzero value requests that the input value repre- 
sents the number of outer filaments to be rolled up 
into the tip vortex filament at a specified azimuth 
position behind the blade. When this model is 
requested, the remaining filaments are implicitly 
rolled up into a remaining or root vortex. See 
figure 4. The value to use should correspond to the 
maximum circulation location. 

Reference rpm for twist increment due to steady 
airloads . 

Propeller rotation speed (rpm). 

Option switch which requests a skewed flow drag 
model. An input value of 1.0 requests this option. 

Freestream speed of sound (fps). 

Fraction of chord measured from the leading edge to 
define the position of the lifting line on each 
blade element, generally the quarter chord line is 
used (0.25). 

Number of inflow stations per blade. Maximum of 15. 
Generally at least 10 are used. 

Exponent for airfoil transition interpolation 
function. 

Input reference blade angle (degrees). This refer- 
ence angle rotates the input twist distribution 
about the pitch axis. If the performance iteration 
is requested this value will be changed internally 
and the input value is the first iteration value. 
Positive leading edge up. 


12 



Input Label 
TRUCIl, TRUCI2 


TRUCTl, TRUCT2 


TYP CAS 


Description 


Azimuth position behind the blades for which the 
root rollup occurs, a zero value assumes rollup 
starts immediately at the blade. Because there is 
generally no root vortex formed, a large* value 
should be input when rollup is requested (degrees). 

Azimuth position behind the blades for which the tip 
rollup occurs, a zero input assumes rollup starts 
immediately at the blades (degrees). 

Option switch to select cascade type. An input of 
0.0 requests no cascade data be used. A value of 

1.0 uses the correlation from reference , while a 

value of 2.0 requests the correlation of reference 


VIM0M1, VIM0M2 


VKTAS 

V0RC0R 


WAKE0P 


WAKNAC 


Input momentum induced velocity, used to define the 
wake geometry. If a performance iteration is 
requested, this value is internally corrected to 
match the resultant performance (fps). 

Freestream flight velocity (knots). 

Fraction of the blade radius to define a vortex core 
for geometric influence coefficient calculations. 
Generally 10 percent of the chord is used. 

Option control for wake model selection. A zero 
value requests the standard wake model (modified 
classical wake) defined by the momentum— induced 
velocity and the radially varying input axial inflow 
velocity distribution be used. A value of 1.0 
requests a wake model defined by the flight velocity 
and momentum input velocity be used (classical 
wake). A value of 2.0 requests that the wake 
geometry be input to the analysis and a value of 3.0 
requests that the wake coefficients for the general- 
ized wake model be input. 

Option control for including the effects of the 
nacelle on the wake geometry through the use of a 
displacement correction to the requested wake model. 
The option generally requires that WAKE0P = 1.0 so 
that double accounting of the nacelle's influence on 
the wake geometry does not occur. An input value of 
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Input Label 


Description 


1.0 requests this option. If the standard wake 
model (WAKE0P = 0.0) is requested with this option, 
the program execution will be terminated because of 
this double accounting of the nacelle's influence. 

If it is actually desired to use the standard wake 
model, this feature can be overridden by adding to 
the program input, immediately following the input 
data, a card with the alphanumeric characters 0VER 
in card columns 1 to 4 . 

ZHUB Nondiraensional (radius of propeller one) displace- 

ment between the propeller disc centers for coaxial 
propellers. A positive value places the second 
propeller behind the first. Must be consistent with 
input for nacelle portion of the analysis. 

Da_ta__S_et__I_I 

The header card label for this data set is BLADE, in card columns 1 
through 5. The first input after the header card must be the integer value 
for the number of blade segment boundaries, free field format. Following this 
input card, the data items are input. For each set of blade segment boundary 
items (STN+1), a labeled header card is input with the alphanumeric labels 
described below, followed by the free field formatted vector (root to tip) 
on the next card for the item in question. For the tip Mach cone definition 
quantities, this format is identical but the vector item is replaced by a 
single value. All of these items should be input. 

Input Label Description 


XSB Input cartesian coordinate vector, X, inboard to 

outboard, to define the blade lifting line segment 
boundaries. Nondimensionalized by the blade segment 
radius boundary value (RADI). Maximum of 16 
boundaries (15 segments). 

YSB Input cartesian coordinate, Y, to define the blade 

lifting line segment boundaries, nondimensionalized 
by the last segment boundary value (RADI). Maximum 
of 16 boundaries (15 segments). 


Free field format consists of a series of numbers (Fortran floating point 
or exponential) separated by commas. If more than 80 card columns are 
needed for an input vector, the vector continues on the following card. 
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Input Label 


Description 


ZSB Input cartesian coordinate, Z, to define the blade 

lifting line segment boundaries. Nondimensionalized 
by the last segment boundary value (RADI). Maximum 
of 16 boundaries (15 segments). 

XMC Input cartesian coordinate, X, to define the blade 

leading edge tip location for the tip Mach cone 
definition. Nondimensionalized as noted above. 

YMC Input cartesian coordinate, Y, to define the blade 

tip location for the tip Mach cone definition. 
Nondimensionalized as noted above. 

ZMC Input cartesian coordinate, Z, to define the blade 

tip location for the tip Mach cone definition. 
Nondimensionalized as noted above. 


I)aJ^a_Set_III 

The header card for this data set is labeled VARDAT. The input inter- 
polation tables for each item in this data set are input with a header card 
with the label for the particular item on it, followed on the next card by the 
integer number of interpolation stations in free field format (minimum of 4, 
maximum of 20). The independent vector (non-dimensional X-wise coordinate) 
for the particular item follows on the next card (root to tip) in free field 
format. The dependent vector values then start on a new card following the 
independent vector in the corresponding order. An example of the format for 
one input item follows: 



Format 

Label 

(Alphanumeric) 

N 

(Integer) 

» X2 9 . • . , x n 

(Floating Point) 

> ^2 9 ^ 3 * • • • * ^n 

(Floating Point) 

The required input items 

are described below. 
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Input 

AIRN 


C0RD 

THET 

DECL 

T0VC 

DENS 

S0UN 

URVO 

V0VO 

BETA 


Label 


Description 


Airfoil type designation number distribution. There 
are only two values for input, 23.0 which requests 
the Manoni airfoil data tables and 24.0 which 
requests the NACA airfoil data tables (reference 1). 
Because the values used internally are computed by 
interpolation from this input vector and then 
converted to integer values, the input values of 
23.1 and 24.1 are generally used to guarantee that 
the integer values of 23 and 24 are used internally. 
These values start at the hub even if the cascade 
data is used. 

Chord distribution in feet. 

Built in blade twist distribution in degrees. 

Leading edge up (direction of positive rotation) is 
positive. 

Design lift coefficient distribution. 

Airfoil thickness to chord ratio distribution. 

Local blade row density to freestream density ratio 
distribution. Overridden if Nacelle portion of the 
analysis is used. 

Local blade row to freestream speed of sound ratio 
distribution. Overridden if Nacelle portion of the 
analysis is used. 

Local blade row radial inflow velocity to freestream 
velocity ratio distribution. Overridden if Nacelle 
portion of the analysis is used. 

Local blade row axial inflow velocity to freestream 
velocity ratio distribution. Overridden if Nacelle 
portion of the analysis is used. 

Dynamic twist distribution, internally scaled by the 
ratio of rpm to reference rpm squared. Incremen- 
tally added to static twist distribution, degrees. 
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Optional Generalized Wake Geometry Input Coefficients (WAKE0P = 3.0) 


In order to maintain flexibility with regard to the generalized wake 
model, the generalized wake geometry equations are included in a separate 
subroutine which requires the input of a set of generalized wake coefficients 
if this model is requested. In this subroutine the input wake parameters are 
applied to the wake equations to compute the wake filament coordinates. The 
input instructions for the wake geometry subroutine (RWZW7) are included 
herein. 


The wake equations for the generalized wake model, containing the input 
generalized wake coefficients, and graphs showing the applicable wake regions 
of the equations are presented in figure 5 (in which program symbols are 
used). The designations r * 0 and r * 1 indicate nondimensional radial 
coordinates at the axis of rotation and at a distance of one propeller radius 
from the axis of rotation, respectively. The wake representation is also 
explained in reference 2; however, for the wake equations therein: (1) AK30 

and AK31 are not included, (2) it is assumed that AK10 is zero and (3) the 
axial coordinate for the tip vortex and the vortex sheet extension to r s 1 is 
relative to the blade tip instead of the propeller hub. 


Input for Generalized Wake Geometry 


Card Program 

No. Column Symbol Description of Input Item 


1 9-10 I0PT Option for the vortex sheet boundary within 

a wake azimuth of 360. /BL and the blade 
(fixed point, right adjusted). Normally, 
set I0PT s 1 to establish a parabolic vortex 
sheet boundary through: (1) the origin of 
the outermost vortex sheet filament at the 
blade, (2) the rolled up tip filament 
coordinates at an azimuth of 360. /BL and (3) 
the intersection of the vortex sheet at an 
azimuth of 360. /BL and the tip vortex 
boundary. IF I0PT =* 0, a linear vortex 
sheet boundary is established between (1) 
and (3) above. See reference 2 for more 
detail . 


11-20 A Curve fit constant. A, in the tip vortex 

radial coordinate equation (see figure 5). 

21-30 LAMBDA Curve fit constant, LAMBDA, in the tip 

vortex radial coordinate equation (see 
figure 5). 
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Column 


Card 

No. 


Program 

Symbol 


31-40 PHINPO 


41-50 PHINP1 


2 1-10 AK1T 


Description of Input Item 

Wake aximuth angle, PHINPO, that separates 
the axial velocity regions AK20 and AK30 for 
the vortex sheet extension to r *= 0, degrees 
(see figure 5). 

Wake azimuth angle, PHINP1, that separates 
the axial velocity regions AK21 and AK31 for 
the vortex sheet extension r * 1, degrees 
(see figure 5). 

Axial velocity of the tip vortex between 
the blade and the passage of the following 
blade at wake azimuth 360, /BL (nondimen- 
sionalized by rotor tip speed; negative 
down) . 


11-20 AK2T 


21-30 AK10 


31-40 AK20 


41-50 AK30 


51-60 AK11 


Axial velocity of the tip vortex after the 
passage of the following blade at the wake 
azimuth 360. /BL ( nondimens ionalized by rotor 
tip speed; negative down). 

Axial velocity of the vortex sheet extension 
to the center of rotation in the wake 
azimuth region between the blade and the 
passage of the following blade at the wake 
azimuth 360. /BL (nondimens ionalized by rotor 
tip speed; negative down). 

Axial velocity of the vortex sheet exten- 
sion to the center of rotation in the wake 
azimuth region between the passage of the 
following blade at the wake azimuth 360. /BL 
and the wake azimuth PHINPO (nondimens ion- 
alized by tip speed; negative down). 

Axial velocity of the vortex sheet extension 
to the center of rotation following the wake 
azimuth PHINPO (nondimens ionalized by tip 
speed ; negat ive down) . 

Axial velocity of the vortex sheet extension 
to r 3 1 in the wake azimuth region between 
the blade and the passage of the following 
blade at the wake azimuth 360. /BL (nondimen- 
s ionalized by rotor tip speed; negative 
down) . 
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Card Program 

No, Column Symbol 

61-70 AK21 


71-80 AK31 


Description of Input Item 

Axial velocity of the vortex sheet extension 
to r * 1 in the wake azimuth region between 
the passage of the following blade at the 
wake azimuth 360. /BL and the wake azimuth 
PHINP1 (nondiraensionalized by tip speed; 
negative down). 

Axial velocity of the vortex sheet extension 
to r = 1 following the wake azimuth PHINP1 
(nondiraensionalized by tip speed; negative 
down) . 


Optional Input Wake Geometry (WAKE0P g 2.0) 


If desired, an arbitrary wake geometry model may be used in the analysis 
by input of the complete wake geometry. The only constraints for the wake 
model are that the description of the geometry is assumed identical for each 
blade (a required assumption in the solution procedure) and that the 
coordinates be input in the cylindrical coordinate system for equally spaced 
wake azimuth positions consistent with the blade azimuth increment and inflow 
station boundaries. This model allows for the most exact description of the 
wake geometry, if known, given the inherent assumptions of the analysis. The 
description of the input format follows. 

The wake geometry is input in separate sets for each trailing wake fila- 
ment, inboard to outboard. Each set contains subsets for each revolution of 
wake geometry requested. Each subset will start with a new card. The radial 
and axial coordinates are paired (radial, axial) for the trailing wake segment 
boundary for each wake azimuth position (wake age) starting at the youngest 
and ending with the oldest wake segment boundary. Thus, the wake azimuth 
position is implicitly assumed to be consistent with the input blade azimuth 
(DPSI) increment and the number of trailing segments must be constant with the 
number of wake revolutions. The FORTRAN format used is 10F8.4 for each card 
of data. There are KT0T (number of trailing filament) sets of cards. Each 
set of cards contains NT0T (number of wake revolutions) subsets, with each 
subset containing JT0T1 (number of wake azimuth positions per revolution + 1) 
pairs of wake coordinates. Because each subset of data will contain a 
complete definition of a revolution of wake geometry, the first coordinate 
pair of each subset will be identical to the last coordinate pair of the 
previous revolution, excepting the first subset. The total number of input 
data pairs is then (KT0T)x(NT0T)x(JT0Tl) . For a typical high speed condition 
using two revolutions of wake, 12 segment boundaries and a blade azimuth 
increment of 15 degrees, the number of pairs would be ( 12)x(2)x(25) = 600 or 
1200 single values. 
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Description of Propeller Solution Output 


The propeller output section can be broken into three distinct portions, 
initial, intermediate and final output. The program user has a large number 
of print options which control the amount of intermediate output. The initial 
and final output are not optional. The descriptions of the output quantities 
are presented in the following sections. A sample printout for selected por- 
tions of the propeller solution portion is presented in Appendix B. It should 
be noted that in the description to follow there is only one propeller and 
propeller position for a single propeller configuration. For coaxial 
propellers the intermediate output is repeated for each propeller. 

Initial^ ()ut_put 

During the reading of the propeller input data, the data as read in is 
immediately printed out. This information is entitled: PRINTOUT OF INITIAL 
DATA AS READ IN, and if the program execution terminates during the reading of 
an input data item, the user will see the item which was last read before 
program termination. This feature has two advantages: first, if an incorrect 

item is attempted to be read in, the user can quickly determine the incorrect 
item; and second, a complete listing of the propeller input data as used in 
the analysis is available for later review if desired. This output section 
always occurs during the input of the propeller data. If the combined 
analysis (propeller and nacelle) mode is being used, this output occurs long 
before (precedes the nacelle inviscid solution) the other portions of the 
propeller output. This feature can be partially suppressed if desired by the 
input control, PRINTI. This option will suppress the vector printout 
quantities if requested. 

The next output for this initial output is a section entitled: PROGRAM 
INPUT SUMMARY, and consists of the input data displayed in a structured 
format. The propeller modeling options used for the particular execution are 
listed in the following form. The integer value of the input modeling option 
is displayed with a brief description of the model used. The freestream 
conditions are listed next, ( VKTAS , S0UND, DENSTY) followed by the propeller 
operating characteristics (PR0PNM, BLADEN, RPM, ZHUB) . The parameters which 
define the wake and blade geometry segmentation are displayed next (STN, 

STACK, DPSI, REV, CNSECT) followed by the propeller characteristics (RADI, 
HUBQ1 , THETA1 , RDCAS1 , VIM0M1) of each propeller. The printout of the 
propeller characteristics includes: the blade lifting line segment boundaries 

coordinates (XSB, YSB, ZSB, BETA); the inflow station coordinates and the 
blade properties at the segment centers as interpolated from the input distri- 
butions (AIRN, C0RD, THET, DECL, T0VC , V0VO, URVO, S0UN, DENS). 
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Intermediate 


The intermediate output is described below. Generally it is limited to 
the minimum amount possible since most of the output is repeated in the final 
output section. The intermediate output is repeated for each iteration in 
blade angle. If there is no performance iteration it is printed only once for 
the input blade angle. This output is entitled: PROGRAM OUTPUT FOR PROPELLER 

PERFORMANCE ITERATION NUMBER X. 

The first output data for this section is not optional; it consists of a 
table of the blade lifting line segment center and boundary coordinates for 
the reference blade angle and the blade angle value in degrees. The coordi- 
nates are listed in cartesian and cylindrical form for the centers and bound- 
aries. Following this table, the coordinates for the definition of the tip 
Mach cone location are listed in cartesian form for the blade angle in ques- 
tion. If no optional printouts are requested, the wake transport velocity 
distribution is printed as a function of blade radial location. If optional 
printouts are requested this print does not follow immediately, but occurs 
later in the output. All other output for the intermediate portion is 
optional. The input option controls (in parentheses) and the respective 
descriptions of the output follow (output labels in parentheses if not 
noted) . 

The output of the trailing wake geometry coordinates (PR0PT) is presented 
in the cylindrical coordinate system and tabulated as a function of wake 
azimuth position and blade radial position. This output is entitled: WAKE 

COORDINATES. The radial coordinates are tabulated first, starting with the 
values at the blade and ending with the oldest element. A table of the axial 
coordinates is then presented in the same format. 

The printout (PRMAT) of the summed geometric influence coefficients for 
each inflow station at each propeller position for each propeller as a func- 
tion of the appropriate inflow station and propeller position of each 
propeller is presented in the cylindrical coordinate system and in the blade 
element coordinate system. This printout is entitled: CYLINDRICAL GEOMETRIC 

INFLUENCE COEFFICIENTS or BLADE ELEMENT GEOMETRIC INFLUENCE COEFFICIENTS. The 
propeller and propeller position indices are so noted on the printout, while 
the inflow station indices are not, since the values are presented for the 
inboard station to the outboard station for each propeller and propeller 
position. 

For detailed intermediate output (DEBUG) the following extensive list of 
items is presented in the same format as noted above for the propeller, pro- 
peller position and inflow station indices. Generally this printout should 
not be used. A section of detailed blade element properties consisting of the 
magnitude (if it applies) and unit vector direction cosines in the cylindrical 
coordinate system is output for each of the following: the blade segment 
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lifting (SB, ALSRAD, ALSPHI, ALSAXL) , input chord (C1NPUT, ALCIRD, ALCIPH, 
ALCIAX), the normalwise unit vector (ALNRAD, ALNPHI, ALNAXL), the blade 
element chord nondimensionalized by blade radius (CH0RD, ALCRAD, ALCPHI, 
ALCAXL), input blade element thickness to chord ratio (T0VERC, ALTIRD, ALTIPH, 
ALTIAX), the blade element thickness to chord ratio magnitude only (THK) and 
the blade element design lift coefficient (DESCLP). This section is entitled: 
DETAILED BLADE ELEMENT 0UPUT . The next section consists of detailed blade 
element velocities and unit vector direction cosines (VT0T, ALVRAD, ALVPHI, 
ALVAXL) in the cylindrical coordinate system along with the direction cosines 
in the blade element coordinate system (VS, VC, VN) and the angle of attack 
(ALPHAN), and inplane aerodynamic skew angle (SKEW), all computed without 
including the propeller induced velocities, entitled: DETAILED VELOCITY 

RELATED OUTPUT (EXCLUDING INDUCED VELOCITY TERMS). The indices associated 
with the internal program "DO LOOPS" for the geometric influence calculations 
are then output, along with the cosine and sine functions for the respective 
inflow stations inplane lag angle and propeller azimuth position and a counter 
rotation flag with each line of output marked: INTERMEDIATE OUTPUT. 

Following this output, the normalwise blade element geometric influence 
coefficients are printed in the circulation matrix form for each propeller at 
each propeller position for each inflow station. The title of the output is: 
GEOMETRIC INFLUENCE COEFFICIENT. Following this output some untitled cascade 
related items are listed. The chord-to-gap ratios (TAU) and gap-to-chord 
ratios (SIGMA) are listed along with the geometric angle between the propeller 
direction of rotation and the local blade element chordwise vector (THETAG) 
which represents the compliment of the cascade stagger angle. Tip Mach cone 
quantities (untitled) are then output. The Mach cone angle is listed and then 
the angle between the blade tip and the location of the specified fraction of 
the blade chord for each inflow station is listed. The tip Mach number value 
is listed next. The station location index (NSTAT) for the intersection of 
the Mach cone and the fraction of the blade chord is listed and the resulting 
Evaard Tip Relief correction factor (XKC0NE) for each blade inflow station is 
presented. The blade element Mach number (SMACH) and the total Mach number 
( CMACH) distributions are listed along with the blade element geometric angle 
of attack (excluding induced terms) distribution in radians (ALPHA). 

Following this output, the linearized lift curve slope (AA) , an aerodynamic 
quantity (D ■ •S— &.) , the matrix co'nstant vector (C0NST) and the blade element 
geometric blade angle (THETA) in degrees are listed. The induced velocity 
component distributions (VIN, VIC, VIS) in the blade element coordinate system 
at each inflow station of each propeller for each propeller position are 
presented due to each propeller, followed by the total of both propellers for 
each inflow station (VINT, VICT, VIST), entitled: DETAILED INDUCED VELOCITY 

OUTPUT. This output is followed by the total velocity magnitude (VT0T), the 
blade element angle of attack (ALPHA), the blade element aerodynamic skew 
angle (SKEW) and blade element inflow angle (PHI) distributions which include 
the induced velocities. It is titled: DETAILED VELOCITY RELATED OUTPUT. 
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The above output starting from the Mach cone correction and ending with 
the inflow angle is repeated for each iteration of the nonlinear circulation 
matrix solution. The intermediate circulation solution output consists of the 
nonlinear correction quantities (C0RPHI, C0RVEL, C0RCL) used in the solution 
technique, the resulting corrected constant vector of the circulation matrix 
(C0NHSD), the actual correction vector (CFDP), the uncorrected constant vector 
(C0NST), the current angle of attack (ALPHA), and previous angle of attack 
(SAVALP), the current lift coefficient (CLSAV), the current circulation (CIRC) 
and previous circulation (SAVCIR), and the current normalwise induced velocity 
(VIN) for each inflow station for each propeller position of each propeller 
for each iteration of the matrix solution. Once the final circulation 
iteration solution is obtained, the final lift, drag and minimum drag coeffi- 
cients are printed (CLSAV, CDSAV, CDO) . Following this output the total blade 
forces are listed in terms of the magnitude and direction cosines (FT0T, 
ALFRAD, ALFPHI, ALFAXL) and the respective lift and drag components of the 
force (FLT0T, ALFLRD, ALFLPH, ALFLAX, FDT0T, ALFDRD, ALFDPH, ALFDAX). This 
output is marked: DETAILED BLADE FORCE SUMMARY. 

Fina^ Outpjjt 

The final output consists of tabled values of many of the output items 
listed in the intermediate printout and integrated performance quantities. 

This output section is entitled: PROPELLER PERFORMANCE. It is repeated for 

each performance iteration and presented for each propeller for each propeller 
position as a function of blade inflow station location (X/R). The descrip- 
tion of each of the tabulated items is included on the printout of Appendix B 
and will not be described here. It is labeled: BLADE SPANWISE VARYING 
QUANTITIES. Only the descriptions of the sections of integrated quantities 
will be presented. The first of these integrated sections is labeled BLADE 
CHARACTERISTICS and contains the blade characteristics for each propeller 
position for each propeller; thrust per blade (lbj), torque per blade (ft- 
lbf), power per blade (ft-lb£/sec) and horsepower per blade (hp) . 

Following this section of integrated quantities, the combined (all blades, 
both propellers) instantaneous values of thrust and power for each propeller 
position are presented. It is titled: INSTANTANEOUS TOTAL PROPELLER PERFOR- 

MANCE FOR PROPELLER POSITION X. This is followed by a section of integrated 
values averaged over all propeller positions for each propeller, entitled: 
INTEGRATED PROPELLER CHARACTERISTICS FOR PROPELLER X. This section contains 
the total thrust (Ib^), thrust coefficient (T/ n 2 D*M, forward velocity 
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(knots), torque (ft-lbf), power coefficient (P/pn 1 2 3 4 D 5 ), advance ratio 
(Vir/ftR), profile torque (ft-lbf), propeller efficiency (CTXJ/Cp), reference 
blade angle (degrees), induced torque (ft-lb^), power ( f t-lb^/sec) , horse- 
power (hp) and the momentum induced velocity (fps). The combined propeller 
performance follows if coaxial propellers are used. This output is followed 
by the nacelle and combined nacelle-propeller quantities. These items for the 
nacelle are the pressure and skin friction drag (lbf), the respective drag 
coefficients and the combined drag and drag coefficients using the same units 
and definitions as used for the propellers. The combined nacelle and 
propeller thrust, thrust coefficient, power and power coefficient and 
efficiency then follow. Following this output, the force components per blade 
per unit span are presented in the cylindrical coordinate system (lbf/ft) 
for each propeller, and labeled: FORCE PER BLADE PER UNIT SPAN. 

Description of Failure Modes 

Generally, if the input data is correct and reasonable for the flight 
condition being investigated, the propeller solution procedure will not fail. 
To help assist the user in running the computer program, certain failures 
which could occur because of incorrect data setup or incorrect data values are 
checked internally by the computer program. If the input is incorrect, 
diagnostic output will occur to inform the user and allow him to make the 
required corrections. 

Q e IL e L a — i n JP u i- Fojnna t 

As noted in the section describing the input data setup, certain label- 
ling formats have been specified for the input data. If these formats are 
violated, explicit output diagnostics will not generally be printed; however, 
program termination will occur immediately with the last item which was 
attempted to be read in printed as the last output. Termination on the input 
of these labels will occur for the following reasons: 

(1) Data set labels not in the required order 

(2) Data set labels mispunched 

(3) Input item labels mispunched 

(4) Missing END labels for the data sets 

Mi_s sjjng^ _Input_ Data 

Assuming all of the input data is read in correctly, the program then 
checks for missing input that is required for successful program execution. 

The following diagnostic messages could occur if certain data is missing. 
Explanations of the messages are noted, if required, for clarity. 
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(1) "PROPELLER DISK DISPLACEMENT NOT INPUT, EXECUTION TERMINATED" 


This message informs the user that the hub displacement between the propellers 
was not input in the coaxial mode of operation. 

(2) "RPM NOT INPUT, EXECUTION TERMINATED" 

(3) "SOUND NOT INPUT, EXECUTION TERMINATED" 

This message informs the user that the freestream value of the speed of sound 
was not input. 

(4) "DENSITY NOT INPUT, EXECUTION TERMINATED" 

This message informs the user that the freestream value of the density of air 
was not input . 

(5) "RADIUS NOT INPUT, EXECUTION TERMINATED" 

This message informs the user that a blade radius input is missing. 

(6) "DPSI NOT INPUT, EXECUTION TERMINATED" 

(7) "NUMBER OF WAKE REVOLUTIONS NOT INPUT, EXECUTION TERMINATED" 

(8) "NUMBER OF BLADES NOT INPUT, EXECUTION TERMINATED" 

In£°rrect_ I)a£a_Inp_ut 

If data is input to the program which is incompatible with the require- 
ments of the computer analysis, diagnostic messages will also occur. The 
messages are listed below along with explanation, if required. 

(1) "POWER COEFFICIENT ITERATION NOT ALLOWED FOR TWO PROPELLERS, EXECU- 
TION TERMINATED" 

(2) "THRUST COEFFICIENT ITERATION NOT ALLOWED FOR TWO PROPELLERS, EXECU- 
TION TERMINATED" 

(3) "COMPRESSIBLE BOUND VORTEX MODEL NOT FUNCTIONAL FOR TWO PROPELLERS, 
EXECUTION TERMINATED" 

This message informs the user that he has requested a combination of modeling 
options which are not compatible. The compressible bound vortex model was not 
derived for coaxial propellers, and thus cannot be used for coaxial propeller 
configurations . 
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(4) "INPUT ERROR 360/DPSI IS NOT A MULTIPLE OF B. WILL STOP PROGRAM. 
JTOT=X, B=X" 

This message informs the user that the requested blade azimuth increment is 
not an integer multiple of 360 degrees. The number of blade azimuth positions 
(JTOT) and the number of blades (B) that were requested are listed in the 
locations marked by X respectively. 

( 5 ) "***BJT0T IS NOT AN INTEGER MULTIPLE OF THE NUMBER OF PROPELLER 
DISKS, EXECUTION TERMINATED" 

This message informs the user that the number of azimuth intervals between 
blades is not an integer multiple of the number of propellers. It checks to 
be sure that for a coaxial configuration, the half blade spacing is an integer 
multiple of the azimuth increment. 

There are also a series of diagnostic messages associated with internal 
program core allocations. If a combination of input quantities exceeds the 
internal dimension limits, self-explanatory messages are output which inform 
the user of the problem, the values input and the allowable limits. Because 
the messages are self-explanatory, they will not be listed here. The required 
corrective action will be clear to the user if they do occur. 
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Nacelle Portion 


This section is intended to describe the general features of the nacelle 
portion of the PANPER program. The technical aspects of this analysis are 
described in reference 1. The first subsection describes what problems can be 
solved and what problems cannot be solved. It also describes any special care 
which should be used in exercising the various options. The second and third 
subsections present a detailed description of the input which is required in 
the operation of the computer program and the interpretation of the printed 
output. Since any complicated computer program may fail due to inconsisten- 
cies in the input or failure of the theory, the computer program is provided 
with self-diagnostics which notify the user of the type of failure. The last 
subsection deals with these program diagnostics as well as helpful hints to 
correct problems which may be encountered. 

Since this computer program is intended for a wide variety of users, some 
note should be made of the nomenclature. The term "duct" refers to any flow 
passage including inlet nozzles, diffusers, or transition ducts or external 
flow problems where the outer wall is replaced with the appropriate boundary 
condition. Typically, such ducts may have struts, compressor or propeller 
blades, inlet guide vanes, or exit guide vanes and these terms are used almost 
interchangeably in the discussion. Depending on the user, the duct wall 
dimensions may be referred to as hub and tip walls or inside diameter (ID) and 
outside diameter (OD) walls respectively. Some users may use the terms 
centerbody and outerbody when referring to ID and OD walls respectively. The 
subscript notation, Fortran symbols, and computer printout generally uses the 
subscript W for either wall without distinction and H and T for hub and tip 
wall. Finally, the term "slot injection" refers to the injection of flow 
tangent to the wall at a discrete axial location, while "mass bleed" refers to 
injection of flow normal to the wall. 
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General Features of the Program 


T^ypes_ o^_FJ_uJ_ds_ 

The fluid may be any compressible gas as defined by its thermodynamic 
properties/?, Cp, C y , p, P RL , p rt* ^ not ot ^ erw ise specified, the gas is 
assumed to be air. The reference conditions for the gas properties must be 
specified at standard sea-level conditions. 

T’ypes^ o^_Fj_ow ^i^ujit_ion_s 

External or internal, transonic, turbulent, swirling or nonswirling flows 
may be calculated, including flows with radial total pressure distortion. 
Two-dimensional flows may be calculated by constructing an annular duct in 
which the inner to outer radius approaches 1.0. 

t£_y_Op t J^on PT 3_) 

The flow through any axisymraetr ic duct may be calculated provided that 
the flow is generally in the axial direction. Duct flows normal to the axis 
of symmetry or which reverse direction cannot be calculated due to logic 
limitations in Subroutine C00R. Ducts with sharp discontinuities, such as a 
step, which produce separation also cannot be calculated. 

Provision is made in the program to either read the duct coordinates from 
input data cards (I0PT3=2), or to calculate the duct coordinates analytically 
(I0PT3>4) from a few input duct shape parameters. If the duct coordinates are 
read from input cards, care should be taken that the input coordinates have 
sufficient smoothness to calculate the first and second derivatives using 
numerical finite-difference equations. When the second option is used 
(I0PT3>4), the user must program his own calculation in Subroutine GDUCT. 
Sample programs (I0PT=1, 3, 4) are given in Subroutine GDUCT for the user’s 
reference. For ducts with no centerbody a zero radius must be specified. 

An important restriction to the computer program is that the inlet and 
exit flow must have no normal pressure gradients produced by streamline curva- 
ture, although it may have normal pressure gradients due to swirl. Many ducts 
do not satisfy this requirement; however, these ducts can still be treated if 
the duct is extended. For curved annular ducts exhausting to atmosphere, the 
exit flow may have curvature. This phenomena may be simulated by extending 
the duct to approximate the curvature of the exit flow. 
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/ 


If the I0PT3=2 option is used, and the number of input points is less 
than the number of specified streamwise stations, the program smooths the 
input data and interpolates the required mesh points, 

_Inlet_ ^i£ w _°£tJ_O£s_(^I0P^TJ_) 

The computer program is provided with two methods to describe the inlet 
flow. When I0PT1=1, the inlet flow is calculated by prescribing the stagnation 
conditions (P 0 ,T Q ) on Card No, 6, the inlet Mach number M , the swirl 
angle aj, , and the boundary layer parameters 6 and n, which are the boundary 
layer displacement thickness and power law velocity profile exponent, on Card 
No. 5, respectively. The core flow is then calculated from isentropic flow 
relations, and boundary layers added using power law velocity profile rela- 
tions. When stagnation conditions are not specified, the calculation assumes 
sea level conditions. 

When I0PT1=2, the inlet flow is prescribed from input data cards which 
specify the stagnation pressure P Q , static pressure P, swirl angle a, and 
stagnation temperature T Q , as a function of the fractional distance across 
the inlet. This data need not be specified at equidistant points since a 
linear interpolation is used to specify the data at the mesh points used in 
the calculation. If experimental data is not used, care should be taken that 
the data is self-consistent and that it satisfies the radial equilibrium equa- 
tion. Since the initial growth of the boundary layer is sensitive to the wall 
shear stress, data describing the boundary layers should be accurately 
specified. When this is not possible, boundary layers may be added to each 
wall by specifying 6 and n. Special care should be exercised in using the 
I0PT1=2 option, with or without the feature of adding in the wall boundary 
layers. If the stress distribution across the duct is not smooth and 
realistic, numerical instabilities might originate in the inlet flow and grow 
rapidly to a point where the calculation is terminated. This may take the 
form of an unrealistically early separation. 

When I0PT1-3, the inlet free stream flow is calculated the same as 
I0PT1*1. The boundary layers on each wall, however, are calculated from 
Coles 1 profiles (reference 3) using Function FC0LES. The I0PT1-4 option is 
the same as the I0PT1*2 option, except that Coles 1 profiles are used for the 
boundary layers. 

For I0PT1=1, 2, 3, or 4 there are no restrictions on 6* other than it 
must be greater than zero and that the transverse grid must be chosen such 
that at least 5 to 10 mesh points exist for 0<^ Y + _< 10. A printout of U + (Y + ) 
is provided by setting I0PT4=O. In absence of other information a value of 
6 of one percent of the inlet height is an adequate approximation for a 
thin initial boundary layer. If the boundary layer thickness is not small 
compared to half-height, the correct input value of 6 must be obtained 
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from other sources such as data correlation, experimental measurements, etc. 
Most zero pressure gradient boundary layers follow a l/7th power law profile 
and it is recommended that this value be used. For I0PT1=3 or 4 in which 
Coles' profiles are used, a shape factor is computed from the input values 
is used 6 and n. This shape factor is used to compute a wake parameter and 
a compatible wall stress for use in Coles' profiles. As shown in reference 3, 
specification of the wake parameter and wall stress uniquely defines the 
Coles' velocity profile. 

Boundary Conditions (T w , m w ) 

Either the adiabatic wall or the heat transfer case may be calculated. 

The program assumes adiabatic walls unless the wall temperature is specified. 
Any wall temperature distribution may be specified, either on input cards when 
the duct coordinates are read, or calculated when the duct coordinates are 
calculated. The case of wall bleed may also be treated in a similar manner; 
wall bleed flow rate is zero, unless otherwise specified. At the present 
state of development of the computer code, only the I0PT3=1 option allows a 
specification of wall temperature as a boundary condition. For all other 
I0PT3 options adiabatic walls are assumed. 

Force Option 

Subroutine F0RCE is provided with two options. For I0PT2 * 0 and N0PPF 
= 1, the blade force is calculated from data taken from the propeller lifting 
line portion of the code. For I0PT2 * 0 and N0PPF = 0 the blade forces are 
read in as input data. 

Failure Modes 


In the event of failure in the calculation, the program prints an error 
message called "diagnostic". These "diagnostics" are in addition to the 
computer diagnostics and are clearly labeled as such. These "diagnostics" 
terminate the calculation only when very serious. A list of these "diagnos- 
tics" appears in a later section. Included with this list is an identifying 
number for the "diagnostic", the location (Subroutine), and the immediate 
cause of the failure. Where possible, suggestions are made to correct the 
calculation. 

Debug, Options UDBGN2 

Auxiliary printout which was originally used to debug the computer 
program is available to the user by setting the appropriate IDBGN option. 
However, the user must refer to the program listing or compilation to deter- 
mine the meaning of this printout. 
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Grid Selection 


The grid selection parameters appear on the third input card and are 
given by DDS, KL, JL, KDS. The number of streamwise stations is divided into 

a coarse grid of JL ^ 100. The number of streamlines including the wall 

boundaries is given by KL _< 100 points and a fine grid of JL*KDS points. The 

solution is numerically stable; however, truncation errors may get large if 

too large a streamwise step size is used. The streamwise step size may be 
made smaller without recalculating the coordinate system by increasing KDS. 

It should be noted that computing time is proportional to JL*KDS. The para- 
meter DDS distorts the normal coordinate by placing more streamlines near the 
wall . 


Mesh Distortion 


The numerical solution of turbulent boundary layers requires accurate 
integration of the mean profile in the turbulent mixing layer. For high 
Reynolds number flows, practical considerations require distributing more mesh 
points near the wall in some systematic manner. This is done using an 
exponential transformation given by 


n( n) 


(c+l/2)exp 

2 In 

(m) 

(rrl/2) - (c-1/2) 

1+exp 

21n 

(m) 

(nrl/2) 


( 1 ) 


where 


0 i n i 1 

0 < n < 1 


The parameter c is chosen so as to place the first mesh point at approximately 
Y* ■ 1. Then for equal increments in An> equation (1) distributes the mesh 
points tn so as to place more mesh points near the wall. 

'Sej>ajration 

The separation point is determined when the streamwise component of wall 
stress goes to zero. However, the calculation can continue past the separa- 
tion point. When the region of reverse flow becomes too large, greater than 
2.0 percent, the calculation stops. 
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Description of Input 


This subsection describes the loading of input data cards for running the 
nacelle portion of the computer program. The input specification follows the 
convention that a blank or zero value for any parameter implies no action by 
the computer program. Numbered cards must be loaded. The remaining cards 
must be loaded only if the proper option is selected. Care should be taken in 
loading the program because of the input changes depending on the options 
chosen in the second data card. Multiple cases can be run simply by stacking 
the cases in order. The last case is followed by two blank cards. 


Card No. 

1: Title 

Card 


Name 

Col. 

Format 

Comments : 

TITLE 

1-72 

12A6* 

Any alphanumeric 

Card No. 

2: Option Card 


Name 

Col. 

Format 

Comments : 

I0PT1 

1-2 

12 

(FL0W1N Option) 


I0PT1®3 The inlet flow is computed by 
specifying the data on card 5. 

I0PT1=4 The inlet flow is read from 2xKLL data 
cards following card 5. 

I0PT1*9 Laminar flow, inlet flow is calculated 
using a Blasius profile. 


I0PT2 3-4 12 (F0RCE Option) 

I0PT2=O No blades or struts exist in the duct 
and these cards are not loaded. 

I0PT2*3 The strut forces are input on cards 
(2xKLL cards following card 3). 

I0PT2*4 The blade forces are calculated from 
lifting line theory if N0PPF=1 . If N0PPF=O 
blade forces are read from data cards. 


I0PT3 5-6 12 (GDUCT option) 

I0PT3=1 Calculate a straight annular duct 
I0PT3=2 Read coordinates 

I0PT3=3 Calculate a straight wall annular 
diffuser 

10PT3=4 Do not use 

I0PT3*5 Calculate curved wall diffuser No. 1. 
I0PT3=8 Straight walled duct 


12A6 UNIVAC system 
18A4 IBM systems 
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Card No. 2: Option Card (Cont'd) 


Name 

Col. 

Format 

I0PT4 

7-8 

12 


I0PT5 9-10 12 


I0PT6 11-12 12 

I0PT7 13-14 12 

I0PT8 15-16 12 

I0PT9 17-18 12 


I0PT1O 19-20 12 

I0PT11 21-22 12 


Comments : 

Print solution every I0PT4 stations. For 
example, if I0PT4=3 every third station is 
printed. If I0PT4=1 every station is printed. 
If I0PT4=-1 additional output at each 
station is printed. 

Strut data input (see I0PT2=3) used to 
calculate strut forces from experimental data 
measured upstream and downstream of strut. 
I0PT5=1 Read in required profiles. 

I0PT5=2 The upstream and downstream strut data 
cards are identical to the inlet and exit flow 
cards and need not be loaded. 

I0PT6=O Strut force plus thickness effects. 
I0PT6=1 Strut thickness effects only. 

Axisymmetric compressible streamline curvature 
corrections . 

0 = No curvature correction 

1 = Curvature correction 


WBLEED option. 

= 0 No Bleed 
= 1 Bleed OD wall 
= 2 Bleed ID wall 
= 3 Bleed OD and ID wall 


I0PT9=O Approximate C00R calculation. 

I0PT9=1 Exact C00R calculation. 

I0PT9=2 Store C00R calculation on mass storage 
device (Unit 9) and stop. 

I0PT9=3 Read C00R calculation from mass 
storage device (Unit 9) and stop. 

For normal running, set I0PT9=1. Subroutine 
C00R is described in a later section. 


I0PT1O= 1 Internal flow problem. 
I0PT1O=O External flow problem. 


I0PT11=1 External flow problem. 
I0PT1 1=0 Internal flow problem. 
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Name 

Col. 

Format 

Comments: 

I0PT12 

23-24 

12 

Not used. 

I0PT13 

25-26 

12 

Not used. 

I0PT14 

27-28 

12 

Not used. 

I0PT15 

29-30 

12 

Start flow calculation at station I0PT15. 
Default =1.0. 

I0PT16 

31-32 

12 

End flow calculation at station I0PT16. 
Default = JL. 

I0PT17 

33-34 

12 

Not used. 

Card No. 

3 : Mesh 

Parameters 


DDS 

1-10 

F10.3 

Mesh distortion parameter, default determined 
internally. 

KL 

11-13 

13 

Number of streamlines including wall, 2 < KL < 
100. 

JL 

14-16 

13 

Number of streamwise stations, JL < 100. 

KDS 

17-19 

13 

Number of steps per streamwise station. 
Default = 2. 

KLL 

20-22 

13 

Number of streamlines of data input (see 
I0PT1, I0PT2). If KLL < KL, inlet flow is 
interpolated from KLL inlet data cards on the 
KL streamlines used for calculating flow. 

KLL £ 31. 

JLAST 

23-25 

13 

Number of C00R records stored on drum. Used 
for tape storage of coordinate functions. Not 
used in this version. 

JLPTS 

26-28 

13 

Number of input duct coordinate points, if 
I0PT3 = 2. Note: If JLPTS < JL, points are 

smoothed and interpolated. Not used in this 
version. 

LFILE 

29-31 

13 

Case stored on tape file LFILE used for tape 


storage of coordinate functions. Not used in 
this version. 
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Card No. 4: GDUCT 


Name Col. Format Comments: 

These input cards are read in subroutine GDUCT 
as programmed by the user. 

The following duct geometries (designated as I0PT3=1 ,2 ,3 , and 5) have been 
programmed (see figure 6). 

Card No. 4: (10PT3=1) Straight Annular Duct 


Name 

Col. 

Format 

Comments : 

Z1 

1-10 

F10.0 

Length (ft) 

RH1 

11-20 

F10.0 

Centerbody radius (ft) 

RT1 

21-30 

F10.0 

Outerbody radius (ft) 

TWH 

31-40 

F10.0 

Centerbody wall temperature (deg R) 

TWT 

41-50 

F10.0 

Outerbody wall temperature (deg R) 

AMWH 

51-60 

F10.0 

Centerbody wall bleed ( lb/ ft 2 sec) 

AMWT 

61-70 

F10.0 

Outerbody wall bleed (lb/ft 2 sec) 

Card No. 

4: (I0PT3“2) Arbitrary Duct Input 

Name 

Col. 

Format 

Comments : 

Z1 

1-10 

F10.0 

Duct length (ft) 

RNOPE 

11-20 

F10.0 

One more than the number of curve fits used 
for smoothing the input geometry. Default 5J 

XNOSE 

21-30 

F10.1 

Distance to nacelle nose (ft) 

Cards Following 

Card No. 4: 

(I0PT3*2) For JLPTS equally spaced points, thus 




Z( J)-Z1 *( J-l ) / ( JLPTS-1 ) 

Name 

Col. 

Format 

Comments: 

R(1,1,J) 

1-80 

8F10.0 

Outerbody radius (ft) 

R(2 ,1 ,J) 

1-80 

8F10.0 

Centerbody radius (ft) 
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Card No. 4: (I0PT3=3) Straight Wall Annular Diffuser 


Name 

Col. 

Format 

Comments: 

Z1 

1-10 

F10.0 

Duct length (ft) 

RH1 

11-20 

F10.0 

Centerbody radius (ft) 

RT1 

21-30 

F10.0 

Outerbody radius (ft) 

ZTHRO 

31-40 

F10.0 

Length of throat (ft) 

ANGH 

41-50 

F10.0 

Centerbody wall angle (deg) 

ANGT 

51-60 

F10.0 

Outerbody wall angle (deg) 

Card No. 

4: (I0PT3=5) Curved Wall Annular Diffuser No. 1 

Name 

Col. 

Format 

Comments : 

Z1 

1-10 

F10.0 

Duct length (ft) 

RT1 

11-20 

F10.0 

Inlet outerbody radius (ft) 

RH1 

21-30 

F10.0 

Inlet centerbody radius (ft) 

RTL 

41-50 

F10.0 

Exit centerbody radius (ft) 

RHL 

51-60 

F10.0 

Exit outerbody radius (ft) 

AT 

61-70 

F10.0 

Power outerbody wall AT 2 

AH 

71-80 

F10.0 

Power centerbody wall AH_2 

Card No. 

4: (I0PT3-6) Not 

Used 


Card No. 5: Inlet Flow Distribution (See figure 7) 


Name 

Col. 

Format 

Comments : 

AMS1 

1-10 

F10.0 

Nominal inlet Mach number 

ALP1 

11-20 

F10.0 

Nominal swirl angle at hub (deg to axis) 

DSH 

21-30 

F10.0 

Boundary layer displacement thickness on hub 
wall (ft) 


36 


Name 

Col . 

Format 

Comments : 

DST 31- 

-40 

F10.0 

Boundary layer displacement thickness on 
tip wall (ft) 

ANH 41- 

-50 

F10.0 

Power law exponent for hub boundary layer 

ANT 51- 

-60 

F10.0 

Power law exponent for tip boundary layer 
For boundary layers of approximately 10 
percent of inlet height, nominal values 
for DSH and DST are 0.0125 times inlet 
height and ANH, ANT equal to 7.0. 

2xKLL Inlet 

Flow Cards Following Card 5 (Only if 10PT1=4) 

Name 

Col. 

Format 

Comments : 

B1NPUT(1 ,K) 

1-10 

F10.0 

Fractional distance across duct Y 

BINPUT(2,K) 

11-J20 

F10.0 

Total pressure (lb/ft 2 abs) P Q 

B1NPUT(3,K) 

21-30 

F10.0 

Static pressure (lb/ft 2 abs) P 

BINPUT(4,K) 

31-40 

F10.0 

Swirl angle to axis (deg) a 

BINPUT(5,K) 

41-50 

F10.0 

Total temperature (deg R) T Q 


The first KLL cards describe the inlet 
flow. The second KLL cards describe the 
exit flow. If the exit flow is not known, 

ic 

KLL blank data cards must be used. If 6«>0 
on Card 5, boundary layers are added 
according to Card 5. 

Blade Row Data Card Following Card 5 (Only if 10PT*O) 

This sequence of data cards is repeated for each blade row. This data must be 
consistent with the propeller input geometry. 


Name 

Col. 

Format 

ZCLI 

1-10 

F10.0 

NBLADE 

11-13 

13 

1 SHAPE 

14-16 

13 


Comments: 

Axial location of blade centerline 
Number of blades 
Blade shape index 
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Name 

Col. 

Format 

Comments: 

NUM 

17-19 

13 

Number of points defining blade segment 
boundaries 

0MEGZ1 

20-29 

F10.0 

Rotational velocity (rpm) 

LROW 

30-32 

13 

Blade row counter 

NROW 

33-35 

13 

Number of blade rows 


Blade Row Geometry Cards 


The blade row geometry cards are read in if the nacelle portion of the 
program is operating without the propeller portion. 

NUM times the blade row geometry cards (root to tip) noted below, follow each 
blade row data card (10PT2*O, N0PPF=O) 


Name 

Col. 

Format 

Comments : 

C0NST1 ( 1 ,K) 

1-10 

F10.0 

Blade 

radius (ft) 

C0NST1 (2 ,K) 

11-20 

F10.0 

Chord 

angle to blade face (deg) 

C0NST1 (3 ,K) 

21-30 

F10.0 

Chord 

(ft) 

CONSTl (4 ,K) 

31-40 

F10.0 

Thickness /Chord 

C0NST1 (6 ,K) 

41-50 

F10.0 

Axial 

location at blade quarter chord (ft) 


Arbitrary Strut Thickness Distribution (XSHAPE * 4) 


Name 

Col. 

Format 

Comments 

KBLADE 

1-10 

110 

No chordwise stations (KBLADE < 50) 


Chordwise location 


Name 

Col. 

Format 

Comments : 

X(K) 

1-80 

8F10.0 

Chordwise location 

Y(K) 

1-80 

8F10.0 

Thickness/maximum thickness distribution 
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2xKLL Strut Data Cards Following (10PT5 f 0, I0PT2 =3) 


Name 

Col. 

Format 

Comments: 

A1NPUT(1 ,K) 

1-10 

F10.0 

Fractional distance across duct Y 

A1NPUT(2,K) 

11-20 

F10.0 

Stagnation pressure (lb/ft 2 abs) P Q 

A1NPUT(3,K) 

21-30 

F10.0 

Static pressure (lb/ft 2 abs) P 

A1NPUT(4,K) 

31-40 

F10.0 

Swirl angle to axis (deg) a 

A1NPUT(5,K) 

41-50 

F10.0 

Stagnation temperature (deg R) T Q 


The first KLL cards describe the inlet 
flow of the strut row. The second KLL 
cards describe the exit flow (see Card 2). 


Card No. 6: Performance Point 

If this card is left blank, the default values shown in parentheses are 


used. 





Name 

Col. 

Format 

Comments : 


PRESO 

1-10 

F10.0 

Inlet stagnation pressure (2117. 
abs) 

lb/ft 2 

TEMPO 

11-20 

F10.0 

Inlet stagnation temperature (519 

. deg R) 

AC1 

21-26 

F6.0 

Clauser constant (0.016) 


AK1 

27-32 

F6.0 

Von Karman constant (0.41) 


API 

33-38 

F6.0 

Van Driest constant (26.0) 


PRT1 

39-44 

F6.0 

Turbulent Prandtl number (0.8) 


PRLI 

45-50 

F6.0 

Laminar Prandtl number (0.9) 


CPR 

51-60 

F10.0 

Specific heat at constant pressure (5997 
ft 2 /sec 2 ) 

CVR 

61-70 

F10.0 

Specific heat at constant volume 
ft 2 /sec 2 ) 

(4283 
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Name 

Col. 

Format 

Comments : 

VISCR 

71-80 

F10.0 

Viscosity (0.37E-06 lb/sec ft^) at 
Standard Conditions. 

Card No, 

7: Wall Bleed Card 


CDISH 

1-10 

F10.0 

Discharge coefficient for holes 
(dimensionless) CDISH <1.0 

AHAS 

11-20 

F10.0 

Ratio of hole area to surface area 

Up 

21-30 

F10.0 

Plenum total temperature (°F) 

pip 

31-40 

F10.0 

Plenum total pressure (psta) 

XBF 

41-50 

F10.0 

Wall distance - start wall bleed (ft) 

XBL 

51-60 

F10.0 

Wall distance - end wall bleed (ft) 

Blade Force Card 



NUM blade 

force cards 

following Card 6 (only if I0PT2*O and N0PPF*O) . Repeat 

for each 

blade row. 



Name 

Col. 

Format 

Comments : 

FRCI(N, 1) 

1-10 

F10.0 

Radial force/span 

FRCI(N,2) 

11-20 

F10.0 

Tangential force/ span 

FRCI(N,3) 

21-30 

F10.0 

Axial force/span 


N=1 , NUM 


40 


Description of Output 


Title Page 

The output presented on this page is self-explanatory except for the 
following variables 


r T 

m l “ / «B P U s dr 

r H 

r T 

a l = / g B dr 

r H 

r T 

p, » — / g B p Vi dr = PRES1 

m l r 

1 r H 
r T 

*1 B — f g B P U 8 (1/2 p U l )dr = DYNP1 

m ] r 
1 r H 


r T 

M. = — / g B p ^Mjdr = MACH1 

m 1 r 

1 r H 

WFL “ 32.2 Ngmj 
USR - m l /p r /a 1 - U r 
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Wall Conditions Page 


REY - r r U r /p r /u r 


This page presents a table of Z, r T , m T> T T , r H , m^, which was calcu- 
lated in Subroutine GDUCT. 

Blade Geometry Page (If I0PT2 * 0) 

This page presents a table of blade geometry properties at each discrete 
point of the lifting line. These properties are radius r CL , stagger angle ctg L , 
chord thickness T/C, and axial location Z q L . Also, the number of blades per 

row and the blade shape are printed. 

Gap Average Inviscid Flow Page 

This page presents the solution for the inviscid flow variables across the 
duct at selected stations depending on I0PT4 . A table of values for Y, P Q , 

P, a > T. # T, M, U, U , U , U , speed of sound a, and P are given. Also printed 

O S Z L 

is the pressure coefficient at the inner wall where 



where * freestream Mach number. 


Inviscid Nacelle Drag Page 

This page prints the nacelle pressure drag calculation of the nacelle. 

Nacelle Wake Corrections Page (If I0PT2 * 0) 

This section indicates the location of the propeller lifting line in the 

(s, n, <f>) coordinate system. The flow conditions of the lifting line are also 

presented. These problems are U r , U 2 , p^, a 2 , sin(T), and cos(T) where 

L Li Ij 

T = 0. 
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Gap Average Viscous Flow Page 


This page presents the solution for the flow variables across the duct at 
selected streamwise stations depending on I0PT4. A table of values for Y, 
u s , U + , a, n, 0, M, H 0 , 0 O , c p are given where 

C p - (n - n(o,o))/q 1 

o 

In addition, the wall values for Z, m, C^, Cf g , Q are printed, where 
and C£ g are defined by 


C f* 


T 


n<t>^l 


c fs “ T ns^l 


The one-dimensional characteristics of the flow are also given: area 
ratio A/A^ , Mach number (isentropic flow) Mj , incompressible and compressible 
flow pressure coefficient CPINC and CPC0MP. 

Wall Surface Conditions Page 

This output page presents a summary of the wall conditions along the 
length of the duct* This xncludes Zjj, Cpjj, ^FH* ^WH* ^SH* ^x* ^px* ^px * 

T wr> A SX> ^st * 

Wall Radiation Summary Page 

This output page presents a summary of information which is useful in 
computing radiation effects. The wall temperature Tyjj, T^, on a differential 
area dA^, dA T , located at point (Z H , r^) , (Z T , r T ) with the sin (wall angle) 
is given. 

Viscous Nacelle Drag Page 

This page prints the nacelle pressure drag, pressure drag coefficient, 
friction drag, and friction drag coefficient of the nacelle. 
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IDBG0 Pages 


Intermediate printouts which were 
by setting the debug options IDBGGl^l . 
input card. The user should refer to 
determine the printout variables. 

used to debug the program may be called 
IDBG0 may be specified on the option 
the program listing in each subroutine to 

Debug Print Out 


Subroutine - Purpose 

IDBG1 


TURB - Debug 

IDBG2 


FC0RCT - Debug 

IDBG3 


FL0WIN - Debug 

IDBG4 


Not Used 

IDBG5 


S0LVI - Debug 

IDBG6 


C00R - Debug 

IDBG7 


F0RCE - Debug 

IDBG8 


MINVRT - Debug 

IDBG9 


SM00TH - Debug 

IDBG10 


GDUCT - Debug 

IDBG11 


Not Used 

IDBG12 


S0LVI - Debug #2 

IDBG13 


CKINPT - Debug 

IDBG14 


Set Number of Streamlines 

IDBG15 


Automatic Step Size Debug, Number of 
Streamlines Calculated (default = 25) 

IDBG16 


Suppress Freestream Instability 

IDBG17 


Not Used 

IDBG18 


GE0MCL - Debug 

IDBG19 


WAKC0R - Debug 

IDBG20 


PERFNA and PERFN2 - Debug 
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Description of Failure Modes 


The nacelle portion of the computer program can diagnose the cause of 
certain failure modes for this portion of the analysis and a printed message 
of the following form is given. 

♦♦DIAGNOSTIC NO. XX FOR ANNULAR DIFFUSER DECK** 

The number XX identifies the type of failure from the list below. 

1) I0PT3 OUTSIDE RANGE OF ALLOWABLE DUCT OPTIONS 

This failure occurs in Subroutine ALTMN. The input option must be 
between 1 _< I0PT3 _< 6. 

2) No solution exists in AMFOR 

This failure occurs in Subroutine AMFOR. This subroutine solves the Mach 
number function 

N * M |l + j 1/2 / (1 + >M 2 ) 

for M given N. This function has a maximum at M * 1. Hence 

N(l) - [2(1 + y)]" 1/2 

Solutions do not exist for values of N > N(l). 

3) MASS FLOW EXCEEDS THE MAXIMUM MASS FLOW POSSIBLE 

This failure occurs in Subroutine AMINLT which solves the Mach number 
function 


N 



Y* .l— 

2(y-l) 


for M given N. This function has a maximum for M * 1 given by 


N(l) 


- 



2(r-D 


corresponding to choked flow. 
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4) Not Used 

5) FOR BEST RESULTS ADD A STRAIGHT ANNULAR CHANNEL INLET 

This diagnostic occurs in Subroutine C00R1. In the construction of the 
duct coordinates, it is assumed that the inlet has no curvature as shown 
in figure 8. This is not a fatal error because small inlet curvatures 
may be tolerated. In order to avoid problems, the best procedure is to 
add a straight annular section to the inlet as shown by the dotted line 
in figure 8. 

6) PROGRAM ASSUMES INLET FLOW HAS CURVATURE 

This diagnostic occurs in Subroutine C00R1. Same as diagnostic 5. 

7) WALL CURVATURE IS TOO LARGE AT STATION X. 

This diagnostic occurs in Subroutine C00R1 usually with bad input data 
describing the duct contour resulting in a numerically discontinuous 
change in wall curvature shown in figure 2. 

8) Not Used 

9) GREATER THAN 1. PERCENT NORMAL PRESSURE GRADIENT ERROR RECALCULATE STATIC 
PRESSURE 

This diagnostic occurs in Subroutine ERPIN. This subroutine integrates 
the radial equilibrium equation 


P T - p H - 



1 


I 


-U2 J2 L u 2 + _q _3R y 2 
V 3n 8 R 3n * 


dj) 

XY 


and compares (P T ~ P H ) to that computed for the input inlet flow (P^-P^) j . 
If the error given by 


E 


l . P T . P H 

(P T ~ p H h 


is greater than 0.01, the input initial static pressure distribution is 
replaced by the above pressure equation and the flow is recalculated. 


10) Not Used 
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11) MASS FLOW REQUIRED EXCEEDS MAXIMUM MASS FLOW POSSIBLE 

This diagnostic occurs in Subroutine CKINPT. Choked flow may exist in 
the duct, and this diagnostic will be printed out. The weight flow must 
be reduced. 

12) PRESSURE RISE EXCEEDS PERMISSIBLE PRESSURE RISE 

This diagnostic occurs in Subroutine CKINPT. This error occurs with the 
failure by the deck to properly set up flow entering duct. Check input 
for errors . 

13) Not Used 

14) BOUNDARY LAYER TOO THIN FOR MESH SPACING 

This diagnostic occurs in Subroutine FL0WIN. The viscous flow calcula- 
tion requires a finite initial boundary layer thickness. In addition, it 
requires enough mesh points to describe the inlet boundary velocity 
profile. The deck assumes arbitrarily that at least five mesh points are 
required. Thus, if this diagnostic occurs, increase the number of mesh 
points, KL, increase the mesh distortion parameter, DDS, or increase the 
assumed inlet boundary layer thickness. Setting DDS * 0 automatically 
sets the mesh distortion parameter for turbulent flow. 

15) TOTAL PRESSURE IS LESS THAN STATIC PRESSURE 

This diagnostic occurs in Subroutine FL0WIN. A check is made on the 
input data for I0PT1 ■ 4, to be sure that P^ > P. 

16) INPUT DATA NOT IN RADIAL EQUILIBRIUM CORRECTIONS APPLIED TO STATIC 
PRESSURE 

This diagnostic occurs in Subroutine FL0WIN. A check is made of the 
input data for I0PT1 ■ 4. If the data is not in radial equilibrium, it 
is assumed that the static pressure is in error, and the other inlet data 
is correct. Then the static pressure is computed from 


.411 = 2 — 

dn T 1 


=± M co8 2 a _ 1_ J& 

XV 3n XR 3n 



1/2 


with the ID wall 


static pressure as a boundary condition. 
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17) INPUT DDS MUST BE SPECIFIED 

This diagnostic occurs in Subroutine FN0RM. At this time there is no 
algorithm to automatically select the mesh distortion parameter DDS for 
laminar flow. 

18) BLADE DATA ERROR IN CKINPT ROUTINE 

This diagnostic occurs in Subroutine CKINPT. Blade data input is 
incorrect. It must be rearranged with increasing Y. 

19) NO UNIQUE SOLUTION FROM MINVRT 

This diagnostic occurs in Subroutine MINVRT. If the matrix set up to 
solve the turbulent flow solution is singular, no solution can be 
obtained. This may occur from numerical truncation error problems. 

20) LEADING OR TRAILING EDGE INDEX OF STRUT OUT OF RANGE 

This diagnostic occurs in Subroutine SLETE. In order to compute blade 
forces, the strut must be wholly contained within the duct length. This 
problem may be eliminated by extending the duct length in a realistic 
manner as shown in figure 10 . 

21) Not Used 

22) Not Used 

23) BOUNDARY LAYER OVERLAP OR TOO LARGE 

This diagnostic occurs in Subroutine FL0WIN. For internal flow, the sum 
of the two boundary layer thicknesses must be less than the duct inlet 
height. Check input data. 

24) SET TOTAL TEMPERATURE, PRESSURE; ANGLE TO VALUE AT EDGE OF BOUNDARY LAYER 
- CORRECTIONS APPLIED 

This diagnostic occurs in Subroutine FL0WIN. For I0PT1 “ 4, calculated 
boundary layer profiles are matched to experimentally measured inlet 
flow. Good matching occurs only if the inlet flow data shows constant 
Pj in the boundary layer region as shown in figure 11 by the dotted 
line. 
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25) TRUNCATION ERROR CANNOT BE REDUCED BY STEP SIZE 

This diagnostic occurs in Subroutine S0LVI, When the step size KDS is 
not specified, it is automatically selected by checking the truncation 
error at each step. When an instability occurs, the program attempts to 
reduce the truncation error by reducing the streamwise step. .If the 
truncation error cannot be reduced below a minimum value, the calculation 
stops with this error message. 

26) NUMERICAL INSTABILITY 

This failure occurs in Subroutine FC0RCT. Temperature and pressure are 
checked for negative values. Calculation stops with this error. 

27) Not Used 

28) Not Used 

29) SOLUTION REQUIRES REVERSE FLOW, INCREASE WFLOW 

This diagnostic occurs in Subroutine CKINPT. For flows with radial 
pressure gradients, there is a minimum weight flow below which reverse 
flow exists. This is corrected by increasing weight flow. 

30) Not Used 

31) Not Used 

32) NORMAL COORDINATE OF LIFTING LINE IS NEGATIVE 
BLADE DATA DOES NOT CORRESPOND TO GEOMETRY OF DUCT 

This diagnostic occurs in Subroutine GE0MCL. Blade geometry was 
inputted incorrectly to program. This will produce a fatal error. 
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DETAILED PROGRAM DOCUMENTATION 


This section is intended to provide sufficient documentation to the user 
so that the internal operation of the program can be related to the analysis 
presented in the technical report (Reference 1). It is assumed that the user 
desiring this information has the required background in aerodynamics and 
computational fluid dynamics or access to the required technical support to 
understand the pertinent aspects of the program code as they relate to the 
theory. 

This section contains two major subsections, the propeller solution 
portion and the nacelle solution portion, respectively. These subsections 
contain, (1) an alphabetic list of the subroutines and external functions and 
a brief desctiption of each, (2) a more detailed description of each sub- 
routine in alphabetical order and, (3) a description of the label common 
blocks and variables used in alphabetical order. Flow charts and figures are 
provided in the subroutine descriptions where deemed necessary to understand 
the program structure and technical features. 

The subroutines and external functions are all described with the same 
format using the name of the subroutine with its argument list given as a 
title. A list of options and FORTRAN symbols used only in the named sub- 
routine are then given. Any special or additional theory used in the subrou- 
tine is presented but well known numerical methods are not described. 
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Propeller Program 


Within this section brief descriptions of the subroutines used in the 
lifting line portion of the analysis are presented followed by the labeled 
common blocks used in the analysis. The objective of a particular subroutine 
is noted, along with a list of symbols which are not in labeled common blocks 
and which are felt to be necessary for the understanding of the particular 
subroutines in question. These lists of symbols have been kept brief. A 
brief explanation of the theory is also included for selected subroutines. 
Generally, descriptions of the options which control the flow of a particular 
subroutine are also included. The subroutines are presented in alphabetical 
order. The labeled common blocks used in the analysis are also listed in 
alphabetical order with brief descriptions of each variable referenced. Brief 
flow diagrams of the major computational subroutines (GCWAKE, PR0P , S0LVEL and 
S0LVEN) are included in the description of each of these subroutines. 

List of Subroutines 


Name 

AFFIDC 

AF65A 

AIRFL 

AIRFLT 

AIRFMN 

AIR23 

AIR24 

ASSOC 

AVECTR 

BILINE 

BLDGE0 


Description 


Calculate blade activity factor and integrated 
design lift coefficient 

Calculate cascade airfoil data 

Control type of airfoil data to be used 

Airfoil package control routine 

Main airfoil package control routine for inter- 
facing interpolation, cascade data and isolated 
airfoil data 

Control lift and drag calculations for airfoil 
type number 23 

Control lift and drag calculations for airfoil 
type number 24 

Test input variable label 

Create a column vector from scalar input values 
Interpolation routine 
Calculate blade geometry 
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Name 

CALCGC 

CALWAK 

CASARF 

CASDAT 

CHKINP 

CL FACT 

C0MBWK 

CPITER 

CR0SSP 

CSCD1 

CTITER 

D0TP 

DRAG24 

DZR0AL 

ELIP2 

FINAIR 

FSQRT 

FVECTR 

GAUSS 

GCB0UN 


Description 

Calculate geometric influence coefficient 

Control flow of wake geometry calculations 

Control routine for analytical cascade correction 

Control selection of type of cascade data 

Check input parameters for obvious errors 

Calculate lift curve slope factor 

Calculate effective displacement of bound vortex 

Calculate blade angle for next power performance 
iteration 

Calculate cross product of two vectors 

Cascade airfoil data subroutine (reference ) 

Calculate blade angle for next thrust performance 
iteration 

Calculate dot product of two vectors 

Calculate drag coefficient for airfoil type 
number 24 

Calculate lift offset due to cascade influence 
Approximate Elliptic Integral of second kind 
Control final airfoil data calculation 
Calculate magnitude of three component vector 
Calculate blade forces 

Solve system of simultaneous linear equations 
(direct method) 

Calculate geometric influence coefficients for 
bound vortex 
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Name 


GCC0RE 

GCFILA 

GCWAKE 

G400LD 

INDVEL 

INTIAL 

IS0AFL 

IS0ARF 

LDDATA 

LIFT24 

LINEAR 
LINTER 
MV MULT 

MC0NE 

NSTAC0 

PAGE 

PCH0UT 

PERF0R 

PERI0D 


Description 

Calculate vortex core model 

Calculate geometric influence coefficient for 
trailing vortex filament 

Control basic flow of geometric influence 
coefficient calculations 

Calculate lift and drag using linearized airfoil 
data 

Calculate induced velocities 

Initialize data and print out selected quantities 
Control selection of isolated airfoil data type 
Control isolated airfoil data calculation 
Read in propeller data 

Calculate lift coefficient for airfoil type 
number 24 

Linear interpolation algorithm 

Control interpolation of data arrays 

Multiply single dimension vector with two 
dimensional vector 

Calculate Evvard Tip Relief Correction 

Calculate Mach cone intersection station index 

Advance output device to new page 

Punch spanwise distributions of aerodynamic 
quantities 

Calculate propeller performance 
Calculate propeller periodicity quantities 
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Name 


Description 


PERPRT 

PHICAL 

PLABEL 

PN 

PR0P 

PRDATA 

PRG400 

PRINTP 

PRTF15 

PRTF16 

PRTGCM 

PRTI15 

PRTI16 

PRTLF 

PRTLI 

PRTRZW 

PRWZW 

RDSCAL 

RDVECT 


Print spanwise distributions of aerodynamic and 
geometric quantities 

Calculate wake azimuth information 

Print a label 

Calculate a special function 
Main propeller subroutine 
Print label and vector 

Write output quantities for aeroelastic response 
analysis (reference 12) 

Print label and vector 

Print label and floating point variable vector, 
maximum of 15 

Print label and floating point variable vector, 
maximum of 16 

Print geometric influence coefficient matrix 

Print label and integer index for maximum of 15 
integers 

Print label and integer index for maximum of 16 
integers 

Print label field and single floating point 
variable 

Print single integer and label 
Print radial or axial wake coordinates 
Print wake geometry 
Controls read of scalar inputs 
Controls read of vector inputs 
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! 


Name 

READWR 

RELAXG 

REDMAT 

RWZWIN 

RWZWl 

RWZW7 

SBFUNC 

SETMAT 

SIEDEL 

S0LVEL 

S0LVEN 

S0LVIT 

SPLIN3 

STARC 

ST0RE 

SWPC0R 

THITER 

TITER 

TITLE 

UNBAR 


Description 

Reads a vector string and outputs it to printer 

Relaxation subroutine for circulation solution 

Read geometric influence coefficient in matrix 
form from disc 

Input wake geometry from cards 

Calculate classical or modified classical wake 

Calculate generalized wake 

Special function subroutine for cascade data, 
reference 11 

Set up geometric influence coefficient matrix 

Solve system of simultaneous equations (indirect 
method) 

Control linearized aerodynamic solution 
Control nonlinear aerodynamic solution 
Control solution procedure 
Interpolate with spline fit 

Convert design lift coefficient to equivalent 
camber angle 

Transfer data from one vector to another 
Calculate conical flow theory tip loss 
Control blade pitch iteration 

Extrapolate or interpolate on blade pitch angle 
versus Gj, or C^ 

Print title information 

Interpolation routine 
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Name 


Description 


UNINT 

VECT0R 

WECTR 

WAKM0D 

WRITGC 

ZER0AL 

ZER0GC 

Descript 

Object 
List of 


Theory 


Interpolat ion routine 

Compute velocity and velocity related quantities 
including induced velocities 

Compute velocity and velocity related quantities 
including induced velocities 

Modify wake geometry due to nacelle influence 
Write geometric influence coefficient to disc 
Calculate zero lift angle 

Set influence coefficient matrices to zero 
ion of the Subroutines Used in the Propeller Portion 


Subroutine AFFIDC (IT0T, DECL, B0D, SC0, AF, CLI, X, XB, R) 

To calculate the blade activity factor and integrated design 
lift coefficient. 

Symbols 

AF Blade activity factor 

CLI Integrated design lift coefficient 

The activity factor is defined as 


“• i'm a f - ft ) 
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where C/D is the blade chord to propeller diameter ratio and x 
is the spanwise location along the blades. 

The integrated design lift coefficient is 


CLI 


4( 1-sco) 


/ 

SCO 



where SC0 is the root cutout and Cj^ is the section design 
lift coefficient . 


Subroutine AF65A (Argument List) 

Object Compute airfoil lift and drag from cascade correlations. 

List of Symbols 

Arugment List 


AMACH 

M 1 

9 

Upstream Mach number 

(INPUT) 

ALP 

a 

y 

Angle of Attack 

(INPUT) 

TM 

t 

y 

Maximum Airfoil Thickness 

(INPUT) 

THETA 

0 

y 

Pitch angle 

(INPUT) 

CB 

^B 

9 

Design Lift Coefficient 

(INPUT) 

S0LD 

g/c 

y 

Cascade Solidity 

(INPUT) 

CL (1) 

C L 

9 

Lift Coefficient 

(OUTPUT) 

CD (1) 

°D 

9 

Drag Coefficient 

(OUTPUT) 

iscade 

Correlation 




ALPS 

a s 

9 

Stagger Angle (degrees) 


ALP1 

“1 

9 

Inlet Air Angle (degrees) 


PHIC 

♦c 

9 

Camber Angle 


AKDELS 

k 6s 

9 

Shape Parameters 


AMSIG 

M o 

9 

Camber Parameter 


B 

b 

9 

Exponent 


DEL 

6 oo 

9 

Deviation Angle, <f» c = 0 


AKDELT 

K 6t 

9 

Thickness Parameter 


DELO 

«o 

9 

Deviation Angle (degrees) 


AIOO 

*oo 

9 

Incidence Angle, $ = 0 (degeees) 

AN 

n 

9 

Power 


AKIT 

K it 

9 

Thickness Parameter 
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AIMO 

1 mo 

f 

Minimum Loss Incidence Angle (degrees) 

D 

D 

9 

Diffusion Parameter 

ZL0SM 

Z sm 

9 

Minimum Loss Coefficient 

A1NC0 

i 

9 

Incidence Angle (degrees) 

ZL0SS 

z s 

9 

Loss Coefficient 

DALST 


9 

Stall Angle Correction 

ALPH2 

a 2 

9 

Exit Air Angle 

RH0CX 

<pU s ) 2 /(P U s>i 

9 

Mass Flow Ratio 


Additional 

Symbols 



(a 2> Z s ) * 

(C L , Cjj) 



T1 

T 0l/ T l 

9 

Upstream Total Static Temperatures 

AMACH2 

M 2 

9 

Downstream Mach Number 

T2 

T 02^ T 2 

9 

Downstream Total Static Temperature 

P02P01 

P 02/ P 01 

9 

Total Pressure Ratio 

P2P1 

P 2 / P l 

9 

Static Pressure Ratio 

FS 

F s 

9 

Streamwise Force Coefficient 

FP 

F p 

9 

Tangential Force Coefficient 

WS 

W s 

9 

Streamwise Induced Velocity 

WP 

w 

p 

9 

Tangential Induced Velocity 

ALIND 

a i 

9 

Induced Flow Angle 

ANG 

a 

9 

Angle of Attack 


Theory See section of reference 1 entitled: "Cascade Airfoil Data". 


Subroutine AIRFL 

Object Controls which type of airfoil data will be used, either 

isolated data or isolated data corrected for cascade effects. 

IDL = 0 print title 

IDL * 0 obtain airfoil data 

ICASDE = 0 obtain isolated airfoil data 

1CASDE = 1 obtain isolated airfoil data corrected for cascade 

effects 
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Subroutine AIRFLT ( IQ , IFQ , IDQ , ICASDQ , ALPHQ , THETAQ , 
TAUBQ , ZMQ ,DECLQ ,H0BQ , CL3Q , CDQ) 


Object 


Options 


Object 


Options 


Argument 

IC 

IFL 

I 

NSTAT 

RADCAS 

RDTRAN 

XMTIP 

RSC 

SMACH 


Control combinations of airfoil data characteristics to be 
obtained . 

IDQ = 0 print title 
IDQ = 1 obtain and 
IDQ = 2 obtain C L only 

Subroutine AIRFMN ( IC, IFL, I , NSTAT, RADCAS, RDTRAN, 

XMTIP , RSC , SMACH ,THK ,THETAG ,DESCL ,SIGMAX ,TAU , 

ALP , CL , CD , FTRAN 1 , FTRAN2 , TAUEXP ) 

Control calculation of airfoil and cascade data and the inter- 
polation between isolated and cascade data when requested. 

NG400 * 0 use linearized airfoil data from aeroelastic response 
analysis (reference 12) 

ICAS = 0 use isolated airfoil data 

ICAS * 0 use cascade data 

If cascade data is used, interpolation between cascade data and 
isolated data is controlled by value of RDTRAN. 

List 


controls lookup of C^ alone or C^ and C^, or output of 
title information 
airfoil type index 
blade station index 

blade station limit for tip loss model correction 
outermost radial location for direct cascade data 
application 

outermost radial location for direct cascade/isolated airfoil 

interpolation procedure application 

tip Mach number 

radial station 

section Mach number 
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THK section thickness ratio 

THETAG section geometric pitch angle with the plane of rotation 

DESCL section design lift coefficient 

SIGMAX section solidity 

TAU section gap to chord ratio 

ALP section angle of attack 

CL section lift coefficient 

CD section drag coefficient 

FTRAN1 section cascade/airfoil interpolation scaling function value 

for c t 

RTRAN1 section cascade/airfoil interpolation scaling function value 

f ° r C d 

TAUEXP exponent used in interpolation function 


Subroutine AIR23 


Object 

Calculate 

tabulated 

Manoni airfoil characteristics from internally 
data bank using transonic simularity rules. 

Options 

IDL = 1 

obtain 


IDL = 2 

obtain 


IDL = 3 

dummy feature 


Theory 

Using transonic simularity rules, empirical data has been 
reduced to a set of tabulated coefficients which can be used to 
reconstruct the airfoil characteristics (C^ and C^) for a 
wide range of parameters. 

Subroutine AIR24 


Object 

Control 
for the 

selection of calculation of lift or drag coefficients 
published NACA data. 

Options 

IDL = 1 

obtain lift 


IDL = 2 

obtain drag 


IDL = 3 

dummy feature 
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Subroutine ASS0C (*,S ,FDUM,S1 ,S2 ,X) 

Object To transfer input value of a dummy parameter to its correct 

allocation if the input label matches one in the argument list. 
If transfer is made, returns to labeled statement in calling 
routine . 

Argument List 

S input label as read 

FDUM input dummy parameter 

SI ,S2 input label list 

X allocation for transfer of dummy parameter 


Subroutine AVECTR (Al,A2,A3,V) 

Object load three scalars into a vector of length three. 

Argument List 

A1,A2,A3 input scalars 
V output vector 


Subroutine BILINE (T ,1 ,XI ,YI ,Z,K) 

Object Bivariant or univariant interpolation on input vectors using 

various interpolation options. 

Options T(l + 1) ■ 0 use first table value 

T(l + 1) = 1 use linear interpolation 

T(l + 1) ■ 2 use third order interpolation 

T(l + 3) = non zero, requests bivariant interpolation 


Argument List 

T = vector with interpolation data 
I * starting location for data table 
XI * input x for interpolation 
YI * input y for interpolation 
Z * output value 
K * error code 
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Theory 

Using either bivariant or univariant data, this subroutine will 
interpolate on the data using standard interpolation schemes as 
noted above - see listing for more detail. 

Subroutine BLDGE0 (IC) 

Object Rotate input lifting line segment geometry to the requested 

blade angle about the pitch axis and compute the inflow station 
geometry for this blade angle. Print the table of the lifting 
line segment and center coordinates. 

List of Symbols 

IC = print control flag 

DR00P ® angle between the coordinate origin (hub center) and 
the axial displacement for the coordinate point in 
question for blade element center. 

DR00PB = angle between the coordinate origin (hub boundary) and 
the axial displacement for the coordinate point in 
question for the blade element boundary. 

Theory Standard geometric operations applied to the input geometry to 

rotate the coordinates to a different blade angle position. 

Subroutine CALCGC (Argument List) 


Object To calculate the geometric influence coefficients for a selected 

field point in the cylindrical coordinate system. 


Options 

NCBWAK 

= 0 

no compressibility effects on the bound vortex 
induced velocity calculation 


NCBWAK 

- 1 

compressibility effects included on the bound 
vortex induced velocity calculations 


NCBWAK 

= 2 

compressibility effects included in the bound 
vortex induced velocity calculation, except for the 
calculation for the particular blade bound vortex 
system on itself 
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IJUNK * 0 
I JUNK = 1 
IDEBUG = 0 
IDEBUG = 1 
NC0MPRS = 0 
NC0MPRS = 1 


no vortex core model 
vortex core model used 
no intermediate printout 
intermediate printout 
no wake compressibility model 
wake compressibility model used 


NCFL0W * 0 wake compressibility applied only on vortex seg- 
ments from inflow sections with Mach numbers 
greater than 1.0 

NCFL0W - 1 relaxes the above restriction 
Argument List 

NC0MP * compressible wake option switch 

IBIP = blade index 

LLINK « blade position index 

NBX = number of blades 

IT0T * number of blade element segments 

LT0T = number of blade positions 

KT0T = number of blade element boundaries 

MT0T * number of filament segments 

KTRUCT 35 number of filaments for tip vortex rollup model 
JTRUCT * wake rollup truncation angle 
JTRUCI = inboard wake truncation angle 
DPSIBR * azimuth interval in radians 

FS * sign (-1.0) for axial induced velocity calculation 

RSCI 33 radial location of field point 

PHICI * lag angle of field point 

ZSCI * axial location of field point 

CPD = cosine of blade element angle 

SPD = sine of blade element angle 

RSBB = blade element boundary radial position 

ZSBB * blade element boundary axial position 

PHIBB * blade element boundary azimuth position 

C0SLB = blade element boundary cosine of azimuth position 

SINLB * blade element boundary sine of azimuth position 

CMACH e local blade element Mach number 

MU = local blade element advance ratio 

DTIPM * blade tip Mach number 

V0RC0R * vortex core radius 
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Theory 


Object 

Options 


Object 


Options 


Using the Biot-Savart relationship for the geometric influence 
coefficients for straight line vortex segments (see Appendix A 
of reference 1), these coefficients are calculated, summed and 
stored in cylindrical system form. See the technical approach 
for the propeller analysis in reference 1. 


Subroutine CALWAK (IWK) 


Control selection of wake models 


IWAK0P = 0 
IWAK0P = 1 
1WAK0P = 2 
IWAK0P * 3 
NACWAK * 0 
1PR0PT * 0 


classical wake geometry 

classical wake geometry 

input geometry 

generalized wake geometry 

modify wake geometry by nacelle influence 

print wake geometry 


Subroutine CASARF 

Controls the calculation procedure for the computation of the 
isolated airfoil data corrected for cascade effects. 

IDL ■ 1 calculate C^ and correct for cascade influence 
IDL = 2 calculate Cp 
IDL = 3 dummy feature 


List of Symbols 

CLKFAC “ cascade correction scaling factor 

THETAZ = geometric blade angle corrected for camber and angle of 
zero lift 


Theory 

The lift coefficient is corrected for cascade effects by 
applying an analytical correction for the cascade influence on 
flat plates. The details of this correction are presented in 
reference 1 . 
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Subroutine CASDAT (Argument List) 

Object Control selection of cascade data source. 

Options ICAS * 0 terminate execution of code 

ICAS = 1 cascade data using correlation from reference 13 
ICAS ® 2 cascade data using model of reference 11 

Argument List 

ICAS = option switch for type of cascade data 
MACH = local section Mach number 

AL * local section angle of attack 

THK = local section thickness ratio 

THET * local section blade angle 

DESCLP = local section design lift coefficient 
SIG = local section solidity ratio 

CL * local section lift coefficient 

CD - local section drag coefficient 

Subroutine CHRINP 

To check input data for correct input values on selected items 
and make sure items that are necessary for successful execution 
are input . 

Subroutine CLFACT ( THETA, TAUB,CLKF AC) 

Obtain tabled value of lift curve slope scaling factor for input 
parameters . 

List of Symbols 

THETA * geometric blade angle corrected for camber and angle of 
zero lift 

TAUB ■ gap-to-chord ratio 

CLVFAC * lift curve slope scaling factor 

Theory Table of data for the lift curve slope scaling factor was 

derived analytically for flat plates by Weinig. The table of 
data was obtained from reference 4. 


Object 


Object 
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Subroutine C0MBWK (Argument List) 

Object To compute an effective axial displacement correction on the 

bound vortex location for incorporating the phase shift of the 
induced influence of the bound wake on an inflow station. 

Argument List 

RA = radial location of the midpoint of a bound vortex 
segment 

RB = radial location of the inflow station 
DZ = axial location of a bound vortex segment 
MU = flight speed divided by tip speed 
DTIPM = tip Mach number 

DPSI = azimuthal position of bound vortex segment 
Z = effective axial displacement 

PHI = effective phase angle 

Theory The first real positive root of a transcendental equation is 

solved by a simple root searching algorithm and a Newton-Raphson 
iteration procedure. This root represents an effective phase 
angle associated with the finite delay time for a signal to 
reach an inflow station point if it originally emanated from a 
bound vortex source inside a zone of silence of a Mach cone. 

Subroutine CPITER (Argument List) 

Object To calculate either a first guess on blade angle or subsequent 

iterations values for the blade angle when the power performance 
iteration feature is requested. 

Options NCP * 1 obtain first guess from tabled values 

NCP = 2 obtain second value from one of two methods 

NCP * 3,..., 10 obtain all subsequent values by linear inter- 
polation or extrapolation from previous iteration 
information 

DPDT =0.0 for NCP = 2 use a new blade angle of ±1.5 degrees 
from the first iteration value 

DPDT * 0.0 for NCP = 2 use DPDT as the linear slope to define 
the next blade angle 
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Argument List 


IT0T = number of blade elements 

BL = number of blades 

ZJI * advance ratio 

RSC = inflow station segment radial centers 

RSB = inflow station boundary radial locations 

DESCL = design lift coefficient 

B0D = chord over diameter ratio 

CPWANT = requested power coefficient 

DPDT = input linear slope of the power coefficient versus 
blade angle relationship 
RAD = blade radius 

Theory Using standard interpolation and extrapolation techniques, the 

blade angle for each iteration of the power performance itera- 
tion is determined. For the first iteration tabled values are 
used if requested and for the second iteration one of two 
methods can be used to determine the new blade angle. All 
subsequent iteration values are obtained using linear extrapola- 
tion or interpolation based on previous iteration information. 

Subroutine CR0SSP (V ,A,R1 ,R2 ,R3) 

Object Calculate cross product of two vectors and place results in 

three scalars. 
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IND = controls selection of airfoil type 

CB = section camber 

TH = section thickness ratio 

S0L = section solidity 

ST = section pitch angle 

RN = Reynolds number (fixed at 500000) 

AM - section angle of attack 

CL = section lift coefficient 

CD = section drag coefficient 


Subroutine CTITER (CTWANT ,DCTDT) 

Object To calculate either the first or all subsequent blade angles 

when the thrust iteration has been requested. 

Options NCP = 1 use a fixed value of 60.0 degeees for the first guess 

NCP = 2 use one of the two methods to determine the second 
value 

NCP = 3,.. ,10 use linear interpolation or extrapolation based 
on the previous iterations to determine the next value 

DCTDT =0.0 for NCP = 2 use a new blade angle of ±1.5 degrees 
from the first iteration value 

DCTDT * 0.0 for NCP = 2 use the value of DCTDT as the linear 
slope to define the next blade angle 

Using standard interpolation and extrapolation techniques, the 
blade angle for each iteration of the thrust iteration cycle is 
determined. For the first iteration, a fixed value is used if 
requested and for the second iteration one of two methods can be 
used to determine the new blade angle. All subsequent itera- 
tions values are obtained using linear interpolation or extrap- 
olation based on previous iteration information. 

Subroutine D0TP (V,A) 

Object Calculate dot product of two vectors, each of length three. 
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Subroutine DRAG24 


Object Calculate drag coefficient from the NACA airfoil data. 

Theory Using linear interpolation techniques applied to a set of 

tabulated airfoil drag coefficients which are functions of Mach 
number, angle of attack, thickness to chord ratio and design 
lift coefficient, the drag coefficient is found for a specified 
combination of the above parameters (reference 1). 

Subroutine DZR0AL (THSTAR,TAUB,THETAT,DELA0L> 

Object Determine increment in angle of zero lift for cascade correction 

procedure. 


Argument List 


THSTAR = effective blade camber angle 

TAUB = gap-to-chord ratio 

THETAT * effective blade angle 

DELA0L = increment in angle of zero lift 

Theory Using trivariate interpolation techniques, the increment in the 

angle of zero lift is determined from a set of tabulated data 
for double circular arc aifoils (reference 5). 

Subroutine ELIP2 (X) 

Object Calculate approximation for elliptic integral of the second 

kind . 


Argument List 


X * argument of approximation function 


Object 

Options 


Subroutine FINAIR 


Control final calculation procedure for airfoil characteristics. 


IDEBUG * 0 
IDEBUG > 0 

ICAS * 0 
ICAS * 0 


no printout 

printout of final values for the airfoil 
characteristics is requested 
no cascade airfoil data is used 

cascade airfoil data is used out to a requested 
radial location 
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Function FSQRT (X,Y,Z) 


Object Calculate magnitude of vector of components X, Y, and Z 

Theory FSQRT = (X 2 + Y 2 + Z 2 ) 1/2 

Argument List 

X,Y,Z » components of vector 
FSQRT * resultant 

Subroutine FVECTR 

Object Calculate components of blade forces at each blade inflow 

station. 

Options IDEBUG * 0 request printout of component forces 

List of Symbols 

ALFLC * lift force in the chordwise direction 

ALFLN 35 lift force in the normalwise direction 

ALFLS * lift force in the spanwise direction 

ALFDC * drag force in the chordwise direction 

ALFDN = drag force in the normalwise direction 

ALFDS = drag force in the spanwise direction 

Theory Using a calculated blade section lift and drag coefficients and 

the appropriate aerodynamic quantities, the forces in the blade 
element coordinate system are calculated and these forces are 
transformed to the cylindrical coordinate system. 

Subroutine GAUSS (Augument List) 

Object Solve a system of simultaneous linear equations in matrix form 

using a direct solution technique. 

Argument List 

NR0WM = row dimension of the coefficient matrix 
N * number of rows and columns used in the matrix solution 
A = coefficient matrix 

B * constant vector and on output contains the solution 
vector 
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Theory 


Object 


Argument List 


Theory 


DET * mantissa of the value of the determinant in base ten 
IDET = integer power to the base ten of the determinant 
LSING = singularity flag 

Using a standard Gauss-Jordan reduction method a system of 
simultaneous linear equations is solved and the determinant of 
the matrix is calculated. 


Subroutine GCB0UN (RSC,CP ,SP ,ZSC,RSB,ZSB,PHIB, 
C0SLB,SINLB,CMACH,MU,TIPM,VC0R) 

Calculate geometric influence coefficients at a specified load 
point for the bound vortex segments which represent the 
propeller blade lifting line. 


RSC * radial position of load point at which induced influence 
is to be calculated 

CP * cosine of azimuthal position of load point 
SP * sine of azimuthal position of load point 
ZSC * axial position of load point 

RSB * radial location of bound vortex segment endpoints 
PHIB = azimuthal position of bound vortex segment endpoints 
C0SLB ■ cosine of azimuthal position of bound vortex segment 
endpoints 

SINLB * sine of azimuthal position of bound vortex segment 
endpoints 

CMACH ® mach number of bound vortex segment endpoints 
MU * advance ratio of propeller disk 
TIPM ■ tip mach number of propeller disk 
VC0R * vortex core radius 

Uses the potential flow solution for the induced influence of a 
finite length straight vortex filament, formulated in a 
cylindrical coordinate system. See Appendix B, Volume 1. 
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Object 


Argument List 


Theory 


Subroutine GCC0RE ( IB , IFILA, IDEBUG , IFLAG , RARB , 

ZAZB , RASQ ,RBSQ , CP , DSCA, DSCA, VC0R) 

Tabulate induced influence of a vortex segment if the load point 
is within the core radius. 


IB 

3S 

type of vortex segment indicator (bound 

or 

trailing) 

IFILA 

= 

vortex segment number 



IDEBUG 

ss 

debug output 

trigger 



IFLAG 

= 

flag 




RARB 

as 

intermediate 

geometric quantity needed 

for 

the 



calculation, 

see Appendix B, Volume I 



ZAZB 

= 

intermediate 

geometric quantity needed 

for 

the 



calculation, 

see Appendix B, Volume I 



RASQ 

= 

intermediate 

geometric quantity needed 

for 

the 



calculation, 

see Appendix B, Volume I 



RBSQ 

= 

intermediate 

geometric quantity needed 

for 

the 



calculation, 

see Appendix B, Volume I 



CP 

= 

intermediate 

geometric quantity needed 

for 

the 



calculation, 

see Appendix B, Volume I 



DSCA 

= 

intermediate 

geometric quantity needed 

for 

the 



calculation, 

see Appendix B, Volume I 



DSCA 

s 

intermediate 

geometric quantity needed 

for 

the 



calculation, 

see Appendix B, Volume I 



VC0R 

= 

vortex core 

radius 




If the load point at which the induced influence due to a vortex 
segment is within a specified core radius, the induced influence 
is modeled by a solid body rotation model. This removes the 
singular behavior due to a purely potential flow vortex. 
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Subroutine GCFILA (K ,NC0MPT , IBIWK ,LBT0T , LLIWK , 

RSCI , CP , SP , ZSCI , KTRUCT , JTRUCT , JTRUCI ,MT0T ,KT0T , 
LT0T,C0SLB,SINLB,DTLPM,V0RC0R) 

Object Calculate geometric influence coefficients at a specified load 

point for the trailing vortex segments of a specified filament 
which represent the wake geometry. 


Options 


NC0MPT 

= 0 

no compressible wake 

NC0MPT 

* 0 

use compressible wake 

KTRUCT 

* 0 

wake rollup model used beyond this filament index 



value 

V0RC0R 

* 0 

vortex core option requested 


Argument List 

K = vortex filament index 

NC0MPT - option control for compressible wake model 

IBIWK = blade index for wake 

LBT0T = blade and rotor position index 

LLIWK = rotor position index of wake 

RSCI * radial position of load point 

CP = cosine of azimuthal position of load point 

SP = sine of azimuthal position of load point 

ZSCI = axial coordinate of load point 

KTRUCT = filament index defining boundary for tip vortex rollup 
model 

JTRUCT = filament azimuth position index for wake truncation 
associated with tip rollup model 

JTRUCI = filament azimuth position index for wake truncation for 
inboard wake 

MT0T * number of segment endpoints 

KT0T * number of filaments 

LT0T * number of rotor positions 

C0SLB * cosine of filament azimuth position at blade 
SINLB * sine of filament azimuth position at blade 
DTIPM ■ reciprocal of blade tip Mach number 
V0RC0R - vortex core radius 

Theory Uses the potential flow solution for finite length straight 

vortex segment, formulated in a cylindrical coordinate system. 
See Appendix B, Volume I. 
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Subroutine GCWAKE 


Object Controls flow of the geometric influence coefficient calcula- 

tions. A flow diagram is presented in figure 12. 

Options IWAK0P 3 0 standard wake model (modified classical wake) 

IWAK0P 3 1 classical wake model 
IWAK0P 3 2 input wake geometry 
IWAK0P 3 3 generalized wake model 

NACWAK * 0 obtain nacelle corrections for wake geometry 
IPR0PT * 0 print wake geometry 

Subroutine G400LD (I , ALPHA, CL , CD) 

Object Calculate lift and drag airfoil characteristics obtained from an 

external source. 

Argument List 

I 3 blade segment index 

ALPHA 3 section angle of attack 

CL 3 section lift coefficient 

CD 3 section drag coefficient 

Theory The lift and drag characteristics for each blade element segment 

are obtained from an external source by curve fitting a 
quadratic polyraonial about the angle of attack at each station. 
The lift and drag are reconstructed in this subroutine for use 
in the solution procedure. 

Subroutine INDVEL 

Object Compute induced velocities. 

Option IDEBUG * 0 print the calculated induced velocities 

Theory Reading the stored geometric influence coefficients in a defined 

order from a disc, the induced velocities are calculated by 
multiplying and summing over the appropriate indices of the 
matrix quantities . 
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Subrout ine INTIAL 


Object To calculate selected data and printout initial input data. 

Subroutine IS0AFL 

Object Control selection of type of isolated airfoil data tables to use 

(Manoni or NACA) . 

Options IDL = 1 lift only 

IDL = 2 drag 
IDL = 3 dummy feature 
IFL * 23 use Manoni data 
IFL = 24 use NACA data 
IFL < 23 or > 24 use Manoni data 

Subroutine IS0ARF 

Object Control the flow of isolated airfoil data module. 

Options IDL * 0 return to calling routine 

IDL * 1 compute lift 
IDL = 2 compute drag 

Subroutine LDDATA 

Object Read in required initial propeller input data and calculate 

selected data from input quantities. 

Subroutine LIFT24 

Object Calculate lift coefficient from tabulated airfoil data tables 

containing the NACA data. 

Theory Using linear interpolation techniques this subroutine computes 

the lift coefficient from a table of NACA airfoil data which is 
a function of Mach number , angle of attack, design lift 
coefficient and thickness to chord ratio (reference 1). 
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Subroutine LINEAR (Argument List) 

Object Interpolate and extrapolate linearly on an input set of data. 

Argument List 

NW - print output file number 

N = number of data points in the interpolation vectors 

XIN * independent data vector 

YIN = dependent data vector 

X0UT = requested interpolation point 

Y0UT * interpolated value 

L0FF - interpolation flag 

Theory Simple linear interpolation algorithm; if off scale on the low 

or high end, the flag is set (1 or 2 respectively) and the 
boundary slope is used to extrapolate to the requested inter- 
polation point. 

Subroutine LINTER (Argument List) 

Object Control the interpolation of selected data arrays from input 

interpolation tables . 

Options NERR s 1 use lower boundary value 

NERR = 2 use upper boundary value 

Argument List 

N - number of data points in the interpolation vector 

K = number of requested interpolation points 

XIN * independent interpolation vector 

YIN ■ dependent interpolation vector 

X = vector of requested interpolation points 

Y 55 vector of interpolated values 
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Subroutine MC0NE (Augument List) 


Object Calculate Eward Tip Relief Correction. 

Options NEVARD =0 no correction 

NEVARD ■ 1 used tabled values 
NEVARD = 2 use equation 

Argument List 

NEVARD = option control for method of calculation 

IP = propeller index 

L = propeller position index 

RSC = radial location of inflow station 

C = chord 

R " propeller radius 

IT0T = number of inflow stations 

CNSECT = fraction of chord 

XMTIP = tip Mach number 

NSTAT = station index for boundary for Mach cone intersection 
XKC0NE = tip relief correction vector 
RT1P = tip radius 

Theory See section of reference 1 entitled: "Eward Tip Relief for 

Propellers". 

Subroutine MV MULT (V ,M,Rl ,R2 ,R3) 

Object Multiply a three by three matrix by a vector of length three. 

Store the result in three separate scalars. 


Argument List 


V “ vector 

M = matrix 

R1,R2,R3 ■ resultant scalars 
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Subroutine NSTAC0 (Argument List) 


Object 
Options 
Argument List 

Theory 


Object 

Argument List 
Object 

Argument List 


Calculate Mach Cone intersection station index, 
IDEBUG > 0 requests printout of selected data 


NSTAT = station index where intersection occurs 
XMTIP = tip Mach number 
IP = propeller index 
L = propeller position index 

Using relative geometry the angle XNETA is computed for each 
station and compared with the Mach cone angle, BETA, until the 
station where the intersection of the Mach cone with the 
specified fraction of the blade chord is determined. 


Subroutine PAGE (NWRITE) 
Print new page. 


NWRITE = print output file number 


Subroutine PCH0UT (NUNIT) 

Output spanwise distributions of aerodynamic and geometric 
quantities to specified output device number. 


NUNIT = output file number 
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Subroutine PERF0R 


Object Compute and print out the propeller performance parameters. 

Theory Uses standard integration techniques to obtain the integrated 

thrust and power from the blades force components in .the 
cylindrical coordinate systems. 


Thrust 


/ 


Tip 
Root 


F z cdr 


Torque 


-/ 


Tip 
Root 


F^rcdr 


where 

c = chord 

F z = axial force per unit area 
F^ = tangential force per unit area 
r = local blade radius 


Object 

Theory 


Subroutine PERI0D 


Calculate propeller disc periodicity and related quantities. 

For single propeller disc configurations, the geometric 
relationship between the wake and blades is fixed. However, for 
a coaxial propeller, the geometric relationships are periodic 
with half-blade spacing if the number of blades and rotational 
speeds are equal. For unequal blades and/or unequal rotational 
speed, the relationship defining the periodicity (t) of the wake 
and blade geometry is 


t 


2ir 


b max^ ®1 


+ fy) 


where and fl 2 are t ^ ie rotational speeds of the two propellers, 
and b max is the maximum number of blades of the two propellers. 
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Subroutine PERPRT (NW,L,IP) 


Object Print spanwise distributions of aerodynamic and geometric 

quantities to a specified output unit. 


Argument List 


NW = output unit number 
L = propeller position index 
IP = propeller index 


Subroutine PHICAL (MT0T,DPSl) 

Object Compute wake azimuth positions and trigonometric relationships 

for all propeller azimuth positions. 


Argument List 


MT0T « number of wake filament segments 
DPSI = azimuthal increment 

Theory Using the requested blade azimuth increment, the wake azimuth 

position and the sine and cosine functions for such are computed 
and stored. The fact that the sine and cosine functions are 
periodic is made use of to reduce the actual number of sine and 
cosine values which are stored. 

Subroutine PLABEL (NWRITE ,NSLB ,NSLF ,NSP ,NL , LABEL) 

Object Print a label field. 

Argument List 


NWRITE = 
NSLB 
NSLF 
NSP 

NL 

LABEL = 


output device number 

number of lines to skip before label is output 
number of lines to skip after label is output 
number of units of 10 spaces to skip before label is 
output 

length of label in increments of 6 
label vector 
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Function PN (N,R,S) 


Object 

Options 

! 

Object 


Options 

! 


Object 


Compute special functions of R and S. 

N * 0 PN * 1 

N = 1 PN = 1 / (R-S) 

N > 2 PN = R N “ 2 


Subroutine PR0P 

Control sequencing of propeller lifting line solution procedure. 
A flow diagram is presented in figure 13. 

CPI * 0 requests power iteration (overrides thrust iteration) 
CT1 * 0 requests thrust iteration 


Subroutine PRDATA (T,IUNIT,N,X) 

Output label vector and floating point vector to specified 
output unit number using 7A6/8E10 .4/8E10 .4 format. 


Argument List 


T * label vector 

1UN1T “ outut unit number 

N = length of floating point vector 

X * floating point vector 


i 

I 
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Subroutine PRG400 


Object Output required velocity quantities to specified unit number in 

a format compatible with an aeroelastic response analysis 
(reference 22 :) . 

tZ 

Subroutine PRINTP ( NWRITE , IFC0DE , N , FDATA ,L ABEL ) 

Object Output label vector and floating point vector with specified 

format of form 7A6,10F9.X. 

Argument List 

NWRITE = output unit number 

IFC0DE * format code index 

N = length of floating point vector 

FDATA = floating point vector 

LABEL = label vector 


Subroutine PRTF15 (NWRITE , IFC0DE ,N, FDATA, LABEL) 

Object Output label vector and floating point vector of specified 

length to a specified unit number using a format of form 
3A6.15F8.X. 

Argument List 


NWRITE “ output unit number 

IFC0DE ■ format code index 

N “ length of floating point vector 

FDATA ■ floating point vector 

LABEL = label vector 


Subroutine PRTF16 ( NWRITE, IFC0DE ,N, FDATA, LABEL) 

Object Output label vector and floating point vector of specified 

length to a specified unit number using a format of form 
3A6.16F7.X. 

Argument List 


NWRITE ” output unit number 

IFC0DE “ format code index 

N = length of floating point vector 

FDATA = floating point vector 

LABEL * label vector 
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I 

I 

I 

I 

j Object 


Subroutine PRTGCM (NWRITE, I0P,NR0W,NC0L,GCMAT) 

Output two-dimensional influence coefficient matrix to specified 
unit number, in one of two formats. 


Argument List 


NWRITE ■ output unit number 

I0P = format option index 

NR0W * row dimension of matrix 

NC0L s column dimension of matrix 

GCMAT * matrix 

Options I0P ** 0 use format 15F8.5 

I0P - 1 use format 15F8.3 


Subroutine PRTI15 (NWRITE, N, LABEL) 

Object Output to specified unit number a label and integer index with 

format of form 2A6,I7,14I8. 


Argument List 

NWRITE = output unit number 
N * length of integer index 

LABEL - label vector 

Subroutine PRTI16 (NWRITE ,N, LABEL) 

Object Output to specified unit number a label vector and integer index 

with format of form 3A6,16,1517. 

Argument List 

NWRITE *= output unit number 
N * length of integer index 

LABEL * label vector 
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Subroutine PRTLF (NWRITE, F, LABEL) 


Object Output label vector and single floating point scalar with format 

of form 5A6.F15.5. 


Argument List 


NWRITE m output unit number 
F * floating point scalar 

LABEL = label vector 


Subroutine PRTLI (NWRITE ,N, LABEL) 

Object Output label vector and single integer value with format of form 

5A6.I10. 


Argument List 


NWRITE * output unit number 
N * integer value 

LABEL * label vector 


Subroutine PRTRZW (NWRITE ,RSB, PHI .RZW.MT0T ,KT0T, IRZ) 

Object Output floating point vectors to specified unit number with 

integer indices with specified format. 

Options IRZ “ 0 use format of form 16F7.3 

IRZ = 1 use format of form 16F7.2 

Argument List 

NWRITE = output unit number 

RSB ■ floating point vector 

PHI = floating point vector 

RZW “ floating point vector 

MT0T * length of column vector 

KT0T “ length of row vector 

IRZ = option for format type 
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I 


Subroutine PRWZW (RSBB) 

Object Control print of wake coordinates. 

Argument List 


RSBB * blade element segment boundary radius 


Subroutine RDSCAL (IVPRNT) 

Object Read scalar input parameters necessary to run program. 

Options NG400 * 0 read data from external unit source for necessary 

scalar input for coupling with aeroelastic response 
analysis (reference 12) 


Argument List 


IVPRNT = print option for vector input listing 


Object 


Options 


Subroutine RDVECT (IVPRNT) 

Read vector inputs necessary to run. 

IVPRNT * 0 no listing of inputs as read in 

NG400 t 0 read vector inputs from external source for coupling 
with aeroelastic response analysis (reference \Z ) 


Argument List 

IVPRNT = option flag to terminate listing of input vector as 
read in 
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Subroutine READWR (*,NREAD,NWRITE,IT,S,ST,X) 
Object Read input vectors if label fields match. 


Argument List 


* 

= 

return 

i 


NREAD 

= 

input 

unit 

number 

NWRITE 

= 

output 

unit 

number 

IT 

= 

length 

of vector 

S 

= 

input 

label 

field 

ST 

* 

test label 

field 

X 

= 

input 

vector 


Subroutine RELAXG ( ITER , IC0V , RELAXF , GAMMA) 

Object Relaxation of circulation solution and convergence flag set in 

this subroutine. 

Argument List 


ITER = iteration index 

IC0V = convergence flag 

RELAXF = relaxation factor 

GAMMA = temporary single dimension variable for solution 
storage 

List of Symbols 

CIRC = current circulation stored in this matrix 

SAVCIR = previous iteration circulation stored in this matrix 

Subroutine REDMAT (Argument List) 

Object Read all geometric influence coefficients from disk and create 

the geometric influence coefficient matrix. 


Argument List 


NM = dimension for the maximum number of rows in the matrix 

MSIZE * number of rows used in the matrix 

GCN * geometric influence coefficient matrix 
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Subroutine RWZWIN (IWK) 


Object Input wake geometry coordinates. 

List of Symbols 

PHI = wake azimuth position 
RW * radial coordinates of wake 
ZW - axial coordinates of wake 
IWK = wake index number 

Subroutine RWZWl (IWK) 

Object Compute wake coordinates for classical or modified classical 

wake model. 

Options IWAK0P = 0 use prescribed inflow distribution and the momentum 

induced velocity to define the wake geometry 
(modified classical wake) 

IWAK0P *= 1 use the freestream velocity and the momentum induced 
velocity to define the wake geometry (classical 
wake) 

List of Symbols 

PHI ** wake azimuth position 

RW “ radial coordinates of wake 

ZW * axial coordinates of wake 
IWK * propeller wake index number 
VWAKE = propeller wake transport velocity 

Theory Classical wake. 

VWAKE * -VKTAS*1 .688+VIM0M 


Modified classical wake 


VWAKE « VZER0B(l)+VIM0M 


Subroutine RWZW7 (IWK) 

Object Compute generalized wake geometry. 

Options I0PT = 0 linear fit for wake geometry between tip filament 

and non-rolled up filaments outboard of the outer 
sheet filament 

I0PT * 1 parabolic fit for filaments as noted above 

List of Symbols 

PHI * wake azimuth position 
RW * radial coordinates of wake 
ZW = axial coordinates of wake 
IWK * wake index number 

Theory See reference 1. 

Function SBFUNC (X) 

Object Calculate value of special function. 

Argument List 


X 28 specified independent parameter 

SBFUNC « function value for specified X 

Theory A special function for the stall bucket (reference 11) is 

tabulated for interpolation on the independent parameter X. 
Beyond the range of tabulated data, the following function is 
used. 

SBFUNC - e X -2.139 
Subroutine SETMAT 

Object Read from a disk in a specified order the geometric influence 

coefficients and combine them on another disk for the matrix 
solution algorithm. 

List of Symbols 

NPR0P * number of propellers 

LT0T * number of propeller positions 

IT0T = number of inflow stations 

MSIZE ■ matrix row size 

GCC ■ chordwise influence coefficients 

GCN « normalwise influence coefficients 

GCS “ spanwise influence coefficients 


88 


Subroutine S0LVEL 


Object Calculate required quantities associated with the propeller 

aerodynamics and control the linearized solution procedure. A 
flow diagram is presented in figure 14. 

Options IDEBUG * 0 print out intermediate quantities 

NEVARD * 0 use Evvard Tip Relief Correction 

ICAS * 0 use cascade airfoil data inboard of requested radial 
station 

1PRMAT * 0 print matrix related quantities 
MATS0L = 0 use direct matrix solution technique 
MATS0L ■ 1 use iteration matrix solution technique 
INPT * 0 print circulation solution 

Theory See section entitled: "Linearized Aerodynamics" of reference 1. 

Subroutine S0LVEN 

Object Calculate required quantities for propeller aerodynamics and 

control the nonlinear matrix solution procedure. A flow diagram 
is presented in figure 15. 

Options ICAS * 0 use cascade airfoil data inboard of requested radial 

solution 

1PNT * 0 print intermediate circulation solutions 
MATS0L = 0 use direct matrix solution technique 
MATS0L * 1 use iterative matrix solution technique 
IDEBUG * 0 print out intermediate quantities 

Theory See section entitled: "Nonlinear Aerodynamics" of reference 1. 


I 
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Subroutine S0LVIT 


Object Control flow of solution procedure* 

Options I TYPES = 0 linear aerodynamic solution only 

1TYPES = 1 nonlinear aerodynamic solution 

Subroutine SPLIN3 (Argument List) 

Object Interpolate using spline fit. 

Options NW0T = controls calculation of derivatives 

Argument List 

XDATA = independent interpolation variable vector 
YDATA * dependent interpolation variable vector 
NDATA = number of interpolation table data points 

XIN - vector of requested interpolation table data points 

Y0UT = interpolation values at the interpolation points 
YPRIME = vector of derivatives at the interpolation points 

NXY = number of requested interpolation points 

NW0T * option control on derivatives 

Theory Standard spline fitting technique* 

Subroutine STARC (Argument List) 

Object Convert design lift coefficient to equivalent camber angle. 

Argument List 

DECL = design lift coefficient 
THSTAR = effective camber angle 

Theory Using second order interpolation technique on a prestored table 

of data, the design lift coefficient is converted to the 
equivalent camber angle for a double circular arc airfoil. 
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Subroutine ST0RE (l,J,K,X,Y) 


Object 

Argument List 


Object 

Argument List 


Theory 


Store the vector X into the vector Y where X and Y can be 
externally dimensioned as three dimensional arrays. 


1 , J ,K “ dimension limits of the X and Y arrays 
X * input vector (array) 

Y * output vector (array) 


Function SWPC0R (S ,C0 ,MACH,LAMLE ,LAMTE ,Y ) 

Calculate tip loss factor for propeller blade using conical flow 
theory. 


S = semispan of wing which is used to approximate swept 

propeller tip 

C0 ■ midspan chord of wing 
MACH “ Mach number 
LAMLE = leading edge sweep angle 
LAMTE * trailing edge sweep angle 

Y = spanwise position on wing measured from midspan 

SWPC0R a scaling result 

Concial flow theory for a thin sweep wing with subsonic leading 
edge and supersonic trailing edge is used to obtain the three- 
dimensional section C^. This solution is divided by the 
equivalent thin wing two-dimensional solution to obtain a 
scaling function to apply to actual two-dimensional tabulated 
C^ data . 
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Subroutine THITER 


Object Control selection of C p or Gj, blade angle iteration. 

Subroutine TITER (N ,T ,DQDT ,QWANT .QCALC ,T0L , IQ0K) 

Object Linearly interpolate or extrapolate on input vector to obtain 

required T at requested Q. 


Argument List 


N * iteration index 
T = output vector 

DQDT * initially assumed slope of Q versus T curve 

QWANT * requested Q 

QCALC * input Q vector 

T0L ■ tolerance or solution 

1Q0K = iteration control flag 


Subroutine UNBAR (Argument List) 
Object Bivariant Interpolation on data vector. 

Argument List 


T = table of bivariant interpolation data (Z vs. X and Y) 

IK = starting location of data 

XIN * requested interpolation point (X) 

YIN * requested interpolation point (Y) 

ZZ * interpolated value 
KK * interpolation flag 

Theory Use standard bivariant interpolation algorithm with degree 

choice internally coded. 
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Subroutine UNINT (Argument List) 


Object Univariant interpolation. 

Argument List 


NW * print output file number 
N * number of interpolation data points 
XA = vector of independent interpolation data 
YA * vector of dependent interpolation data 
X * requested interpolation point 
Y = resulting interpolated data values 
Z 85 interpolation flag 

Theory Interpolates over a four point interval using a variation of a 

2nd degree interpolation to produce a continuity of slope 
between adjacent intervals. 


Subroutine VECTOR 

Object Compute velocity related quantities required for propeller aero- 

dynamics. 

Options IDEBUG * 0 printout of selected intermediate quantities 

Theory Using vector algebra, the direction cosines of the local 

velocity vectors at the blade are calculated neglecting induced 
velocity terms, along with other related quantities. 

Subroutine VVECTR 

Object Compute velocity related quantities for the propeller aero- 

dynamics . 

Options IDEBUG * 0 printout intermediate quantities. 

Theory Using vector algebra, the direction cosines of the local 

velocity vectors at the blade are compiled including effects 
along with other related velocity quantities. 
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Subroutine WAKM0D (KT0T , MT0T , IDEBUG ) 


Object 


Read wake displacement corrections from disk and modify the 
internally calculated wake geometry. 


Options IDEBUG * 0 printout selected quantities 

List of Symbols 

PHI = wake azimuth position 

RW = radial wake coordinate 

ZW = axial wake coordinate 

KT0T = number of wake filaments per blade 

NT0T = number of wake revolutions 

JT0T = number of wake azimuth positions per wake revolution 
JT0T1 = JT0T+1 
IDEBUG = print option 


Theory 


See section entitled: Nacelle Influence on Wake Geometry of 

reference 1. 


Object 


Subroutine WRITGC (Argument List) 

Convert vectors of geometric influence coefficients to the blade 
coordinate system in the order in which they are calculated and 
write them to disk for later retrieval. 


Options 
Argument List 


IPRMAT * 0 printout geometric influence coefficients 


Theory 


IX = local blade station index specifying the particular 
blade element station to transform the influence 
coefficients vectors from cylindrical to blade element 
coordinate system 
LL « propeller position index 

IWK = propeller wake index 

IBB * blade index 

IP = propeller index 

IWRITE « disk number to write influence coefficients onto 

Standard geometric transformation applied to the cylindrical 
coordinate system influence coefficients to transform them to 
blade element coordinate system, before storing the disk for 
later retrieval. 



Subroutine ZER0GC (MK ,N2 ,N3 ,GC) 


Object Set an externally dimensioned three dimensional array to zero. 

Argument List 


MK,N2,N3 * external dimensions of array GC 
GC * input array 


Subroutine ZER0AL (THSTAR,APZL) 

Object Calculate angle of zero lift for isolated airfoils. 


Argument List 

THSTAR = effective camber angle 
APZL c angle of zero lift 

Theory The angle of zero lift is calculated for the requested isolated 

airfoil type by computing the linear lift curve slope near zero 
angle of attack and solving for the intercept of the straight 
line with the resultant linear lift curve slope. 
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Labeled Common Blocks used in the Propeller Portion 


Included herein is a list of the labeled common blocks in alphabetical 
order used in the propeller portion of the analysis and a description of each 
variable used in them (NDN designates a non-dimensional number). 

Variable 

Common Block Name (Object) Names Description of Variables 

AERDAT (Store 

Miscellaneous Aerodynamic 

Quantities) 


AA 

Linearized Lift curve slope 
(per radian) 


DIAG 

Diagonal element vector of 
geometric influence 
coefficient matrix (per ft.) 


C0NST 

Constant vector of circulation 
solution matrix (ft 2 /sec) 


CMACH 

Local total Mach number (NDN) 


SMACH 

Local section Mach number (NDN, 
normal to lifting line) 


CFDP 

Constant vector correction term 
( ft 2 / sec) 


SKEW 

Aerodynamic skew angle (degrees) 

AFDATX (Store 

Isolated Airfoil Data Package Quantities) 


I 

Inflow station index 


IFL 

Airfoil type flag 


IDL 

C L , C D calculation flag 


ICASDE 

Cascade correction flag 


ALPHA 

Local angle of attack (degrees) 


THET 

Local blade angles (degrees) 


TAUB 

Chord-to-gap ratio (NDN) 


ZM 

Local Mach number (NDN) 


DECL 

Design lift coefficient (NDN) 


H0B 

Thickness-to-chord ratio (NDN) 


ZMCR0M 

Critical Mach number (NDN) 


CL2 

Lift coefficient (NDN) 


CL3 

Temporary lift coefficient (NDN) 


CD 

Drag coefficient (NDN) 


DCDCL 

Change in drag coefficient with 
lift coefficient (NDN) 
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Common Block Name 


Variable 

Names 


Description of Variables 


BIGMAT (Large Storage Matrix, used primarily for wake geometry and matrix 
coefficient storage) 


PHI 

C0SPHI 

SINPHI 

RW 

ZW 


Wake segment azimuth angles 
(degrees) 

Cosine of wake azimuth angles 
Sine of wake azimuth angles 
Radial coordinate of wake 
segments (NDN) 

Axial coordinate of wake 
segments (NDN) 


CASDT (Store Cascade and Airfoil Related Data) 


SIGMAX 

THETAG 

THETAB 

i 

TAUB 

CDPER (Store Nacelle Drag Quantities) 

DFR 

DPR 

CC0M (Store Section Chord Data) 

CH0RD 

ALCRAD 

I 

i 

! ALCPHI 

ALCAXL 


Section cascade solidity, gap- 
to-chord ratio (NDN) 

Geometric angle between section 
chordwise vector and rotation 
plane (degrees) 

Section angle of attack 
neglecting induced terms 
(degrees) 

Section chord-to-gap ratio (NDN) 


Nacelle skin friction drag (lb) 
Nacelle pressure drag (lb) 


Blade element section chord 
length (feet) 

Blade element section chord 

length radial direction cosine 
(NDN) 

Blade element section chord 
length tangential direction 
cosine (NDN) 

Blade element section chord 
length axial direction cosine 
(NDN) 


i 


i 
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Common Block Name 


Variable 

Names 


Description of Variables 


CINC0M (Store Input Chord Data) 

CINPUT Input chord length (feet) 

CLCDDT (Store Circulation Solution and Airfoil Characteristics) 


CLSAV 

CDSAV 

ALPHA 

PHINSD 


CIRC 

SAVGIR 

FT RAN 

C0NSTI (Store Input Data) 

RPM 

S0UND 

DENSTY 

VIM0M 

BL 

R 

STN 

THETAO 

HUBQ 

DPSI 

REV 

CPI 

T0L 

CNSECT 

SC0 

RADCAS 


Section lift coefficient (NDN) 
Section drag coefficient (NDN) 
Section angle of attack 
(degrees) 

Section inflow angle (degrees) 
Section minimum drag coefficient 
(NDN) 

Section circulation (ft 2 /sec) 
Section circulation from 

previous iteration (ft 2 /sec) 
Interpolation function (NDN) 


Propeller rotational velocity 
(rrap) 

Freestream speed of sound (fps) 
Freestream density (slugs/ft 3 ) 
Momentum induced velocity (fps) 
Number of propeller blades per 
propeller 

Blade radius (feet) 

Number of inflow stations 
Blade angle (degrees) 

Hub torque (ft-lb^) 

Blade azimuth increment 
(degrees) 

Number of wake revolutions 
Requested power coefficient 
(NDN) 

Matrix solution tolerance (NDN) 
Fraction of blade chord measured 
from leading edge (NDN) 

Not used 

Blade radius denoting the end of 
the cascade region (NDN) 
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Variable 

Common Block Name Names 


V0RC0R 

DCPDT 

STACK 


CTI 

DCTDT 

ZHUB 

VKTAS 

TIPM 

RDTRAN 

RPMREF 

DPSIB 

OMEGA 

DTIME 


C0NST1 (Store Internal Constants) 

PI 

RC 

R4PI 

0N04PI 

C0NST2 (Store Miscellaneous Quantities) 

DIA 

0MGR 

ZMSQ 

UAX 

ZJI 

MU 


Description of Variables 

Vortex core (NDN) 

Change in power coefficient with 
blade angle (per degree) 
Position of lifting line as a 
fraction of the blade chord 
measured from the leading edge 
(NDN) 

Requested thrust coefficient 
(NDN) 

Change in thrust coefficient 
with blade angle (per degree) 
Coaxial hub displacement (NDN) 
Freestream velocity (knots) 

Tip Mach number (NDN) 
Interpolation radius limit (NDN) 
Reference rpm for steady load 
induced twist (rpm) 

Blade spacing azimuth interval 
(degrees) 

Propeller rotational speed 
(rad/sec) 

Periodic blade/wake geometry 
time interval (sec) 


ir 

ir/180 (radial degree) 

4ir x blade radius (feet) 
1/R4PI (per ft.) 


Propeller diameter (feet) 
Propeller tip speed (fps) 
Freestream Mach number squared 
(NDN) 

Freestream velocity (fps) 
Freestream advance ratio (NDN) 
Freestream velocity/tip speed 
(NDN) 
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Variable 

Common Block Name Names Description of Variables 

C0NST3 (Store Input Option Scalars) 

PR0PMN See description for scalar 

inputs 

PRMAT " 

PR0PT " 

DEBUG " 

PCHPLT " 

WAKE0P " 

WAKNAC " 

C0MPRS " 

EVAARD " 

SKIN0P " 

CAS CAD " 

TYP CAS " 

CBWAKE " 

C0FL0W " 

TAUEXP " 

CPTHET (Store Performance Related Quantities) 

NCP Performance iteration counter 

ICP0K Performance iteration control 

flag 

THETO Blade angle storage vector 

(degrees) 

CPCALC Power coefficient storage vector 

(NDN) 

CTCALC Thrust coefficient storage 

vector (NDN) 

FC0M (Store Force Data) 

FT0T Total force per unit area 

(lb f /ft 2 ) 

ALFRAD Total force per unit area radial 

direction cosine (NDN) 

ALFPHI Total force per unit area 

tangential direction cosine 
(NDN) 

ALFAXL Total force per unit area axial 

direction cosine (NDN) 
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Common Block Name 


Variable 

Names 


Description of Variables 


FLT0T Lift force per unit area 

( lb f /ft 2 ) 

ALFLRD Lift force per unit area radial 

direction cosine (NDN) 

ALFLPH Lift force per unit area 

tangential direction cosine 
(NDN) 

ALFAX Lift force per unit area axial 

direction cosine (NDN) 

FDT0T Drag force per unit area 

(H ef /ft 2 ) 

ALFDRD Drag force per unit area radial 

direction cosine (NDN) 

ALFDPH Drag force per unit area 

tangential direction cosine 
(NDN) 

ALFDAX Drag force per unit area axial 

direction cosine (NDN) 

FLIGHT (Store Freestream Quantities) 

ZMO Local station rotational Mach 

number (NDN) 

ZJO Local station advance ratio 

(NDN) 

FL0WDT (Store Inflow Distribution Data) 

DENS Section density ratio (NDN) 

S0UN Section speed of sound ratio 

(NDN) 

V0NVO Section axial inflow velocity 

ratio (NDN) 

UR0NVO Section radial inflow velocity 

ratio (NDN) 

VZER0 Section center axial inflow 

velocity (fps) 

VZER0B Section boundary axial inflow 

velocity (fps) 

GCDIMD (Store Geometric Influence Coefficients) 

GCDIM Radial, tangential and axial 

geometric influence 
coefficients (per ft) 
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Common Block Name 


Variable 

Names 


Description of Variables 


GCKDAT (Store individual trailing filament and bound vortex influence 
coefficients ) 


GCKRTZ Radial, tangential and axial 

trailing vortex filament 
influence coefficients 

GCRTZB Radial, tangential and axial 

bound vortex influence 
coefficients 


GE0DAT (Store Blade Geometry Quantities) 

RSB 

ZSB 

YSB 

RSBB 

ZSBB 

YSBB 

PHIBB 

RSC 

RSCC 

ZSCC 

YSCC 

PHICC 

XSBB 

XSCC 

COSLB 

SINLB 


Input x-wise segment boundary 
coordinate (NDN) 

Input axial segment boundary 
coordinate (NDN) 

Input y-wise segment boundary 
coordinate (NDN) 

Segment boundary radius (NDN) 
Segment boundary droop (axial) 
displacement (NDN) 

Segment boundary lag 
displacement (NDN) 

Segment boundary lag angle 
(degrees) 

Input segment center radial 
coordinate (NDN) 

Input segment center radial 
coordinate (NDN) after blade 
angle rotation 
Input segment center axial 

displacement (NDN) after blade 
angle rotation 
Input segment center lag 

displacement (NDN) after blade 
angle rotation 

Input segment center lag angle 
(degrees) 

Segment boundary x-wise location 
(NDN) 

Segment center x-wise location 
(NDN) 

Cosine of blade segment boundary 
lag angle 

Sine of blade segment boundary 
lag angle 
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Variable 

Common Block Name Names Description of Variables 


GE0INP (Store Input Interpolation Arrays) 


V0NV0X 

CX 

DTHETX 

T0VERX 

BETAB 

BETAC 

GE0M0D (Store Design Lift Coefficient) 

DESCLP 

GE00UT (Store Blade Section Properties) 

DTHETA 

AF0IL 

DESCL 

T0VERC 

B0D 


Axial inflow ratio distribution 
(NDN) 

Input chord distribution (feet) 

Input pitch angle distribution 
(degrees) 

Input thickness-to-chord ratio 
distribution (NDN) 

Load induced twist increment at 
blade segment boundary points 
(degrees) 

Load induced twist increment at 
blade segment boundary points 
(degrees) 


Section design lift coefficient 
(NDN) 


Input section pitch angle 
(degrees) 

Section airfoil type (NDN) 

Input section design lift 
coefficients (NDN) 

Section thickness-to-chord ratio 
(NDN) 

Section chord-to-diameter ratio 
(NDN) 
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Variable 

Common Block Name Names Description of Variables 

G400DT (Store Quantities from Aeroelastic Response Analysis) 


NG400 

NSEG 

XYZCG 


G400CL 


G400CD 


I0UNIT (Store Standard Input/Output Unit 

NREAD 

NWRITE 

NPUNCH 

INTDT1 (Store Indexing Limits) 

IT0T 

KT0T 

JT0T 

NREV 

NBL00P 

LT0T 

NPR0P 

HT0T 

MSIZE 

NBCALC 

IFL 


Option flag 

Number of blade stations in the 
response analysis 
Segment center coordinates of 
blade stations used in the 
response analysis 
Quadratic coefficients for 
segment from the response 
analysis 

Quadratic coefficients for 
segment Cq from the response 
analysis 

Numbers) 

Standard input unit 
Standard print unit 
Standard punch unit 


Number of inflow stations 
Number of inflow station 
boundaries 

Number of blade azimuth stations 
Number of revolutions of wake 
geometry 

Number of blade loops 
Number of propeller positions 
Number of propellers 
Number of segment endpoints 
LT0T*JT0T*JT0T 

Number of blade calculations 
Airfoil type index vector 
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Variable 

Common Block Name Names Description of Variables 

INTDT2 (Store Program Option Flags) 


NC0MPR 

Wake compressibility flag 

NEVARD 

Evvard Tip Relief flag 

IPRMAT 

Matrix print flag 

IPR0PT 

Wake geometry print flag 

NCFL0W 

Section Mach number test flag 

IWAK0P 

Wake model flag 

NACWAK 

Nacelle wake correction flag 

IV0RT 

Vortex core model flag 

ISKIN 

Skewed flow skin friction drag 
addition flag 

ICASDE 

Analytical cascade correction 
flag 

I CAS 

Cascade airfoil flag 

IDEBUG 

Intermediate print flag 

NCBWAK 

Compressible bound wake flag 

IPCH 

Card punch option flag 

ITYPCS 

Cascade option flag 

I BLN 

Blade number control flag 


1UNITD (Store Disc Unit Numbers) 


1UNIT 

MUNIT 

MC0NED (Store Tip Mach Cone Coordinates) 

XMC 

YMC 

ZMC 

XMCT 

YMCT 

ZMCT 


Disc unit number storage vector 
for geometric influence 
coefficients 

Disc unit number for matrix 
solution coefficients 


Input tip Mach cone definition 
coordinate x (NDN) 

Input tip Mach cone definition 
coordinate y (NDN) 

Input tip Mach cone definition 
coordinate z (NDN) 

Tip Mach cone definition 
coordinate x (NDN) 

Tip Mach cone definition 
coordinate y (NDN) 

Tip Mach cone definition 
coordinate z (NDN) 
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Common Block Name 


Variable 

Names 


Description of Variables 


MHC0NE (Store Mach Cone Correction Quantities) 


NSTAT 

XKC0NE 

N0RC0M (Store Section Normal Data) 

ALNRAD 

ALNPH1 

ALNAXL 


Inflow station index, where tip 
Mach cone intersects specified 
fraction of chord line 
Evaard Tip Relief correction 
factor (NDN) 


Blade element section normal 
radial direction cosine (NDN) 
Blade element section normal 
tangential direction cosine 
(NDN) 

Blade element section normal 
axial direction cosine (NDN) 


PH1C0M (Store Nacelle Induced Wake Azimuthal Distortion Quantities) 


NPHI 

DELPHI 

COSDPH 

SINDPH 

R0LL (Storage of Wake Rollup Quantities) 

TRUNCT 

TRUNCI 

R0LLUP 


Number of azimuth increments 
affected by nacelle 
Incremental azimuthal distortion 
(degrees) 

Cosine of distortion angle 
Sine of distortion angle 


Tip filament rollup truncation 
angle (degrees) 

Inboard filaments truncation 
angle (degrees) 

Number of filaments for rollup 
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Variable 

Common Block Name Names 

STAC0M (Store Section Span Data) 

STABAR 

ALSRAD 

ALSPH1 

ALSAXL 

THICKD (Store Thickness Data) 

THK 

U1C0M (Store Induced Velocity) 

UIR 

U1T 

U1Z 


Description of Variables 


Blade element section length 
(NDN) 

Blade element section length 
radial direction cosine (NDN) 

Blade element section length 
tangential direction cosine 
(NDN) 

Blade element section length 
axial direction cosine (NDN) 


Blade element thickness (NDN) 


Radial induced velocity (fps) 
Tangential induced velocity 
( fps) 

Axial induced velocity (fps) 


UUC0M (Store Input Noninduced Velocity Data) 

UR Radial noninduced velocity (fps) 

UT Tangential noninduced velocity 

(fps) 

UZ Axial noninduced velocity (fps) 
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Variable 

Common Block Name Names 

VC0M (Store Noninduced Velocity Data) 

ALPHAN 

VS 

VC 

VN 

VT0T 

ALVRAD 

ALVPHI 

ALVAXL 

VIDAT (Store Induced Velocity) 

VIS 

VIC 

VIN 

WAKDAT (Store Wake Rollup Data) 

KTRUCT 

JTRUCT 

JTRUCI 


Description of Variables 


Noninduced angle of attack 
(radians) 

Spanwise noninduced velocity 
direction cosine (NDN) 
Chordwise noninduced velocity 
direction cosine (NDN) 
Normalwise noninduced velocity 
direction cosine (NDN) 

Total noninduced velocity (fps) 
Radial noninduced velocity 
direction cosine (NDN) 
Tangential noninduced velocity 
direction cosine (NDN) 

Axial noninduced velocity 
direction cosine (NDN) 


Spanwise induced velocity (fps) 
Chordwise induced velocity (fps) 
Normalwise induced velocity 
(fps) 


Filament station index for 
rollup 

Wake azimuth position index for 
tip rollup 

Wake azimuth position index for 
root rollup 
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Nacelle Program 


A detailed description of the nacelle portion of the computer program is 
given in this section. The subroutines and external functions are described 
individually in alphabetical order. The labeled common blocks are briefly 
described along with the FORTRAN variables used in them. Flow charts and 
figures are provided whenever necessary to understand the objectives and 
theory for the subroutines or external functions. 


List of Subroutines and External Functions 


Name 

ALTMN 

AMF 

AMFL0 

AMU 

BATCH 

BILINE 

BLDGEO 

BLD0UT 

BLKDAT 

BLKRED 

BLPARM 

BPLUSR 

CALDRM 

CALINV 

CDS 

CKINPT 

C00R 

C00RST 

C00R1 

C00R3 

C00R4 

C00R3 

CPLX1 

DAMU 


Object 

Control I/O and calculation flow 
Compute isentropic nozzle flow 
Calculate Mach number from area ratio 
Compute molecular viscosity 
Main routine 

Calculate neighboring points on output line 

Locate blade centerline 

Store blade parameters on drum 

Load block data 

Read data records from mass storage device 

Compute boundary layer parameters 

Compute law of wall integration constant 

Read inviscid or viscous solution from drum 

Calculate inviscid flow field 

Calculate Roberts' mesh distortion parameter 

Check input data for radial equilibrium 

Interpolate coordinates 

Control flow of coordinate calculation 

Compute approximate coordinates 

Compute coordinate functions 

Compute Schwartz-Christof fel parameters 

Interpolate wall curvature at station 5 

Evaluate Schwartz-Christof fel transform 

Find derivati r e of molecular viscosity 


no 



Name 


DR0BRT 

DR0UT 

DRUTPE 

ERPIN 

FAMACH 

FAVER2 

FC0LES 

FC0RCT 

FCPLX 

FETA 

FINTG 

FL0WIN 

FN0RM 

F0RCE 

F0RCL 

FTHIK 

G BLADE 

GDUCT 

GE0MCL 

INITQ 

INTFRE 

L0ADRR 

MINVRT 

MYT1ME 


Object 

Compute derivative of Roberts' transformation 
Drum I/O routine 

Transfer data drum to tape 

Check normal pressure gradient 

Calculate Mach number from velocity 

Compute mean flow 

Compute Coles velocity profile 

Correct truncation error 

Evaluate complex functions 

Calculate distorted mesh 

Integrate complex functions 

Set inlet flow 

Normalize input variables 

Compute blade forces 

Compute local blade force 

Compute blade thickness 

Compute blade geometry 

Compute duct shape 

Calculate coordinates of lifting line 

Initialize data file parameters for Q array 

Initialize freestream conditions 

Loader formatted input 

Inverts block matrix 

Dummy time trap 
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Name 


0UTPUT 

PERFNA 

PERFN2 

P01S 

P0ISCF 

P0IS0N 

QINTER 

READPF 

READPG 

R0BRTS 

R0UND 

SCURVA 

SLETE 

SM00TH 

S0LVI 

SPLIN3 

STRESI 

STRT 

TPRINT 

TURB 

UBLAS 

UC0LES 

WAKC0R 


Object 

Print title page 
Compute viscous nacelle drag 
Compute inviscid nacelle drag 
Solve Poisson equation 

Set initial quantities in solution procedure 
Calculate axisymmetric streamline curvature 
Interpolate curvature 
Read P and F files 

Read variable for curvature calculation 

Compute distorted mesh using Roberts' transformation 

Round corners on straight wall ducts 

Calculate curvature from potential flow solution 

Find blade control surfaces 

Smooth duct wall contour 

Integrate equations of state 
NASA Spline Fit routine 

Compute initial stress distributions 
Find inlet flow locations 
Call CPU time 

Compute turbulent viscosity 

Calculate velocity ratio according to Blasius 
solution 

Compute Coles friction velocity 
Compute nacelle wake corrections 
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Name 


WBLEED 

WRITPF 

XH 

XT 


Object 

Calculate perforated wall bleed 
Store updated potential flow solution 
Calculate wall length on ID wall 
Calculate wall length on OD wall 


Description of Subroutines and External Functions 

This section describes the subroutines and external functions used in the 
nacelle portion of the analysis. The source name for the main control routine 
is called BATCH and has the following two entry points: ALTMN and (JFFLNE. 

The main program determines the entry point ALTMN or 0FFLNE. 
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Subroutine ALTMN 


Obj ect Controls I/O and calculation of flow. 
Options 

IREADR=0 Card formated input 
IREADR=1 Loader formated input 
All I0PT0 and IDBG0 options 

Entry Points 

ALTMN Initial loading of first case 


0FFLNE 

Not used 



List of Symbols 




AMACHE 

= M i 

9 

Average inlet Mach number (dimensionless) 

BBO 

- B o 

9 

Inlet blockage (dimensionless) 

DZ 

= AZ 

9 

Increment in axial length (ft) 

IREADR 

= 

9 

Loader format flag 

KDSH 

= KDS 

9 

Temporary storage for KDS 

PI 

" P 1 

9 

Average inlet static pressure (psf) 

REYH 

II 

25 

9 

Reynolds number based on duct height (dimen- 




sionless) 

T1 

= T i 

9 

Average inlet static temperature (deg R) 

Z 

= z 

9 

Axial length (ft) 


Theory 

This subroutine controls I/O and the calculation of flow depending on the 
options selected. A flow chart is shown in figure 16. 
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Subr ou t ine AMF ( AA , AMI , AMG ,AM,ACPC, ACPI) 


Obj ect Compute isentropic nozzle flow 
Opt ions 

AMG < 1 Subsonic solution 
AMG > 1 Supersonic solution 


List of 

Symbols 



AA 

* 

A/Ai 

Area ratio 

AAI 

= 

(A/ Al ) V , 

V t ^ 1 guess for area ratio 

ACPC 

S 

C P 

Pressure coefficient 

ACPI 

= 

C pl 

Incompressible pressure coefficient 

AM 


M , 

Mach number 

AMF 

= 

F(M) 

Function name 

AMG 

= 

m ( g) 

1st guess for Mach number 

AMI 

= 

M2 

V*”* 1 guess for Mach number 

AMI 

= 

Mi 

Mach number at station 1 

AP2 

= 

V P 1 

Static pressure at station 1 

AP2 

= 

P o /P 2 

Static pressure at station 2 

Cl 

= 

Y-l/2 

\ 

C2 

= 

1 Y+l 

2 

/m) 

C3 

= 

M i/<i +J£ r M i 

?\ 2 Y-l / \ 

) \ Constants in Iteration 


Constants in Iteration 



C6 

C7 

DAAI 

DAMI 

ITER 

Theory 


= y/y-1 


= Mi 


(1 


(1 


(iS) 

+ "T“ [ 

+ ^ J 


Constants in Iteration 


= (dA/dM) 
= AM 
= V 


Derivative 

Correction to Mach number 
Iteration counter 


Fiven the area ratio A/A^ and the inlet Mach number M^, find the exit 
Mach number M, 'the pressure coefficient Cp, and the incompressible coefficient 
Cpj. The exit Mach number M is determined from the one-dimensional isentropic 
flow relations using Newton’s method for determining roots of nonlinear equa- 
tions. Thus, we setup the iteration cycle 


( a^-.£(, + 4V,) + ^ 

/_A\ V r c 4 m>/ j_i 
“ U,/ l c 5 ~ J 


d(A/A,) 

V 

dM . 

\ A i ) 


AM = (A/A, - (A/A, ) v )/ -£-(A/A,) 

CIM 

+ AM 


( 1 ) 

( 2 ) 

(3) 

(A) 


When | AM | < 10 -5 the iteration has converged and the pressure coefficient may 
be computed. 


5l = ( 

r , + r :' 


pi ‘ 

i 2 '/ 

(5) 

O 

TJ 

•H 

II 


(6) 

Co = 

(P/P 0 -F|/P o ) 


p 

1 - p, /p 0 

(7) 
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Subroutine AMFL0 (AA, AMI, AMG, AM, ACPC, ACPI) 
Object Calculate Mach number from area ratio 


Variables 



AA 

A/A 

Area ratio 

ACPC 

C PC 

Compressible pressure coefficient 

ACPI 

Cpi 

Incompressible pressure coefficient 

AM 

M 

Mach number 

AMG 


Flag 

AMI 

Theory 

M 1 

Mach number at station 1 


The Mach number can be calculated from 


A 

A 


1 


m i/ 1 +I T m2 


Y + 1 

2(Y-1) 


M 


1 + 1 r M i 2 


using Newton's iteration. With M known, 



( 1 ) 


( 2 ) 


P 


T 


Pi 


[1 + ^r m i 2] 


Y 

Y-l 



AMG < Subsonic root of (1) 
AMG > Supersonic root of (1) 


(3) 

(4) 

(5) 
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Function AMU (T) 


Object Compute molecular viscosity 


Options 

None 


List of Symbols 

AMU = y/y r , Ratio of molecular viscosity (dimensionless) 

T =0 , Static temperature ratio (dimensionless) 

Theory 

The molecular viscosity is computed according to Sutherland's formula 
(Ref. 6). The working fluid is assumed to be air. Accordingly, 


Wr 


r © 3/2 1 + '98.0/ T r 

©+I98.0/T r 


Subroutine BLD0UT 

Object Write (or read) blade parameters on drum 

Entry Points BLDIN 
BLD0UT 
BLDIND 
BLD0UD 
BLDBK 


( 1 ) 


List of Symbols 

STRI , STR3, STR5, STR7 , ISTR8 blade parameters (dimensional) 

STR2, STR4, STR6 blade parameters (nondimens ional) 

Theory 

Blade parameters for each row of blades are stored on a drum. This 
routine will read or write onto that drum. 
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Main Program BATCH 


Object Main program 
Variables None 
Theory 

The main programs calls only subroutines TITLE and ALTMN 
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Subroutine BILINE (L ,RB1 ,ZB1 ,RB2 ,ZB2) 


Object 

Calculate neighboring points on output line 
Options 

None 

Variables 


L 

, Point 

number 

RBI ,ZB1 = Rj.Zj 

, Point 

at L-l 

RB2 , ZB2 = R 2 ,Z 2 

, Point 

at L 

Theory 




The points are read from an input table. 
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Subroutine BLDGE0 


Object 

Locate blade centerline in (n,s) coordinates 
Options 

I0PT2 =0 No blades in duct 

= 1 Blades in duct 

List of Symbols 
COMMON BL0CK Variables 
Theory 

This subroutine uses subroutine FLINE to find the centerline in the (n,s) 
coordinates. Thus the input data from the blade stacking plane is transformed 
to the duct plane using subroutine TRBLD and stored in the input data line data 
block BLNE (1,2) used by subroutine FLINE. The output (n,s) coordinates are then 
stored in the blade data array C0NST(I,L). At the completion of this calculation, 
the location of the upstream and downstream blade force calculation surfaces are 
determined by calling subroutine SLETE. 
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Subroutine BLKDAT 


Obj ect Load fixed constants 
Options 

None 


List of 

Symbols 


ACHI 

= 

X 

AKAPPA 

= 

K 

APLUS 

= 

A + 

CPR 

= 

C P r 

CVR 

= 

C 

vr 

EP 

= 

e 

GAMMA 


Y 

GASR 

= 

ft 

GRAVR 

= 

g 

PI 

= 

IT 

PRESR 


P r 

PRL 

= 

P rL 

PRT 

= 

P r t 

RH0R 

= 

Pr 

SNDR 

= 

C r 

TEMPR 

= 

T r 

TI 

= 

t 

VISCR 

= 

y r 


to program. 


0.016 (dimensionless) 

0.41 (dimensionless) 

26.0 (dimensionless) 

5997. (f t 2 /sec 2 /deg R) 

A283 . (ft 2 /sec 2 /deg R) 

2 . 7182818 (dimensionless) 
1.4 (dimensionless) 

1714. (/t 2 /sec 2 /deg R) 

32.2 (ft/sec 2 ) 

3.1415926 

2117. (ft 2 /sec 2 /deg R) 

0.72 (dimensionless) 

0.90 (dimensionless) 
0.00238 (slugs/ft 3 ) 

1116.0 (ft/sec) 

519.0 (deg R) 

0.01745329 (dimensionless) 
0.37 x 10" 6 (slug/ft/sec) 
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Values of parameter's defined in C0MM0N/PARAM/ are also included for 
IBM and CDC computer programs. 


Subroutine BLKRED (UNIT, RECSIZ, ADDR, BEGREC, NRECS) 


Object 

Reads NREC 'records' from file 'UNIT' beginning with record BEGREC. NRECS 
records are stored as a single block, beginning at ADDR. The Univac 1100 library 
I/O routine NTRAN is used in this subroutine. 

Variables 


UNIT 

logical unit# 

(Integer) 

RECSIZ - 

record size in words 

(Integer) 

BEGREC - 

first record to read 

(Integer) 

NRECS = 

# of logical records to read 

(Integer) 

ADDR 

beginning address to store the NRECS & 
records read 

RECSIZ 


Theory 

BLKRED (with entry BLKWRT) was developed to allow NTRAN compatibility with ANSI 
standard DEFINE FILE I/O operations. In particular, a call to BLKRED with NRECS = 1 
is identical to a random access fortran read. 

In order to simulate DEFINE FILE I/O, it is necessary for BLKRED to maintain 
a list of pointers into the various disk files. The pointer list, DSKLOC, is in a 
common block | UNITS | which must be allocated in a static (root) segment. The loca- 
tion pointer and read size are used to position the disk for I/O access. After the 
I/O access, the pointer is positioned accordingly. 

BLKRED will issue a diagnostic message and cause program termination if either 
of two abnormal conditions are detected. 

1) the record # is negative 

2) NTRAN returns on error status less than zero (see UNIVAC FORTRAN V 
library routine NTRAN description on UNIVAC ASCII FORTRAN routine 
NTRAN 0 description) 
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Subroutine BLP ARM 


Object Calculate boundary layer parameters from viscous solution. 

Options II * 1 Calculate hub boundary layer 
II = 2 Calculate tip boundary layer 


List of Symbols 


BLC(1, I) 

3= 

u 

oo 

BLD(2 , I) 

= 

n 

BLP (3 , I) 

= 

T o 

BLP (4 , I) 

32 

T 

BLP(5 , 1) 

32 

M 

BLP ( 6 , 1) 

= 

U s 

BLP (7, I) 

= 

P 

BLP (8, I) 

= 

Y 

BLP (9, I) 

= 

A*/A 1 

BLP(10, I) 

= 

e*/e 1 

BLP(11, I) 


A*/0* 

BLP (12, I) 

= 

N e 

CF 

S3 

c f 


Freestream velocity 
Static pressure 
Total temperature 
Static temperature 
Mach number 
Streamwise velocity 
Density 

Distance from hub 
Displacement thickness ratio 
Momentum thickness ratio 
Shape factor 
Reynolds number 
Wall Friction Coefficient 


Theory 

This subroutines 
thickness A* where 


determines the momentum thickness 0* and the displacement 

00 6o 

O 

e* «fj? (— - ^ ) dy 
J p u u 7 
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Function BPLUSR (AKPLUS) 


Object Compute the law of wall integration constant for initial profile, 


Options None 


None 


List of Symbols 


AKPLUS 


, Roughness Reynolds number 


BPLUSR 


B + (K + ) , Integration constant 


Theory 

The inner layer turbulence model given in FCOLES was integrated to get 
the constant of integration. A data correlation for B (K^ ) is then given by 


B + (0) = 2.2 


K*< 4.1270 


B + (Kg) = “ 0.81486 - 1.2070 * (in K s - 3.91538) , K j > 4.1270 


Equation (1) indicates that for K < 4.127, the smooth wall model applies, 

o 


Subroutine CALDRM 


Object Read inviscid or viscous solution from (drum) 


Options 


N0PT8 = 0 Read inviscid solution 
N0PT8 « 1 Read viscous solution 


List of Symbols 


Viscous flow parameters 


Inviscid flow parameters 
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Subroutine CALINV 


Obj ect Calculate inviscid flow field solution. 
Options 

Calculate flow from J = JFIRS, JLAS 
IF (I0PT15.NE.O) JFIRS=I0PT15 
IF (I0PT16.NE.O) JLAS=I0PT16 
Calculate only for I0PT1 = 3 or 4 
N0PT5 4 0 Error exit 


List of Symbols 




Same as CKINPT 




Variable store on drum 

K=1 ,KL 



BINP(1, K) 

P 

o 

9 

Total pressure (psf) 

BINP (2, K) 

P 

* 

Static pressure (psf) 

BINP (3, K) 

a 

9 

Swirl angle (deg) 

BINP (4, K) 

T 

o 

9 

Total temperature (deg R) 

Additional Variables 




ITERAL 

V 

a 

9 

Swirl angle iteration number 

ERRA 

E o 

9 

Local error in swirl angle 

ERRAM 

E nd 

9 

Maximum error in swirl angle 

Theory 





Given the swirl angle a, the analysis is identical to subroutine 
CKINPT. The calculation of the inviscid flow field requires also the solution 
of the angular momentum equation which is given by 


RU * aR i u *j (i) 

where RU^ is the angular momentum at an arbitrary station and R^U ^ is the in- 
let angular momentum which is given. An outer iteration loop is then pro- 
grammed to solve eqt. (1) to get the swirl angle a. With a known, the 
inner iteration loop is the same as subroutine CKINPT. This solution is 
obtained for each streamwise station J ■ JFIRS, JLAS. The computed solution 
BINP (4, KL) is stored on a drum. 
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Function CDS (DDS, DETA) 


Obj ect 

Calculate Roberts 1 mesh distortion parameter 
Options 

None 

Input Variables 

DDS * Ratio of mesh distortion at wall An/An 

DETA = An , Mesh size at boundary - uniform mesh 

An , Distorted mesh size at wall 

Output Variables 


CDS C 


Roberts* mesh parameter 


Theory 

Let C = 1/2 + e 

Then Roberts* transformation can be written 



(2A17-O 


4n = (!_+«>»-« 

i + 9 

An = A17/DDS 


( 1 ) 

( 2 ) 

(3) 

(4) 


Eq. (2) through (4) can be solved iteratively for e as follows: 

< = 0 


j An + e 
* = l+e-An 


(5) 

( 6 ) 




( 7 ) 
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Subroutine CDS (Cont'd) 


Convergence occurs when 


v + 1 v 
€ -f 


< IE-04 


( 8 ) 
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i 


Subroutine CKINPT 


Obj ect 

Check input 

data for radial equilibrium 

Options 




I0PT1 ^ A 

Do not calculate 


N0PT5 J 0 

Error exit 


I0PT5 = 2 

F0RCE data equals FL0WIN data 


IJ = 1 

Inlet flow data 


IB = 2 

Force data 


INEX = 1 

Upstream data 


INEX = 2 

Downstream data 


IDBG13 = 1 

Debug printout 


WFL0W = 0 

Static pressure check only 


WFL0V7 > 0 

Pressure check and weight flow iteration 


List of Symbols 


ERR 


e 

, Error in interaction 

EPS 


E o 

, Minimum error 

FG(2, 

K) 

<t> (2) 

, Equation 4 

FG(4 , 

K) 

*(2) 

f Equation 5 

ITER 


V 

, Iteration number 

PHMAX 



, Maximum pressure possible 

PSI1, 

PSI2 

V V 

Upper and lower bound air stream function 

WF 


«v 2 

, Weight flow v th iteration 

WFL0V7 


w 

, Input weight flow 

WMAX 


w 

max 

, Maximum weight flow possible 
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WMIN 

W min 

, Minimum weight flow possible 

XL 

X L 

, Lower bound on X 

XM 

X M 

, X for choked flow 

XU 

X U 

, Upper bound on X 

XI 

Y u Y u 
X 1 » 2 

, Iterative values for X 

PSIHT,PSIT 

u 

^ T . 'Pj 

, Value of stream function 

Theory 




Input data for the total pressure, static pressure, swirl angle, and 
total temperature must satisfy the continuity equation, and the radial momentum 
equation. If these equations are not satisfied, the static pressure is ad- 
justed. The solution of these equations can be obtained by a transformation 
of variables. Let 



(i) 


and 


0 ( 7 ?) « 2 


f i av 

(XV dn 


cos 2 a 


I 

XR 


dR . 2 \ 
drT sin a ( 


( 2 ) 


then the radial momentum equation becomes 


ax 

dt 


+ 


MU * ^ a7< lnn o) 


X = a(rj) 


(3) 


which is an ordinary first order ,, linear M equation. The solution is given by 

(4) 


*Oti=expj/o(Y)dY + Til 


•l 

♦ t 7 ?) -J a(Y) 4>{y) d y 


X = (X 0 +<J>(r?))/^( f7 ) 


(5) 

( 6 ) 
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where is the hub static pressure ratio. The continuity equation becomes 


dj£ _ n 0 G cos a v yT] 


d tj M r vy®^ 


X'-' (i-x) 


w0?) = 2 7r p r C r r r 2 g ^(TJ) 

The constant X Q is determined by the boundary condition 


W(l) = w 


using an iteration scheme described below. 


First let us examine the function 


f(X)= X ' (i-X) 


f (X) has a maximum at 


X * X 


"> r + 1 


Hence from equation (1) 


I!* 

n rt \r+i/ 


2 \ X-l 


Equation (12) is precisely the condition for choked flow when M = 1. Sub- 
stitution of equation (12) into equation (10) yields an a priori condition 
for the maximum weight flow possible (i.e., the choked flow condition). For 
a subsonic solution we then have the condition 


<X <l 


Furthermore, it is noted that we can find a priori X Q such that 


7+7 <X L < X 0 <X V < I 
f(X L ) < f(X M ) 
f(X 0 )>0 
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by substituting equation (6) into equation (10). Thus X^ and are the 
bounds for choosing subsonic solutions with no reverse flow. The iteration 
scheme then consists of narrowing the bounds of X L and x u until convergence 
occurs. This procedure is illustrated in figure 17. 


,V*I 


x’ = x; + 


*-*, v 


- x ;) 


(15) 


if(^<v) xr=x v j <=^ 

/ V (16) 

If(^ V+, »J/) X 2 V =X V j^2 = ^ 

v v+l 

and ip is obtained by integrating equation (7) with X Q = X substituted into 

equation (6). Convergence occurs when 



(17) 


Once X is known, the static pressure is obtained from (1) and substituted 
for the input static pressure. 
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Subroutine C00R(JS,KS) 


Object Controls logic in determining coordinates and interpolates 
coordinates streamwise. 

Options 

I0PT9=O Computes approximate coordinates 
I0PT9^O Coordinates stored on drum 

I0PT9=1 Compute exact coordinate functions and store on drum 
I0PT9=3 Read coordinates from drum 


List of 

Symbols 


DSTEP 

= 

AS 

, Streamwise step size (dimensionless) 

DX 

= 

AX 

, Interpolation between JS and JS-1 stations 

DXNEXT 


ax n 

, Interpolation for next step 

JDRUM 

= 


, Drum unit number 

JS 

= 


, JSth station stored on drum 

KS 

= 


, KS t ^ 1 station interpolated between JS 

ZHUB 

= 

Z H 

, Axial station hub (dimensionless) 

ZNEXT 

= 

Z N 

, Next axial station (dimensionless) 

ZTIP 

= 

Zt 

, Axial station Tip (dimensionless) 

Theory 





Let the streamwise coordinate S be given by 

s = AS(JS-|) + dS-(KS-l) 


( 1 ) 


where 


AS = S l AJLAST-|) 
dS = AS/KDS 


( 2 ) 

(3) 


Then if station 1 is at JS-1 and station 2 at JS, a simple linear interpolation 
of the coordinates may be made from those stored on the drum. ZNEXT is the 
axial location of the KS+1 station. 
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Subroutine C00F5 (XH, XT, CURVH, CURVT) 


Object 

Interpolate wall curvature at Station S 
Options 


None 




Variables 




XH 

x H ( s ) 

* 

Arc length ID wall (dimensionless) 

XT 

X T (S) 

9 

Arc length OD wall (dimensionless) 

CURVH 

V s) 

9 

Curvature ID wall (dimensionless) 

CURVT 

Kj(s) 

* 

Curvature OD wall (dimensionless) 

RM(3, J) 

V J) 

9 

Table of X^ (dimensionless) 

RM(4 , J) 

K HI (J) 

9 

Table of (dimensionless) 

RM(7,J) 

( J) 

9 

Table of Xj. (dimensionless) 

RM(8, J) 

Kjj ( J) 

9 

Table of ICp (dimensionless) 

Theory 





The values of X^S) and X T (S) are input. Then the tables are searched and the 
values of fC (S) , K T (S) are calculated by linear interpolation. 
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Subroutine C00RST 


Ojbect 

Options 

List of 

DSTEP 

DX 

DY1 

DY2 

ICK 

II 

ISLT1, 

ISLT2 

IW1.IW2 

JC0UNT 

KLH0LD 

KN 

MSL0T 

NSL0T 

RACl.J) 

X 


Controls flow of coordinate calculation 


None 

Symbols 


» AS = DS 
= An 
= <SY 
= <$ 2 Y 
= 0,1 
= 1,2 


, Streamwise step size (dimensionless) 

, Normal coordinate step size (dimensionless) 

, First derivative (dimensionless) 

, Second derivative (dimensionless) 

, Flag; no overlap, overlap 

, Flag; start integration, continue integration 


f First slot number, second slot number 


- R(I,1, J) 
= X 


, First slot wall, second slot wall 
, Streamwise station counter 
, Number of streamlines to interpolate 
, Number of streamlines to integrate 
, Slot counter 
9 Total number of slots 
, Temporary storage for wall coordinate 
, Normal coordinate (dimensionless) 
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XM 


X 


, Interpolation distance (dimensionless) 


X2 * x 2 » Midpoint of three point difference (dimensionless) 

Y = Y , Function to be interpolated (dimensionless) 

, Known values of Y (dimensionless) 

, Known values of Y (dimensionless) 

, Known values of Y (dimensionless) 

, Axial distance (dimensionless) 

, Location of adjacent slots (dimensionless) 


This subroutine controls the calculation flow for the coordinates 
according to flow chart figure 18. The basic calculation scheme with slots 
is to calculate the streamlines through successively larger ducts and 
storing only those coordinates satisfying the condition (Z^ZZ^Zy) as shown 
in figure 19. 

In addition, it was determined that only KN streamlines need be 
calculated by integration, the remainder up to KL streamlines may be 
calculated using a parabolic interpolation. 


Y1 

= Y(X r )=Y 1 

Y2 

= y(x 2 )=y 2 

Y3 

= y(x 3 )=y 3 

Z 

= Z 

zz L ,zz 2 

= z r z 2 

Theory 
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Subroutine C00R1 (KSS,JSS) 


Object Compute approximate coordinate functions. 

Theory 

This subroutine provides an initial guess from the iterative method 
described in reference 6 which is used to approximate the coordinate 
functions and hence determine the numerical grid structure used in the 
flow analysis. 

Subroutine C00R3(KSS ,JSS,L0P) 

Obj ect Compute coordinate functions (I0PT9=1 ,2) . 

Options 

L0P=1 Compute initial constants 
L0P=2 Integrate one step 
L0P=3 Do not integrate 


List of 

Symbols 



AO 

= 

A 0 

9 

Inlet height W plane (dimensionless) 

ALO 


ao 

9 

Inlet angle W plane (dimensionless) 

ANO 

= 

n o 

9 

Inlet height Z plane (dimensionless) 

SLO 

= 

b L 

9 

Initial coordinate length (dimensionless) 

RO 


r o 

9 

Inlet radius Z plane (dimensionless) 

VO 

= . 

v 0 

9 

Inlet metric scale coefficient (dimensionless) 

SO 

ss 

s o 

9 

Inlet streamwise coordinate (dimensionless) 

XBO 

= 

'k. 

A o 

9 

Real part of C-^ (dimensionless) 

YBO 

= 

Yo 

9 

Imaginary part of C^ (dimensionless) 

ZETAO 

= 

Ko 

9 

Constant of integration (real) (dimensionless) 

ETAO 

Theory 


ho 

9 

Constant of integration (imaginary) (dimensionless) 


Using procedures described in reference (6), one can integrate along 
streamlines to get R, Z, V, S, N at the location of the poles b^. 
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Subroutine C00R4 


Obj ect Find Schwartz-Christof f el parameters . 
Options 


None 





List of Symbols 




AX(1, J) 

ss 


9 

Distance along hub wall (dimensionless) 

AX(2 , J) 

« 

x r 

9 

Distance along tip wall (dimensionless) 

AX(3, J) 

= 

“H 

9 

Wall angle hub (deg) 

AX(4,J) 

= 


9 

Wall angle tip (deg) 

AY(1,I, J) 

= 


9 

Streamwise coordinate hub (dimensionless) 

AY(2,I,J) 



9 

Distance along hub wall (dimensionless) 

AY(3,I, J) 

= 


9 

Metric scale coefficient hub (dimensionless) 

AY(4,I,J) 

= 

*H 

9 

Radius hub (dimensionless) 

AY(5,I,J) 

= 

s^ 

b T 

9 

Streamwise coordinate tip (dimensionless) 

AY(6,I,J) 


Xj 

9 

Distance along tip wall (dimensionless) 

AYC7,1,J) 

= 


9 

Metric scale coefficient tip (dimensionless) 

AY(8,I, J) 

ss 


9 

Radius tip (dimensionless) 

AERR(1,J) 

as 

v+1 v , 

S H “Sh v 

9 

Error in S^ (dimensionless) 

AERR(2,J) 

= 

v+1 v 

9 

Error in S^ (dimensionless) 

AO 

= 

h o 

9 

Inlet height W plane (dimensionless) 

ALO 

= 

a o 

9 

Inlet angle W plane (dimensionless) 

ANO 

= 

n o 

9 

Inlet height Z plane (dimensionless) 

SLO 

= 

cO 

S L 

9 

Initial guess of coordinate length (dimensionless) 
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RO 

* 

r 

0 

9 

Inlet radius in Z plane (dimensionless) 

VO 

- 

V 0 

9 

Inlet metric scale coefficient (dimensionless) 

SO 

- 

S o 

9 

Inlet streamwise coordinate (dimensionless) 

XBO 

as 

X o 

9 

Real part of constant (dimensionless) 

YBO 

- 

'X/ 

Y 

o 

9 

Imaginary part of constant (dimensionless) 

ZETAO 

= 


9 

Constant of integration (real) (dimensionless) 

ETAO 

= 

% 

9 

Constant of integration (imaginary) (dimensionless) 

Theory 






The theory to this subroutine is described in reference (1) . 
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Subroutine CPLXl ( XI , Y 1 , XB 1 , YB 1 , N 1 , N2 , L0P ) 


Ojbect Evaluates Schwatz-Christof f el transform. 
Options 


L0P=1 W z - X + iy - dW/dZ 



=2 W zz = 

X + iy = d 2 W/dZ 2 


List of Symbols 



N1 

= N X 

, Product index 


N2 

" N 2 

, Product index 


N3 

= N]+l 

, Product index 


XI 

= X 

, X coordinate in Z plane 


Y1 

= Y 

, Y coordinate in Z plane 


XB1 

= X 

, X coordinate in W z plane or W zz 

plane 


a. 



XB2 

= Y 

, Y coordinate in plane or W^ z 

plane 

Theory 





This subroutine evaluates the real and imaginary parts of the Schwartz- 
Christoffel transform. 

Function DAHU 

Object - Find derivative of molecular viscosity with respect to temperature 
List of Symbols 

T Temperature 

Derivative 

Theory 

The derivative of the molecular viscosity may be expressed as a 
function of temperature. 
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D = A2 * |^1.5/T - 1/ (T+ Al) j (1) 

where 

Al = 198. /T r (2) 


A2 = T**l. 5 * (1+A1) / (T+Al) 


(3) 
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Function DR0BRT (C, ETA, L0P) 


Object 

Compute derivative of Roberts' transformation 
Options 

L0P = 0 wall - wall boundary 

L0P = 1 wall-freestream boundary 
I0P = -1 freestream-wall boundary 


Input Variables 
C = C 

ETA = n , 

L0P 

Output Variable 
DROBRT = 9n/3n 


Distortion parameter 
Input variable 
Option 

Output variable 


Theory 


The transform of the Roberts' Stretching for a distorted mesh is given by 


d T) 

Tn 



4 c In 


/ C + 1/2 \ ft 1-i 
\ C - 1/2 / (I +ft) 2 J 


(1) 


ft = exp 


[« 


in 


/ C-M/2 
V C - 1/2 


) (V - •/*)] 


where the options are 


( 2 ) 


n'.rf } L ^ P = 0 


V = i7/2 
n = 2n' 



r)' = (\ + r))/ 2 •» 

n = 2 o' - I 


We note that 


0 < rj < 1.0 


0 < n' < i.O 


(3) 
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Subroutine DR0UT (UNIT, ADDR, BL0CK) 

Obj ect Tape/Drum Read/Write data 
Options - Entry Points 
DR0UT (UNIT , ADDR , BL0CK) 

Write on drum unit UNIT data contained in ADDR containing BL0CK 
number of words . 

DRMIN (UNIT , ADDR , BL0CK) 

Read from drum unit UNIT data contained in ADDR containing BL0CK 
number of words. 

DREWND (UNIT) 

Rewind drum unit UNIT. 

DRMBK1 (UNIT, ADDR, BL0CK) 

Back space drum one BL0CK and read as in DRMIN. 

DRMBK2 (UNIT , ADDR , BL0CK) 

Back space drum two BL0CK and read as in DRMIN. 

DR1IE0UF(UNIT) 

Write end of file on unit UNIT (tape only) . 

TPFILE (UNIT , LFILE) 

Find file number LFILE on tape unit UNIT. 

Theory 


This subroutine is an I/O routine to facilitate conversion of 
UNIVAC NTRAN I/O routines to other computer systems. 
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Subroutine DRUTPE (LFILE , JBL0CK,L0PT) 


Object Transfers data Drum/Tape or Tape/Drum . 
Options 


L0PT=1 Transfer data tape to drum 
L0PT=2 Transfer data drum to tape 

List of Symbols 


JBL0CK 

JDRUM 

JTAPE 

LFILE 


, Number of records to copy 
, Drum unit number 
, Tape unit number 
, Tape file number 


Theory 


This subroutine transfers data using subroutine DR0UT. 


Subroutine ERPIN(II) 

Obj ect Checks normal pressure gradient for viscous flow 


Options 

11*1 Check inlet station 
11=2 Check J-l station 
11=3 Check J station 


List of Symbols 

ERPIN * £ , Error in normal pressure gradient 

Theory 

At any station in the duct, the normal momentum equation must be 
satisfied. This subroutine checks the error. 


E = 


[ni')-nio)l-(Hi -f£] pu « + [^-flr] pu »} J « 

[mo - nco)] 


( 1 ) 
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Function FAMACH (AN) 

Object Calculate Mach number from velocity 
Variables 

AN V/C 0 Velocity/stagnation speed of sound 

FAMACH M Mach number 

Theory 

The Mach number can be calculated from 

7- - M/(l + M 2 ) (D 

L o z 

using Newton’s method. 
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Subroutine FAVER2 


Object Solve for mass flow weighted average flow and output parameters. 
Options 

Ncne 


List of 

Symbols 



AAR 


= 

AR 

, Area ratio (dimensionless) 

AHA 


= 

M 

, Local Mach number (dimensionless) 

AMACHE 


* 

M 

, Area average Mach number (dimensionless) 

AMFH 


= 

(GM/V) h 

, Wall bleed hub (dimensionless) 

AMFT 


= 

(GM/V) t 

, Wall bleed tip (dimensionless) 

AMG 


= 

m mid 

, Midpoint Mach number (dimensionless) 

AMM 


= 

m max 

, Maximum Mach number (dimensionless) 

ASH 


= 

a sh 

, Surface area hub wall (dimensionless) 

AST 


= 

a st 

, Surface area tip wall (dimensionless) 

ATFH 



(GQ/V) h 

, Heat flux/length hub (dimensionless) 

ATFT 


- 

(GQ/V) t 

, Heat flux/length tip (dimensionless) 

CFH0 


= 

C FH J 

, Wall friction coefficient hub (dimensionless) 

CFT<t> 


= 

C FT J 

, Wall friction coefficient tip (dimensionless) 


* 

= 

1,2,3 for 

J+l-4> station 

CPC 


= 

Pr/q], 

, Normalizing factor for C p dimensionless 

CPC0MP 


= 

C PC 

, Pressure coefficient (compressible) (dimensionless) 

CPJNC 


= 

C PI 

, Pressure coefficient (incompressible) (dimensionless) 

DASH1,DASH2 

as 

m sh 

, Area increment hub (dimensionless) 
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DAST1,DAST2 

= 

AA S t 

, Area increment tip (dimensionless) 

DENTP1, 

r= 

AI 

, Change in entropy (dimensionless) 

DENTP2 




DISP 

= 

i 

, Dissipation function (dimensionless) 

DPSTI1, 

as 


, Increment in mass flow (dimensionless) 

DPSTI2 




DQSH1, 

- 

aq t 

, Increment in heat flow hub (dimensionless) 

DQSH2 


1 H 


DQST1, 

= 

a^t 

, Increment in heat flow tip (dimensionless) 

DQST2 




DTHE01 

S 

A0 O 

, Increment in total temperature (dimensionless) 

DTHE02 




PIO 

B 

n 

o 

, Average total pressure (dimensionless) 

PIOO 

- 

n 0 o 

, Average inlet total pressure (dimensionless) 

PSI 

« 

$ 

, Mass flow (dimensionless) 

PSIO 

- 

»0 

, Initial mass flow (dimensionless) 

QM 

= 

1/2 (.PU 2 )^ 

, Free stream dynamic pressure (dimensionless) 

QSH 

B 

'V 

Qh 

, Total heat flow hub (dimensionless) 

QST 

8 

q t 

, Total heat flow tip (dimensionless) 

R0M 

- 

Pmax 

, Free stream density (dimensionless) 

RUM 

- 

(pU )max 

, Maximum momentum (dimensionless) 

THEO 

B 

So 

, Average total temperature (dimensionless) 

THEOO 

- 

®oo 

, Inlet average total temperature (dimensionless) 

UM 

- 

U MAX 

, Free stream velocity (demens ionless) 

Theory 





The mass flow weighted average quantity <j> is defined by 
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( 1 ) 


* ■ +1 


r '^ ♦«" 


where 


* • / 


' GPU S 


dn 


( 2 ) 


The area is given by 


*/f* 


(3) 


It is noted that the mass flow weighted averages satisfy certain conditions 

dt . /GM\ 4 
dT ' \ V/T \ V /H 


£ (V®0> 
as 


[iiM, . [lit]. - j,. ipavm.1 <» 


(4) 


'6M® 


The entropy is related to the dissipation by 
h . /G&I\ /Gftl \ 


f -6_ 

J o ©V 2 

[^ns^ns 

♦ 

E n^ + ^R]dn 

r r 

/GQ_\ _ 

(GQ_\ 

, (6) 

y- 1 l 

\V@/ T 

V v© 'h 

1 1 


Equation (4) states that the change in mass flow is the net flow crossing the wall 
boundary. Equation (5) states that the change in total energy flux is the net 
crossing the boundary walls. Thus, for an adiabatic wall, the mass flow weighted 
average total temperature is constant. Finally, equation (6) states that 
the change in entropy flux is the net crossing the boundary plus the change 
due to the dissipation function and heat fluxes. Then using the definition 
of entropy we have _ y 

(7) 


t * 


>01 

and the loss coefficient is given by 


'PL 


Pq. ~ ^02 

P 0I - P, 


V 


I - p,/ p, 


01 


1 JWno, 

i-ff, /So, 


( 8 ) 
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If we now defined the mass flow weighted averages 



n* / 

7r M r 

r~> t \-2L 
— M ) y -' 

dn 

(?) 

©0 = °r / 

•to 

~ 77 -y@r m 

/ r _| 2 ' 

r 2 

}dn 

(10) 

*1 

It 

o 

[' 6 z 
J 0 v /© 

M dn 


(11) 


This subroutine also calculates some quantities printed in the summary at the 
end of the viscous calculation* 


The pressure coefficient is defined by 


n - n, 

n 0 , n i 


( 12 ) 


An effective area A can be defined as the geometrical area minus the block- 
age caused by the boundary layer. If we define the freestreain as the point; 
of maximum velocity across the duct we have by definition 


u - A eff = ™ * (* T “ n i 


( 13 ) 


Hence, the blockage B is defined as 


B = I - A/A eff = I - m/t^U.A) 


( 11 *) 


The area averaged (effective Mach number) may then be defined by the isentropic 
flow relations. Thus, 


A eff 

A 


M eff 

M 



I + 



1 111 
2 Y " 1 


( 15 ) 


The effectiveness n of the diffuser is based on an ideal isentropic flow 
with the mass flow weighted average Mach numbers. Thus, the ideal pressure 
coefficient is given by 


149 



06) 


C 


pi 


n - n, 


n 0 ,-n, 


where 


Then 




Mi 

M 



I + 






C P^ C PI 


(IT) 

( 18 ) 


( 19 ) 


The wall friction coefficient is defined as 


c, * W{z PU ‘) 


MAX 


(2D) 


The wall surface area and heat flow are determined by integrating the equations 

• s 


A S - i — dS 


( 21 ) 




( 22 ) 


using the trapezoid rule. 
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Function FC0LES (Argument List) 

Obj ect Compute Coles velocity profile for boundary layer at initial station. 

Options 

None 

List of Symbols 


AK 

= 

K 

> 

Von Karman constant 

APLS 

= 

A + 

> 

Van Driest constant 

API 

= 

M 

9 

Coles shape factor 

AKPL 

= 

V 

9 

Roughness Reynolds number 

DELT 

- 


9 

Boundary layer thickness 

DAMP 

= 

D 

9 

Damping factor 

PI 

= 

■n 

9 

3.14159 

Y 

= 

Y 

9 

Distance from wall 

YPLUS 

= 

Y + 

9 

Universal distance 

UPLUS 

= 

U + 

9 

Universal velocity 

FC0LES 

Theory 

s s 

U + 
c 

9 

Coles velocity (output) 


This subroutine integrates the differential equation for the inner layer 
given by 

dU*~ £ 

dy + l+[l + 4K 2 y +2 0 2 ] 1 ^ (1) 

with the damping factor, references (7, 8), given by 

o.,-e,p(-ii) + ( l + ^l)exp(-23^) (2) 
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and adds Coles' wake function (reference 3), given by 



Subroutine FC0RCT 


Obj ect Correct truncation error 
Options None 
List of Symbols 


AMUW 

=5 

pw/M 

9 

Wall value of viscosity (dimensionless) 

DU 

' = 

AU 

9 

Velocity difference (dimensionless) 

EC0R 

= 

( P u) / ( P u) 

m 

9 

Mass flux ratio (dimensionless 

EMB 

= 

E 

m 

9 

Error in Mach number (dimensionless) 

EPB 


E 

P 

9 

Error in static pressure (dimensionless) 

ERB 

= 

e r 

9 

Error in density (dimensionless) 

ETB 

ss 

e t 

9 

Error in static temperature (dimensionless) 

EUB 

= 

E u 

9 

Error in velocity (dimensionless) 

EUP 

= 

E U0 

9 

Error in swirl velocity (dimensionless) 

EUS 

= 

^S 

9 

Error in streamwise velocity (dimensionless) 

EW 

= 

E w 

9 

Strain 

Theory 






The point-to-point instability described in reference 6 is minimized by 
recalculating the stresses and heat flux using central differences rather 
than centered differences. 
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Subroutine FCPLX 


Object Evaluates complex functions for exact coordinate calculation. 


Options 



L0PT 

= 1 

Compute and store functions 

L0PT 

= 2 

Compute only derivatives 

List of Symbol 

s (Note subscript notation 

N1,N2 

= 

N 1’ N 2 

XS 

- 

S 

XN 

- 

n 

XSX 

- 

S = n 
x y 

XSY 

S 

S « n 
y x 

XD 

= 

D 

XZETS 


- n 
S n 

XETAS 


n = -£ 
s n 

XXS 

= 

x s " 

XYS 

- 

y = x 
J s n 

XB1 

= 

* - «* - "y 

YB1 

= 

% 

Y - -5 = n 
y x 

XI 

“ 

X 

Y1 


y 

XDB 

= 

D 

XNZET 

= 

n ? 

XNETA 

a 

n 

n 

XV 

a 

V 
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xsxx 


S = n “ s 

xx xx yy 

XSXY 

= 

S = n “ n 
xy yy xx 

xxss 

S5 

X ss “ ^ns - x nn 

XXSN 

= 

X sn ' y nn ' ' y ss 

XB2 

= 

X _ ^xx “ ^yx “ - ^xy 

YB2 

* 

X “ ~^xy “ -n yy "\x 

XZSS 

= 

^ss =n sn “ ~^nn 

XZSN 

S 

^sn _T1 nn n ss 

XDBN 

as 


XDBS 

**“ 

5 s 

XESDN 

88 

(n s /5) n 

XESDN 


<V 5) . 

XZSDN 

= 

<5 s /5 >„ 

XZSDS 

as 

« s / 5 >s 

XVN 

ss 

u n 

XVS 

ss 

V 

S 

XESS 

* 

’"’ss 

XESN 

as 

^sn 

XD1S 

= 

D s 

XDIN 

as 

D n 
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Theory 


The theory for evaluating the complex functions and all derivatives is 
derived in reference (6) . With the use of orthogonality relations which are 
implicit in the theory of complex function, the functions and derivatives may 
be evaluated. It is noted that this subroutine was programmed to accept 
multiple sources in the z plane, although only one is used in the present 
calculation. The derived functions calculated in this subroutine are listed 
as follows: 


NS A, r 

S = T ln [(* _ bl ) 2 + y z ] 

(1) 

NS 

n « X A, ton"' [ y/ (x — bj)] 

(2) 

S„ . £ A. <*->!> 

1,1 [ix-b/V] 

(3) 

NS At y 

S, = X — 

[lX-b I ) ! + y*] 

(A) 

s.. = r A ’ (, ) 

Ir| [ (X * b I )2+ y e ] ( [(X-b I )*+y 2 ]| 

( 5 ) 

c . A i 2y(x-b x ) 

^yy “ “ L 7 

(6) 

D = '(S x * + S y 2 ) 

( 7 ) 

X s = - S x /D 

(8) 
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Y S = “ Sy/D 

(9) 


£s = £x X s + £ y Y s 

(10) 


= Vx X S + Vy Y S 

(ID 


v = ! 

r r e r z V' 2 
Ks+Cn] 

(12) 

D s 8 

~ [2 Sy (Sy X X S + SyyYg) + 2 S X ( S XX X, + S X yY S ) j 

(13) 

D n 8 

“[2Sy(Sy X X n + SyyY r ) + 2 S R ( Sj J X + S g yYp ) j 

(14) 

*s$ 

r S xx X s + S xy Y s S x D s -I 

' L D D 2 J 

(15) 

x sn 

. r (S xx X n + S xv Y n ) S x D n I 

L D D 2 J 

(16) 

*SS 

£x X ss + ^xx x s + ^£ X y Y s *s + £y\s + £yy^is 

(17) 

e 8 n 

S C x X sn + Cy Y sn + (Cyy- £j V* ♦ ^y (Xj - Y # # ) 

(18) 

V s 

s -v 8 /[e $ Css^„ e„ s ] 

(19) 

V n 

” V /[C s £$n + C n Cnn] 

(20) 

Numerical accuracy can be significantly improved by ordering the way in which 
sums and products are made. As an example, the first equation, equation (1) , 
may be written 

NS 

s 8 1 

i *i 



NS 

= I 
1 = 1 

M'4*( x ; b, )T} , » ,> l x - b -l 
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Thus the square root of the sum of squares of 0(1) and S = 0(|x-bj|). This 
rule has been applied to all equations by extracting the order of magnitude of 
the term from each calculation. 


Subroutine FETA (B, ETA, AN, DEDN, D2EDN) 

Object 

Calculate distorted mesh to be used in Subroutine P0IS 

Options 

B = 0 
B > 0 

Variables 

B 

ETA n 

AN n 

DEDN 9n/9n 

D2EDN 9 2 n/9n 2 

Theory 

This subroutine calculates the distorted mesh that will be used in calculation 
by subroutine P0IS. The transformation is given by 

n ■ tanh [Bn] / tanh [B] B > 0 ) 

( 1 ) 

n = n B = 0 I 


Uniform mesh (no stretching) 

Tanh stretching 

, Constant 

, Transformed Normal Coordinate 
, Normal Coordinate 
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Subroutine FINTG (IKL) 

Object Integrate equations for potential flow. 
Options 

IKL = Number of streamlines 


List of Symbols 

IKL , Number of streamlines 

Theory 

Four simultaneous ordinary differential equations 
third order Runge-Kutta numerical integration method, 
given in FCPLX and denoted as 


dx 

dS 

dY 

dS 

dc 

dS 

d Tj 

"ds 


x s (x,y) 
Y s (x ,Y) 

e s (x ,y) 

1? S (X,Y) 


The Runge-Kutta formulas applied to the first equation 

B n = ASX S (X,Y) I 

B i2 s ASX s (x + B m /2 ,Y + B ei /2) ( 

B„ = ASX S (X+2B„-B„ .V + BB.j-B,,) ( 

x(S+AS) = X(S)+(B„+ 4B„+ B„)/6 


are integrated using a 
These equations are 


( 1 ) 


are 


( 2 ) 


(3) 


159 


Subroutine FL0WIN 


Obj ect Setup inlet flow. 

Options 

I0PT1 = 1 Compute inlet flow 

If (T =0) P =P and T -T 
o o r or 

I0PT1 = 2 Read inlet flow 

Xf (x o >10.) Normalize with T r and P r 


List of Symbols 

BINP (I , J,K) , Interpolated from BINPUT(IH, J,L) 

Theory I0PT1=1 

For this option the freestream flow is assumed to be isentropic with a 
constant freestream Mach number H g and wall boundary layers defined by power - 
law velocity profiles. Since the swirling flow must be in radial equilibrium, 
the normal momentum equations must be satisfied together with isentropic 
relationships of pressure ratio and temperature ratio. 


Neglecting curvature in the meridional plane 

in . JL « n A 

d n R an r 


where 


2 2 2 
= M - M s 


M 2 * 


a- 1 




Equation (1) can be integrated with 


( 1 ) 


( 2 ) 

(3) 


M(o) s M s /COS (a H ) 


(4) 


and 



(5) 


n 0 _ r 
nio ) 5 L 


i + 


Y~ 1 
2 



as initial conditions using a Runge-Kutta method. 


For a given displacement thickness and velocity profile power law, wall 
boundary layers can be added, assuming collateral boundary layers, such that 
a is unchanged. Then 


A* (l + n) A* 



(6) 

jjl . ( y y /n 2 
u© [ a ; 



(7) 

and 





( u y-1, e w -e„ 

\ UflD / ®C0 

[-Tfe] 

(8) 

Finally, the inlet mass flow and 
follows 

reference velocity are 

determined as 


. N e r'J.p u £2. 
“r - T Jo v P ‘ x 


(9) 


w = 9 Pr u r °r A (10) 


Theory I0PT1=2 

For this option, the input flow is calculated from experimental input 
data. The input variables selected are spanwise location, total pressure, 
static pressure, flow angle, and total temperature, since these are the pri- 
mary measured variables. A simple linear interpolation is used so that for 
any variable <f>. 


<£(y (7j)) = <£(y,) + [<£(y 2 )-<£(y,)][y(77)-y,] 


(ID 
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The flow variables are calculated from equations 1 through 7. 

If (T_ > 10) , it is assumed that pressure and temperature are given in 
(psf) and deg R, respectively, and the flow is normalized accordingly. If 6* 
is given, it is assumed that boundary layers should be added accordingly to 
the velocity profile power law above. Finally, the weight flow and reference 
velocity are determined from equations 9 and 10. A flow chart of this sub- 
routine is presented in figure 20. 
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Subroutine FN0RM 


Object Normalize input variables . 
Options 


None 

List of Symbols 

All variables in C0MM0N blocks 


Theory 

All input variables are normalized according to the List of Symbols. 
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Subroutine F0RCE/PR0P 


Obj ect Compute blade forces 

Options 


N0PPF = 1 
N0PPF = 0 

List of Symbols 
FRCI(N, 1) 
FRCI(N, 2) 
FRCI(N, 3) 
FRC(N, 1) 

FRC(N, 2) 

FRC(N, 3) 
F0RC(3, K) 
F0RC(4, K) 
FL0C(5 , N) 

FL0C ( 6 , N) 

K 

N 


Radial, axial, and swirl blade forces are defined 
in propeller portion of program. 

Radial, axial, and swirl blade forces are read 
from input. 


Blade force 
Blade force 
Blade force 
Blade force 
Blade force 
Blade force 
Blade force 
Blade force 
cosine of 6 


(radial direction) /span 
(phi direction) /span (dimensionless) 
(axial direction) /span (nondimension) 
(radial direction) /span (nondimension) 
(phi direction) /span (nondimension) 

(axial direction) /span (nondimens ional) 
(stream direction) /volume (nondimensionai; 
(swirl) /volume (nondimensional) 


sine of 0 

Index of n coordinate 

Index of points along propeller centerline 


Theory 

This subroutine is used to calculate the streamwise and swirl blade com- 
ponents for each streamline.lt is assumed that the normal component is small. 
The process has two steps. 

(1st) Calculation of radial, axial, and phi blade force components for each 
streamline by using linear interpolation of blade points. 

(2nd) Converting radial, axial, phi blade force components to streamwise and 
swirl blade forces. 
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Subroutine F0RCL 


Obj ect Compute local blade force. 
Options None 


List of Symbols 
GAP = G 

ZLE « Z T 


LE 


ZTE 

Theory 


= Z 


TE 


Strut gap (dimensionless) 

Location of leading edge (dimensionless) 
Location of trailing edge (dimensionless) 


The chordwise distribution of blade loading (force/volume) is defined by 






V-£t(C d T) 


pT* ds p 




( 1 ) 


( 2 ) 


I-I r = C P In (^-) - In (-£-) 


(3) 


In this subroutine, the chordwise loading is assumed uniform. 
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Function FTHIK (Z,IS,L0P) 


Obj ect Compute blade thickness distribution. 

Options 

IS=1 NASA 5 digit series distribution, 

IS=4 Input thickness distribution, 

IS=7 65-A series thickness distribution. 

L0P=1 FTHIK = t/t 
- max 

L0P=2 FTHIK = d (t/t max )/d(x/c) 

List of Symbols 

Z = x/c , Fractional chordwise distance, 

FTHIK = t/tmax , Ratio of thickness to maximum thickness of blade. 

Theory 

A thickness distribution of the form 

t/l: max = i* 4845 ft - .63 Z - 1.758Z 2 + 1.4215 Z 3 - .5075Z 4 (1) 

is used to represents a NASA 5 digit series distribution (IS=1). 


For a NASA 65A series distribution or for any arbitrary distribution 
(IS“4) table data is used to represent the distribution. 
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Subroutine GBLADE 


Obj ect 

Compute Blade Geometry 


Options 

None 



List of 

Symbols 



ALPHS 

a 

9 

Stagger angle to axis (deg) 

CHRD 

B 

9 

(Local) blade chord (dimensionless) 

GAP 

IRT 

ISHAPE 

G 

9 

9 

9 

Gap between blades (dimensionless) 
Index Counter 
Blade shape (Option) 

NBLADE 


9 

Number of blades in one row 

NRW 


9 

Number of blade rows 

NUM 


9 

Number of points along blade centerline 

PHIC 


9 

Blade camber 

RCLH 

^LH 

9 

Hub radius of blade centerline 

RCLT 

r clt 

9 

Tip radius of blade centerline 

SLD 

a 

9 

Solidity 

THICK 

t 

9 

Local blade thickness 

THICKN 

t n 

9 

Blade thickness 

THIKM 

t/B 

9 

Maximum thickness/chord 

ZBAR 

Z 

9 

Axial position with respect to blade 

ZCL 

Z CL 

9 

Blade axial centerline 

ZCLX 

Z CLX 

9 

Axial distance to blade centerline 

ZKK 


9 

Fractional distance along chordline 

ZLE 

Z LE 

9 

Blade leading edge 
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ZTE 


Z TE 

9 

Blade trailing edge 

CONST (1, 

L) 

r l 

9 

Radius (input data) 

CONST (2, 

L) 

“l 

9 

Stagger angle (input data) 

CONST (3, 

L) 

b l 

9 

Chord (input data) 

CONST (A, 

L) 

(t/B> L 

9 

Thickness to chord ratio (input data) 

CONST (5, 

L) 

♦l. 

9 

Camber angle (input data) 

CONST (6, 

L) 

Z CL 

9 

Axial distance (input data) 

Q(2, L) 


R 

9 

Radius 

Q(13, L) 


G 

9 

Gap 

Q(14, L) 


3G/3n 

9 

Normal derivative of blade surface 

Q(15, L) 
Theory 


9G/3S 

9 

Streamwise derivative of blade surface 


This module interpolates blade data with respect to the radius of blades 
in order to obtain local blade information. The following parameters are 
interpolated; a, B, t/B, <J> C , and Z^,^. The blade thickness is then calcu- 
lated in FTHIK and used to determine G, 9G/9n and 9G/9s. A flow chart is 
given in Fig. 21. 
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Obj ect Compute Duct shape 
Options 


Subroutine GDUCT 


I0PT3=1 Straight annular duct 
I0PT3=2 Read duct shape 
I0PT3=3 Straight-wall diffuser 
IOPT3=4 NACA curved -wall diffuser 

List of Symbols 

(As needed by user) 

Theory 

This subroutine is used to prescribe the duct shape fjj(Z) , r-p(Z) , 
wall bleed m H (Z), m T (Z), and wall temperature T H (Z), T T (Z), as required. 
Since these functions are input, the programmer may write a subroutine 
for this purpose or read the required information according to I0PT3. In 
addition, the subroutine computes the reference radius and normalizes the 
variables r and T. The variable m is normalized in Subroutine FL0WIN when 
U r is calculated. 

Input/Output 

The user may program any duct shape and wall boundary conditions as 
required. The output of this subroutine must be (R(I,K,J), 1=1,3; K=l,2, 
J=1,JL) and Zl. Note that all variables are normalized as shown in the 
sample subprogram described in the Subroutine GDUCT listing and that equally 
spaced spanwise stations are used. The flowchart (figure 22) should be 
followed in programming. 
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Subroutine GE0MCL 


Obj ect Calculate coordinates of lifting line. 
Options 


None 







List of Symbols 






AL 

= 

a l 


9 

Slope of lifting line 


AM 

= 

a m 


9 

Slope of coordinate line 


AN1, AN 2 

3S 

N l» 

n 2 

9 

N coordinate of grid 


AN1, ANI1 

= 

Nl, 

Nl-1 

9 

N coordinate of lifting line 

intersection with grid 

ANL 

- 

n l 


9 

N coordinate of input point of lifting line 

Rl, R2 

- 

R r 

R 2 

9 

R coordinate of grid line 


RI, RI1 

= 

R i« 

R I-1 

9 

R coordinate of lifting line 

intersection with grid 

RL 

= 

r l 


9 

R coordinate of lifting line 

input 

SI, S2 

- 

s i» 

S 2 

9 

S coordinate of grid 


SI, SI1 

as 

s i» 

S I-1 

9 

S coordinate of lifting line 

intersection with grid 

SSL 

= 

S L 


9 

S coordinate of lifting line 

input point 

Zl, Z2 

- 

Z l* 

Z 2 

9 

Z coordinate of grid line 


ZI, ZI1 

= 

Z I»' 

Z I-1 

9 

Z coordinate of lifting line 

intersection with grid 

C0NS(1,L) 

= 

r l 


9 

Lifting line radius (Local) 


C0NS(2,L) 

= 

“L 


9 

Stagger angle (Local) 


CONS (3 ,L) 

= 

b l 


9 

Chord (Local) 


C0NS(4,L) 


(t/B) L 

9 

Thickness /chord (local) 


CONS(6,L) 

= 

Z CL 


9 

Lifting Line Axial Distance (local) 
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C0NS(8,L) n L 

C0NS(9,L) S L 


Lifting line normal locations 
Lifting line streamwise location 


Theory 

The equation for the lifting line passing through the points (L, L-l) 
is approximated by a line given by 

Z = Z L _, + A l (R-R l _,) (u 


where the slope is given by 


( 2 ) 


The streamline coordinate system forms a box around the L-l mesh point 
(see figure 23) given by the points (1), (2), (3), (4). If it is assumed 
that point (5) is known; the object is to find point (6) by successively 
checking each side of the mesh to determine if the lifting line crosses. 
Let the index I = 1, 4 represent each side of the box. Then a straight 
line is passed through the pair of points 


I 

Points 


1 

(2), 

(1) 

n^ - S varies 

2 

(3), 

(2) 

S 2 - n varies 

3 

(A), 

(3) 

n£ - S varies 

4 

(1), 

(4) 

St - n varies 


The equation of the straight line approximating the lifting line is given by 

Z = Z, + A m -(R-R,) 


(3) 


A M = ( Z 2 - Z |)/( R 2 - R .) 
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Then the intersection of the lifting line with the coordinate line is 
given by 


R I " ( Z L-I “ A L* R L-I _Z I +A M R |)/ ( A M“ A l) 

Z I =[( 2 L-rA L - R l . ( ) • A M -(Z, - a m - R| ) . a l ]/(a m -a l ) 

provided 

a M“ A l^O 


(4) 

(5) 

( 6 ) 


From figure (23) R,Z will intersect inside the mesh if; 

R, < or r 2 Sr x <r ( 


(7) 


and 


Z,< z z <l z or Z 2 <Z z <Zi 


( 8 ) 


For 1=1 and 3, n is constant and S varies. Thus 

Sj = S, + (R i -R,)/(R 2 -R | ) (S 2 -S,) 

(9) 

n x = n, if I=i 
= n 2 if 1 = 3 

Likewise for 1=2 and 4, S is constant and n varies. Thus 

n i s n, +(R i -R,)/(R 2 -R,) (n 2 — n ,) 

S I = S 2 ifI = 2 < 10 > 

= S, if I =4 

Then Rj, Z r , Sj, nj is known for point (6) 

From figure 18, it can be seen that if 

R I-I- R L“ R I (11) 
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R l lies within the grid (1), (2), (3), (4), and L should be advanced, 
n and S vary along the lifting line. If (R^, lies within the grid 

then n^, can be determined by the relations 

S L = Sj + (S I ~Sj_|)/ (Rj~ Rj-|) (R(_“ ^l-0 

n L = n i-i + ( n i~ n i-,)/ ( R i - R I-|) (R l - R l _,) 

This procedure can then be repeated for L = 2, NUM 


Both 


( 12 ) 
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Subroutine INITQ 


Object 

Initialize data file parameters for Q array 
Option 


None 




Variables 




BLOCK(l) 

JGTEP 

» 

Block (record number) 

Q(1,K) 

R 


Radius (dimensionless) 

0(2, K) 

Z 

> 

Axial distance (dimensionless) 

Q(3 ,K) 

9R/9n 

» 

Derivative 

0(4, K) 

9R/9s 

y 

Derivative 

Q(5,K) 

(cos0)^ . 
axi 

y 

Axisymmetric flow angle 

Q(6,K) 

V 

y 

Metric coefficient (dimensionless) 

0(7, K) 

3V/9n = (K + AK ) 
s s 

y 

Curvature of streamline 

Q(8,K) 

3V/9s 

y 

Curvature of potential line 

Q(9,K) 

X 

y 

Distance along streamline (dimensionless) 

Q(10,K) 

Y 

y 

Duct height (dimensionless) 

Q(11,K) 

y/y t 

y 

Normalized duct height 

Q(12,K) 

A 

y 

Duct Area (dimensionless) 

Q(13,K) 

G 

y 

Gap (dimensionless) 

Q(14,K) 

3G/9n 

y 

Derivative 

Q(15,K) 

9G/9s 

y 

Derivative 

Q(16,K) 

3n/9n 

y 

Transformation of normal coordinate 
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Subroutine INITQ (Cont'd) 


Variables (Cont’d) 


Q(17,K) 


, Not used 

Q(18,K) 

n 

, Normal coordinate (dimensionless) 

Q(19,K) 

n 

, Transformed coordinate 


II 


QPARM(l) 

r 

r 

, Reference radius, (ft) 

QPARM(2) 


, Not used 

OPARM (3) 

JL 

, Number of streamwise steps 

QPARM(4) 

KL 

, Number of streamlines 

Theory 




This subroutine initializes the independent variable array BLOCK which is 
stored on a disc file and sets all parameters QPARM required by the calculation. 

* Note that Q(7,K) stores either K g or K g +AK S . 
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Subroutine INTFRE 


Obj ect 

Initialize freestream conditions 
Options 

None 

List of Symbols 

See C0MM0N BL0CKS 
Theory 

Data read from file NDRUM is used to set up the freestream conditions for 
subroutine P0IS. 
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Subroutine L0ADRR 


Obj ect Loader formated input 
Options 

N0PT8=O Continue reading input 
N0PT8=1 Stop reading input 

List of Symbols 

See EQUIVALENCE ARRAYS in subroutine listings. 

Theory 

This method of reading input permits the changing of one or more input 
variables. The remaining input variables remain the same as the previous 
case. 
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Subroutine MINVRT(A,B,N) 


Obj ect Invert NxN matrix. 

Options 
None 

List of Symbols 
A = A 

B « B 

N 
M 

Theory 

The A matrix is inverted using the Gauss-Jordan elimination procedure. 
First the augmented ^(N, M)is formed including the identity matrix, 

A =(A I) (1) 


, Augmented A matrix 
, Augmented B metric (A~l) matrix 
, Number of equations (rows) 

, Number of columns 


Then the following revision formula is used 


= a iw-°w o v /o i l 


Vj-I ^ a g /Q U 


1 < J £ M 


( 2 ) 

( 3 ) 


Note that the I matrix has one less column than the X matrix. Then the 

'U 

substitution is made 


a 


ij 


'ij 


I < T< N 
I < J < M -I 


And repeated until 


the B matrix is an NxN or the A' 


matrix. 
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Subroutine MYTIME 


Object Dummy time trap routine 


Subroutine 0UTPUT 


Object 

Print title page 
Option 
None 

List of Symbols 
None 
Theory 

This subroutine prints the title page which records all modifications, dates, 
and references to changes incorporated into the ADD code. 
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Subroutine PERFNA 


Object Compute viscous nacelle drag 
Options 
None 


List of 

Symbols 





AREAM 

= 

A 

max 

J 

Maximum cross-sectional area 


CDFR 

= 

C DF 

> 

Friction drag coefficient 


CDPR 

= 

C DP 

9 

Pressure drag coefficient 


DFR 

= 

d f 

9 

Friction drag (lb) 


QREF 

= 

Qr 

9 

Reference dynamic pressure (psf) 


RMAX 

= 

R 

max 

9 

Maximum radius of nacelle 


AVE(9,1) 

= 

p i 

9 

Initial mass flow average density 


AVE(4 ,1) 

= 

u i 

9 

Initial mass flow average velocity 


P 

= 

p 

9 

Static pressure 


P1,P2 

= 

* ^2 

9 

Static pressure (P/P r ) at prescribed 

stations 

R 

= 

R 

9 

Local radius 


R1,R2 

= 

R 1 >r 2 

9 

Nacelle radius (r^/r r ) of prescribed 

stations 

S 


S 

9 

Streamwise coordinates 


S 1» S 2 

= 

S 1 ,S 2 

9 

Streamwise coordinate at prescribed stations 

Theory 







In terms of dimensionless variables used in the analysis, the friction and 
pressure drag are given by: 

D f - 2tt p x u r 2 J*R £ns ds (1) 

S. 
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where 

A MAX = w r r 2 *MAX 

O r = \/Zp r u}ppf 


( 5 ) 

( 6 ) 
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Subroutine PERFN2 


Obj ect Compute inviscid nacelle drag 
Options 
None 

List of Symbols 


DPR 

P 

P1,P2 
Rl, R2 
S 

S 1» S 2 

Theory 


P 1’ P 2 


Rl , R£ 

S 

t S2 


, Pressure drag 
, Static Pressure 

, Static pressure (P/P r ) at prescribed stations 
, Nacelle radius (ri/r r ) at prescribed stations. 
, Streamwise coordinate. 

, Streamwise coordinate at prescribed station. 


In terms of the dimensionless variables used in the analysis, the 
pressure drag is given by; 


D p = 2 7T 



ds + n, 


Rf 

~2 


-n$j 


( 1 ) 
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Subroutine P0IS (RESM, ITER) 


Object 

Solve Poisson equation 
Option 

IDBGP =0 No debug printout 

= 1 Printout residuals 

= 2 Print solution 

List of Symbols 


P(K,J) 

= 


, Stream function (dimensionless) 

F(K,J) 

= 


f Coefficient (1/P G) (dimensionless) 

PSI(K) 

= 

Or 

, Iterative guess for J 

ITER 

= 

V 

» Iteration counter 

RLX 



* Relaxation factor 

RESMAX 

= 

e MAX 

, Maximum residual accepted 

RESDM 

= 

e M 

, Maximum residual/ J station 

RESM 

= 

e 

» Maximum residual/sweep 

Theory 





The solution algorithm is described in Reference 1. 
References 


1. Anderson, 0. L. and D. E. Edwards: Extension to an Analysis of Turbulent 

Swirling Compressible Flow in Axisymmetric Ducts, NASA Contract NAS3-21853, 
1981, UTRC Report R81-914720. 
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Subroutine P0ISCF 


Object 

(1) Set of coefficients of = 0 

(2) Set boundary conditions on 
(2) Set initial guess for ^ 

Options 



IDBG17 

= 0 

Compressible Flow 



= 1 

Incompressible Flow 


IDBGT 

= 0 

No debug test 

case 



= 1 

Debug test case 


IDBGP 

= 0 

No debug test 

printout 



= 1 

Debug printout 

Variables 




Q(I,K) 



, Coordinate functions 

FIV(I 

,L,K) 



, Dependent variables for inviscid flow 

P(K) 

= 

K 


, Stream function at station J 

A(K) 

= 

GP/V/(dn/dn) 

.Coefficient for dip/dr) 

F(K) 

= 

i/g/p 

.Coefficient of ^ 

G(K) 

= 

v/g/P 

.Coefficient for velocity calculation 

R(K) 

= 

P 


.Density ratio (P/P r ) 

T(K) 


0 


.Temperature ratio (T/T r ) 

V(K) 

= 

V 


.Metric coefficient dimensionless 

GAP 

= 

G 


.Gap (g/r r ) 

VMET 

= 

V 


.Metric coefficient (dimensionless) 

RH0 

= 

P 


.Density ratio (p/P r ) 

TEH 

= 

0 


temperature ratio (T/T r ) 

USO 

= 

U 

so 


.Upstream constant velocity (u 0 /u r ) 


184 



Subroutine P0ISCF (Cont'd) 


US INF 

U 

s°° 

>Free 

stream 

axial velocity ( u soo / u r) 

UPINF 

V 

,Free 

stream 

tangential velocity u^/uj. 

PSIKL 


,Free 

stream 

stream function (dimensionless) 

TEMINF 


,Free 

stream 

static temperature ratio (T/T r ) 

RHOINF 

K 

,Free 

stream 

density ratio (P/b r ) 

AMINF 


,Free 

stream 

Mach number 

PTINF 

"co 

,Free 

stream 

total pressure ratio (P /P r ) 

BLK 


,See COMMON/ SPCGD/ 

Theory 






This subroutine does the following steps 

(1) Reads coordinate Q file and solution FIV file 

(2) Interpolates the solution to the (n,S) grid 

(3) Calculates the coefficient F = 1/PG 

(4) Calculates coefficients BLK for streamline curvature calculation 

(5) Sets boundary condition on ¥ 

(6) Calculates initial guess for ^ 

(7) Stores F, BLK, P on disk files 

The initial guess is given by the inviscid solution obtained from CALINV. 
boundary conditions are given by: 


^(0,s) = 0 

(1) 

V|/( 1 ,s) = 

(2) 

1 dr, 

0 

(3) 


(4) 


o 


The 
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Subroutine P0IS0N 


Object 

Calculate axi symmetric streamline curvature 
Options 


I0PT7 = 0 

* 1 
IDBG15 = 0 

> 0 


No curvature corrections 
Curvature correction 
Use input KL streamlines 
Use IDBG15 streamlines 


List of Symbols 


IRH0 

.Density 

iteration counter 

ITERL 

.Maximum 

number of iterations 

KH0LD 

.No. ADD 

code streamlines 

KL 

.No. SCURVA streamlines 

RESMAX 

, Maximum 

residual for convergence 


Theory 

This subroutine is a calling subroutine for subroutines INTFRE, P0ISCF, 
P0IS , and SCURVA. 
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Subroutine QINTER 


O bjec t 

Interpolate curvature from P0IS mesh to S0LVI mesh 
O ptions 

None 

L ist of Symbols 

Q(J,K) Coordinate functions 

Theory 


After the curvature and flow angle has been calculated from the potential 
flow solution,, this subroutine interpolates to obtain values at the numerical 
grid points which will be used in the S0LVI calculation. 
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Subroutine READPF(J.JJ) 


Object 

Read P and F files in NIST word blocks 


Option 

None 

List of Symbols 

J J 


JJ 

N 

F 

P 

NST 

NBK 

NIST 

NFDRM 

NL 

NBIST 

NMOVE 


JJ 

N 

F 

25 


23 


.Record number 

.Record number in block N 

.Block number 

2 

.Coefficients of V i|> *= 0 
.Stream function (dimensionless) 
.Number of records per block 
.Number of words to move pointer 
.Number of words per block 
.Coefficient file number F array 
.Last block number 
.Number of words for two records 
.Number of words to move pointer 


Theory 

The entire F and P arrays cannot be kept in core at the same time so that the 
I/O is arranged to keep fixed blocks in core (Fig. 1). Let (J.K) be a point on 
the computational mesh and (JJ.KK) the corresponding point in core. Let each 
record be the Jth line with the number of words in the record given by 

K * l.IST 

If there are NST records per block, then these are NIST words per block. 


NIST = NST x 1ST 
188 
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Subroutine READPF(J,JJ) (Cont'd) 

Theory (Cont f d) 

The solution algorithm requires overlapping blocks as shown on Fig. 1. 
Hence we have the block number 

N - (J-2) / (NST-2) +_1 
and the JJ point in core is given by 

JJ + J-(N-l) x (NST-2) 

This subroutine is coded so that a new block is ready only when N = NL. 
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Subroutine READPG(J,JJ) 


Object 

Read variable for curvature calculation 
Options 
None 

List of Symbols 


J 

= J 

, Record number 

JJ 

= JJ 

, Record number in block N 

N 


, Block number 

G 


.Streamline coordinate data (see COMMON /SPCFD/) 

P 


.Stream function 

NST 

= 25 

.Number of records per block 

NBK 


.Number of words to move pointer 

NIST 


.Number of words per P block 

NGDRM 

* 25 

.Unit number of G array 

NPDRM 

= 24 

.Unit number for P array 

NL 


.Last P block number 

NBIST 


,No. words for 2 P records 

NGIST 


,No. words for 2 G records 

NGL 


.Last G record number 

Theory 



This subroutine 
reads only the Jth G 

reads the P file according to Subroutine READPF but 
record which is kept in core. 
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Function R03RTS(C,ETA,L0P) 


Object 

Compute distorted mesh using Roberts' transformation 
Options 

L0P = 0 Wall - wall boundary 

= 1 Wall-free stream boundary 
=-l Free stream-wall boundary 

List of Symbols 

C = C .Distortion parameter 

ETA = n .Input variable (uniform mesh) 

L0P = .Option 


Output Variable 

R0BRTS = n .Output variable 


Theory 

The Roberts’ transformation for a distorted mesh on both sides is given by 

(c+ 1 / 2 ) exp [2 in ( ' %” ]/ f ~ ) (i 7 '-i/ 2 )]-(c-l /2 

' + e *p[2 ln (f^|-) (V- '/2>] 

where 


0 <rj < 1.0 


0 < n < 1.0 


For the different options we have 

17 ^ } L0P = 0 

n = n J 

t7’ = t?/2 ^ 

' , f L0 P = 1 

n = 2n J 

V = (l + 7))/2 


n = 2n - 1 


} L0P = -1 


C'l> 
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(2) 

(3) 

(A) 

( 3 ) 



Subroutine R0UND 


Obj ect Round corners on straight wall ducts. 

List of Symbols 

XM , Axial location 

IWALL , Indicates Hub or Tip Wall 

DX , Stepsize in x direction 

Theory 

In subroutine GDUCT, if the option I0PT3 * 8 is used then a straight 
wall is constructed between initial data points. Thus in order to remove 
discontinuity between wall segments this subroutine is used to round or 
smooth out the discrete representations of the wall in order for it to 
appear smooth. 

Subroutine SLETE(KSSLE,KSSTE) 

Object Find blade control surfaces. 

Options 
None. 

List of Symbols 
KSSLE.KSSTE 
SLE,STE 
ZLEH,ZLET 

ZTEH.ZTET 

Theory 

The intersection of the leading and trailing edge of the blade with the 
hub and tip casing are obtained from Subroutine GBLADE. Then the coordinates 
of the hub and tip boundaries are searched until the proper value of streamwise 


, Leading edge and trailing edge index 

, Leading and trailing edge coordinates (dimensionless) 

, Axial distance hub leading and trailing edge 
(dimensionless) 

, Axial distance tip leading and trailing edge 
(dimensionless) 
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coordinates for the leading edge and trailing edge of the blade are found. 
The coordinate index KSSLE is located just upstream of the blade and the 
coordinate KSSTE is located just downstream of the blade. 
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Subroutine SCURVA (IDBU.KH0LD) 


Object 

Calculate curvature from potential flow solution 
Options 

IDBU = 0 Update density 

> 0 Print SCURVA solution, update curvature 

Variables 


KH0LD 


, No. ADD code streamlines 

KL 


,No. SCURVA streamlines 

Q(J,K) 


, Coordinate functions 

POO 

*£ 

, Stream function a station J (dimensionless) 

GOO 

G 

Coefficient for velocity calculation 

R(K) 

P 

, Density Ratio (P/P r ) (dimensionless) 

US 

U s 

,Streamwise velocity (dimensionless) 

UN 

U n 

, Normal velocity (dimensionless) 

U 

U 

, Total velocity (dimensionless) 

COSTH 

CD 

S-X 

(Nl 

(0 

o 

o 

2 

, (Cosine) of flow angle 0 

CURV 

9V 

3n 

, Curvature of streamline (dimensionless) 


Theory 

Once the potential flow solution has been obtained from subroutine P0IS, 
this subroutine will calculate the flow angle and streamline curvature according 
to Ref. (1). 

References 


1. Anderson, O.L. and D. E. Edwards, Extension to an Analysis of Turbulent 
Swirling Compressible Flow in Asixymmetric Ducts, NASA Contract No. 
NAS3-21851, 1981, UTRC Report R81-914720. 
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List of S> 

(HTlbols 


JX 


, Number of input points 

JXB 


, Number of output points 

JXK 


, Number of spline knots 

X(J) 

II 

c* 

, Input points abscissa 

Y (J) 

" Y J 

, Input points ordinate 

XB(J) 

II 

C*"' 

, Output points abscissa 

YB(J) 

1 

II 

, Output points ordinate 

YPP (J) 

11 

, Second derivative of Y(X) 

CK(I, J) 

= C(I,J) 


YK(I) 

II 

7 ** 

, Spline coefficients 

A(I) ,B (I) 

“ V B i 

, Constants of Integration 

Theory 




This subroutine computes the second derivative of the input vector Y(X). 
With the use of the standard math package ISML routines, (reference 9) it 
then fits a least square spline to the second derivative with JXK movable 
knots. The spline equations are then integrated analytically to obtain the 
output solution vector Y(X) at JXB points. Subroutine SM00TH uses ISML 
routines ICSFKU, ICSFKV, UERTST. 
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Subroutine S0LVI 



List of Sy 

mbols 

AA(I,J) 

= 

a IJ 

AB(I,J) 

s 

b IJ 

AC(I, J) 

= 

C IJ 

AD(I , J) 

s 

d IJ 

ADI(I,J) 

= 

(«u\ 

AE(I,J,K) 

- 

e IJ 

AQ(I) 

- 

<1 

AZ(I,K) 

* 

Z K 

C T. 

CFPH 

= 

C f <J>H 

CFPT 

= 

C f <|>T 

CFST 

as 

C fSH 

CFST 

= 

C f ST 

DAYE 


1/4 (p U x 2 ) 

EENTP 

= 

E(I) 

ENREF 

as 

f l 

EPRES 

- 

E(tt) 

ER0TH 

as 

e(P0) 

ER0US 

= 

e(PU s ) 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 




Element of A matrix 

as 

Element of B matrix 

Element of C matrix 

Element of D matrix 
s-1 

Element of D matrix 
Element of E matrix 
Element of Q matrix 

Sir 

Element of Z matrix 

- - 2 

Stress coefficient hub (2 W< p i u i » 

Stress coefficient tip (-2e n ^/(P^U^^)) 

— _ O 

Stress coefficient hub (2e ns /(P )) 

Stress coefficient tip (-2e ns /(PjU-^)) 

Mean inlet dynamic pressure (dimensionless) 
Truncation error (it) 

Mean inlet entropy 
Truncation error (ir) 

Truncation error (P0) 

Truncation error (PU S ) 
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EUPUP 

«=> 

, 2 n 
E(U0 ) 

t 

, 2 n 

Truncation error (U ( / ) 

EUSUS 

= 

E (u s 2 ) 

9 

2 

Truncation error (U s ) 

PIREF 

= 

iri 

9 

Mean inlet reference pressure (dimensionless) 

PIO 

= 

n 

o 

9 

Total pressure (dimensionless) 

PRCEF 

= 

C P 

9 

Local pressure coefficient (dimensionless) 

PSIH1, 

PSIH2 


'J’H.’i'H 1 

9 

Wall stream function (hub) (dimensionless) 

PSIT1, 

PSIT2 

= 

^,^-i 

9 

Wall stream function (tip) (dimensionless) 

QAVE 

= 

PU s (tt 0 -5) 

9 

Inlet energy flux (dimensionless) 

QWALH 

= 

Ph u h * 3 

9 

Energy flux (hub) (dimensionless) 

QWALT 


p t u t * 3 

9 

Energy flux (tip) (dimensionless) 

QPLUS 

= 

Q/(PU*3) 

9 

Universal heat flux (dimensionless) 

SIG 


e 

9 

Stress (dimensionless) 

SIGWH 

= 

e WH 

9 

Wall stress (hub) (dimensionless) 

S1GWT 

= 

Sjt 

9 

Wall stress (tip) (dimensionless) 

ST AH 


S tH 

9 

Stanton number (hub) Qfl/ [PU g (6 o -0) ] 

STAT 

\ 

= 

S tT 

9 

Stanton number (tip) Q^/[PU s (6 o -0) ] 

THETAO 

= 

0 o 

9 

Total temperature (dimensionless) 

TPLUS 

= 

T + 

9 

Universal stress (dimensionless) 

TWH,TWT 

S 

®H» 0 T 

9 

Wall temperature (hub, tip) (dimensionless) 

U 

= 

U 

9 

Magnitude of velocity (dimensionless) 

UPLUS 

= 

U+ 

9 

Universal velocity (dimensionless) 

USH.UST 

= 

uh,ut 

9 

Friction velocity (hub, tip) (dimensionless) 
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XMACH = M 

YPLUS * Y + 

Z = Z 

ZZ « Z 


, Mach number (dimensionless) 

, Universal distance (dimensionless) 

, Axial distance (dimensionless) 

, Axial distance to next slot (dimensionless) 


Theory 

The equations of motion are solved using the method derived in 
reference 6. A flow diagram is shown in Fig. 24. 
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Subroutine STRESI 


Obj ect Compute initial stress distribution 
Options 

L0P-1 Store inlet flow stress 
L0P#1 Do not store inlet flow stress 
List of Symbols 

See C0MM0N block variables in subroutine listings. 

Theory 

The initial stress and heat flux distribution is computed from 



Subroutine STRT 


Object Find inlet flow location 
Option 

None 

List of Symbols 

ZINLET , Axial inlet flow location 

SINLET , Streamwise inlet flow location 

INLET , Counter 

Theory 

Once the axial location of the inlet flow is determined the streamwise 

location may be found since the axial location at each point in the coordinate 

system is paired with its streamwise coordinate. 

Subroutine TPRINT 

Object Calls CPU time. 
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Subroutine TURB 


Obj ect Compute turbulent viscosity 
Options 

N0PT=O Initial turbulence model 
N0PT=1 Subsequent turbulence model 


List of Symbols 


AMUE 

AMUER(K) 

AMUM 

AMUW 

AMUWK 

AMUO 

DELO 

DU 

DUDZ 

E 

EM 

EMH.EMT 

ENP 

ENS 

EW 

PHI 

RH0M 


- Pj/P r 


(P T / P_) 

X oo X 


= vjv r 


% /V T 




(p T /P ) 

loo » 


= A 


AU 

dU/dZ 

E 

Eco 

E coH» E «T 


n<f> 


= E. 


ns 


w 


= p. 


, Turbulent viscosity (dimensionless) 

, Freestream turbulent viscosity (dimensionless) 

, Freestream molecular viscosity (dimensionless) 

, Wall value of molecular viscosity (dimensionless) 

, Molecular viscosity (dimensionless) 

, Maximum freestream viscosity (dimensionless) 

, Displacement thickness (dimensionless) 

, Velocity finite-difference (dimensionless) 

, Velocity derivative (dimensionless) 

, Rate of strain (dimensionless) 

, Strain fretestream (dimensionless) 

, Strain hub, tip, edge of inner layer (dimen- 
sionless) 

, Swirl rate of strain (dimensionless) 

, Streamwise rate of strain (dimensionless) 

, Wall rate of strain (dimensionless) 

, Turbulence model function (dimensionless) 

, Density freestream (dimensionless) 
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SIGWH 


, Wall stress (hub) (dimensionless) 


SIGWK 

SIGWT 

TPLUS1 , 
TPLUS2 

UK(K) 

UM 

USTARH, 
U START 


— T+ ■p"*’ 
A 1» A 2 


■ Vi 


, Wall stress (inner layer) (dimensionless) 

, Wall stress (tip) (dimensionless) 

, Universal wall stress (dimensionless) 

, Magnitude of velocity (dimensionless) 

, Freestream velocity (dimensionless) 

, Friction velocity (hub , tip) (dimensionless) 


Y = Y , Distance across duct (dimensionless) 


YK 

II 

$ 

, Distance across duct (dimensionless) 

YMH.YMT 

- y h> y t 

, Distance to inner layer (hub, tip) (dimensionless) 

YPLUS1 , 



YPLUS2 

= Y Y 
1’ 2 

, Universal distance (dimensionless) 

Theory 




The turbulence model is described in reference (1) and the resulting 
equations are described below. Let the eddy viscosity be described by a con- 
tinuous function 
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(outer layer) 


(4) 


4> 


j[W* 

Em 


'V 

where y is the distance from the wall 


y = |y-y w 


(5) 


A matching point for the inner layer and outer layer exists for each wall 
denoted and Y T and with a corresponding strain and E^. Then for the 
outer layer 


E„ * Eh 


* 

yT yH 


( 6 ) 


The turbulent viscosity and thermal conductivity is given by 


Jfl JL + 

M r Mr 


Me 

Mr 


RE 


rL 


Mr 


J_ ±1 

J RT Mr 


(7) 


( 8 ) 


The turbulent flow properties can be calculated at station J-l because the 
flow field is known. For station J, it is noted that the turbulent viscosity 
is a strong function of stress, thus from equation (1), 


(-£)’■ N »* z 


(9) 


where 


2 1 yF„; + 


( 10 ) 


Hence, 


/Ml' \ J . 


r * 

/ Me' 

\ 

J-l 

/ 

f ~ J - J-!\ 


'Mr/ 

KF7) * 

|_az 

iMr > 

) 

1 

, z - z j 

(11) 

1 Mr y . 

(' JtL.) 

J-l 

+ 

Laz ' 

r Me^ 1 


(12) 

P RE Mr > ' 

' P RE 1 


' Mr /. 

\ J 
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Finally, 
is not known. 


it is noted that at the initial station t (Y ) in equation 3 
therefore, t + = 1 is assumed. 



Function UBLAS 


Ob j ect Calculate velocity ratio according to the Blasius solution for 

the initial flow. 

Options 


None 

Theory 

The velocity ratio is determined by interpolations of data block con- 
taining Blasius solution to the flow past an ax i symmetric shape (reference 10). 
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Function UC0LES (Argument List) 

j ec t Find friction velocity for the initial profile from Coles law 


Options 


List of Symbols 


ITERL 


= dE/dU* 


Von Karman constant (dimensionless) 

2 

Kinematic viscosity (ft /sec) 

Boundary layer thickness (ft) 
Displacement thickness (ft) 

Slope of error function (dimensionless) 
Error function (dimensionless) 
Convergence criteria (dimensionless) 
Iterate 

Maximum number of iterations 
Freestream velocity (ft/sec) 

Guess for friction velocity (ft/sec) 
Initial guess for U* (ft/sec) 


Theory 


The friction velocity is obtained from Coles law using Newton's method 


(U*) M 


tSH) «> 


V wn* 


2 S* _.\ U» 

S I (u* 2 )^ 
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( U *)"* 1 = ( u *)' 1 . i */( 

' Vdu* 


Convergence occurs when 


M< 
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Subroutine WAKC0R 


Obj ect Compute nacelle wake corrections. 
Options 

None 


List of Symbols 
AM1.AM2 = M k _ ls M k 

9 

Inviscid flow Mach number 

DRL 

= 

AR l 

9 

Radial wake correction 

DZL 

= 

az l 

9 

Axial wake correction 

ETA 


Tl 

1 

Normal coordinate 

G(1,J,L) 

= 


9 

Wake distance (radians) 

G(2 , J ,L) 

= 

Alp 

9 

Tangential wake correction 

G(3, J,L) 

= 

d\p/ds 

9 

Partial derivative 

G(4, J,L) 

= 

dAip/ds 

9 

Partial derivative 

G(5, J,L) 


AR 

9 

Radial wake correction 

G(6,J,L) 

= 

AZ 

9 

Axial wake correction 

RL 

= 

r l 

9 

Radial location of lifting line streamline 

SL 

= 

S L 

9 

S coordinate of lifting line streamline 

S2,S1 

= 


9 

S coordinate of mesh 

T0T1,T0T2 

= 

(T o /T) K -1’ (T ° /T) K 

9 

Total to static temperature ratio 

Tl, T2 

= 

T K-1 * t k 

9 

Static temperature 

Ul, U2 

= 

U K-1> U K 

9 

Inviscid flow velocity 

US1, US2 

- 

Us K-l’ Us K 

9 

Inviscid flow streamwise velocity 

UP1, UP 2 

= 

U <t»K-l’ U< f> K 

9 

Inviscid flow tangential velocity 

VL 

= 

V L 

9 

tn 

Metric scale coefficient of L streamline 
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ZL 


, Axial distance L t ^ 1 streamline 


J 

K 


, Index of S coordinate 


, Index of N coordinate 
th 

, Index of L streamline 


L 

Theory 

To obtain the corrections to the wake geometry due to the nacelle's presence 
in the flow field, this subroutine integrates the following equations, 
(formulted in reference 1). 


s j 


^ J u 


ds 


S L ’ S 


( 1 ) 


/ <J® nc 

RU S “v* 

S L 

using the trapezoid rule along the streamlines passing through the point 

(R , Z ) or (S , n ). The remaining wake corrections are given by 
L L L 


( 2 ) 


AR l = R l (Sj)- r l 


^ Z L= Z L (Sj)-Z l (S l )- 


(3) 

(A) 


Labeled Common Blocks Used in The Nacelle Portion 


Included herein is an alphabetical list of the labeled common blocks used 
in the nacelle portion of the analysis and a description of each variable used 
in them. 


COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


AC0NS (Blade Data - C0NSTI (1, 
Dimensionalized ) 

C0NSTI (2, 
C0NSTI (3, 
C0NSTI (A, 
C0NSTI (5, 
C0NSTI (6, 


L) 

r cl 

Radius of propeller lifting 
line 

L) 

= a 

s 

Stagger angle 

L) 

= B 

Chord 

L) 

= t/ B 

Thickness 

L) 


Not used 

L) 

- Z CL 

Axial location of propeller 
lifting line 
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Subroutine WBLEED 


Object 

Calculate perforated wall bleed 
Options 


I0PT18 

= 

0 

No wall bleed 


= 

1 

Tip wall bleed 


= 

2 

Hub wall bleed 


= 

3 

Tip/hub wall bleed 

PCHEK 

> 

1.0 

Flow enters tunnel 


< 

1.0 

Flow leaves tunnel 

Input Variables 


AHAS 

= 

V A s 

Ratio of hole area to surface area 

CDISH 

= 

c 

Discharge coefficient 

pip 


P TP 

Plenum total pressure 

HP 


T x 

Plenum total temperature 

Internal Variables 


AMTU 

= 

M 

TU 

Tunnel Mach number 

GAMMA 

= 

Y 

Ratio of specific heats 

GASR 

= 

R 

Gas constant 

PS 

= 

P 

Static pressure (psfa) 

PT 

= 

P T 

Total pressure (psfa) 

PSTU 

= 

P TU 

Tunnel static pressure (psfa) 

PTTU 

= 

P 

TTU 

Tunnel total pressure (psfa) 

TSTU 

= 

T 

TU 

Tunnel static temperature (deg R) 

TTTU 

= 

T 

TTU 

Tunnel total temperature (deg R) 

RH0R 

= 

r 

2 

Reference density (slug/ft ) 

USR 

= 

u 

r 

Reference velocity (ft/sec) 

PRESR 

= 

P 

r 

Reference pressure (psfa) 
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Internal Variables (Cont'd) 


TEMPR = T 

r 

SGN * 1 

Output Variables 

RH ( 9 , J ) » (PU n ) H 

RT(9,J) = (P.U ) T 
n 1 


Reference temperature (deg R) 

Sign convention 

2 

Mass bleed hub wall (slugs/ft /sec) 

2 

Mass bleed tip wall (slugs/ft /sec) 


Theory 


If one treats a single hole in a perforated wall as an orifice, then the mass 
flow can be derived in terms of the plenum stagnation conditions and the local static 
pressure inside the tunnel Holman (Ref. 1). Then an expression for the mass flow 
added to the tunnel flow is given by 


</>U n ) W 


A h y p t 

A s yr*T T 



(i) 


where and T T are the plenum conditions, P is the local tunnel static pressure, 
A^/A s is the ratio of the hole area to surface area, and C the effective discharge 
coefficient which is a property of the perforated wall. If the tunnel static pres- 
sure is greater than the plenum total pressure, the mass flow bleed is out of the 
tunnel. Under these conditions, and Trp are taken from the wind tunnel condi- 
tions, and P is the plenum pressure which is assumed known. 

The mass flow bleed is related to the stream function by 

d* _ G (f u n)w ... 

ds " V pr u r (2) 


Equations (1) and (2) provide the boundary condition for a perforated wall relating 
two dependent variables ip and P in terms of the characteristics of the perforated 
wall and the plenum conditions. 

The program checks the options according to the table below. 
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O.D. Wall 


I.D. Wall 


PCHEK = P ,,/P^,, > 

1.0 set 


= 

p mT , 




TP STU 


T 


TP 






P 

= 

P 

STU 

SGN = 

-1.0 

+1.0 



T 

_ 

T 






T 


TP 




PCHEK = P„„/P„ mT , < 

1.0 set 


— 

P _ ___ 




TP STU 


T 


STU 






P 

= 

P TP 

SGN = 

1.0 

1 

• 

O 



T 


T 






T 


STU 





Reference 

1. Holman, J. P.: Experimental Methods for Engineers. McGraw-Hill Book Co., 

New York. 1966. 
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r 


\ 


{ 

Subroutine WRITPF(JJ) 


Object 

Store updated potential flow solution. 
Options 


None 


List of Symbols 
JJ 
MST 
NIST 

NPDRK = 24 
P 


,JJth station in core 
,No. records per block 
,No. words per block 
,Unit number 
, Stream function 


Theory 

The stream function array P(JJ,KK) is arranged in core as described in 
subroutine READPF, Fig. 1. When an iterative sweep of one block is 
complete, the new updated solution is written on a disk file. This occurs 
when JJ = NST-1. 
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Function XH(J) 


Ob,j ect Calculate wall length on ID wall 
Variables 


J 

XH 

Theory 


Wall point no. 

AX r Wall length 


4x h ■ <(r hj 


- Vi )2 + “ 2) 


1/2 


Function XT(J) 


Object Calculate wall length on OD wall 
Variables 


J 

XT 

Theory 


AX„ 


Wall point no. 
Wall length 


2 2 
AX_ * ( (R T - R T ,) + (Z mT - Z rnT ,) ) 

T TJ TJ-1 TJ TJ-1 


1/2 


HU 



COMMON BLOCK NAME 
(OBJECT) 


AC0NS (Blade Data - 
Dimensional) 


AC0NX (Blade Data - 
Nondimens ional) 


ADPS (Coordinate 
for Slot Calcula- 
tions) 

AINV (Store Inviscid 
Flow Variables) 


VARIABLE NAME 
C0NSTI (7, L) 
0MEGZI = fl 

NUM 

L = 1, 2 


DESCRIPTION OF 
VARIABLES 


Not used 

Propeller rotational velocity 
Number of blade rows 


C0NST (1, L) 

C0NST (2, L) 
C0NST (3, L) 
C0NST (4, L) 
C0NST (5, L) 
C0NST (6, L) 

C0NST (7, L) 
C0NST (8, L) 
C0NST (9, L) 


Radius of propeller lifting 
line (dimensionless) 

Stagger angle (dimensionless) 

Chord (dimensionless) 

Thickness 

Not used 

Axial location of lifting line 
(dimensionless) 

Not used 

Normal location of lifting line 

Streamwise location of lifting 
line 


L = 1,2 
0MEGZ 
NUMX 
DPSI(K) 

K = 1, KL 
CINP (1, K) 
CINP (2, K) 
CINP (3, K) 


Rotational velocity (dimensionless) 

Number of propellers 

Radial coordinate (dimensionless) 

Total pressure inviscid flow 
Static pressure inviscid flow 
Swirl angle 
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COMMON BLOCK NAME 
(OBJECT) 

VARIABLE NAME 

DESCRIPTION OF 
VARIABLES 


CINP (4, K) 
K - 1, KL 

Total temperature inviscid 
flow 

AKK (Complex Coor- 

Ml 

Flags for DEBUG printout 

dinate) 


M2 

See Subroutine FCPLX 

AMATRX (Matrix Inver- 

AD(I, J) = d, T 

Element of D matrix 

sion) 


ADI(I,J) = djj 

Element of D matrix 

APL0T (Store Variables 

W(l, J) - z R 

Axial distance (hub) (dimen- 

for Plotting) 

sionless) 


W(2, J) = Z T 

Axial distance (tip) (dimension- 
less) 


W(3, J) = C pH 

Pressure coefficient (hub) 
(dimensionless) 


W(4, J) = C pT 

Pressure coefficient (tip) 
(dimensionless) 


W(5 ’ J) = C FSH 

Friction coefficient (hub) 
(dimensionless) 


W(6, J) = C F$T 

Friction coefficient (tip) 
(dimensionless) 

BCPLX (Complex Vari- 

A(l, I) = A ± 

Source strength (dimensionless) 

ables) 


A(2,I) = b ± 

Location of pole (dimensionless) 


A(3,I) - a. 

Wall angle change (deg) 


A(4,I) = r i 

Radius in z plane (dimension- 
less) 


A(5,I) * <t> ± 

Angle in z plane (radians) 


A(6,I) * x ± 

Relative x distance in z plane 
(dimensionless) 
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COMMON BLOCK NAME 
(OBJECT) 


BCPLX (Complex Vari- 
bles) 


DESCRIPTION OF 

VARIABLE NAME VARIABLES 


A(7,I) = 

Y i 

Relative Y distance in z plane 
(dimensionless) 

B(1,I,K) 

= ASX s 

Change in coordinate x (dimen- 
sionless) 

B(2,J,K) 

" ASY s 

Change in coordinate y (dimen- 
sionless) 

B(3,I,K) 

" AS ^s 

Change in coordinate £ (dimen- 
sionless) 

B(4,I,K) 

= ASn s 

Change in coordinate n (dimen- 
sionless) 

X(1,K) = 

s 

Streamwise coordinate (dimen- 
sionless) 

X(2,K) = 

n 

Normal coordinate (dimension- 
less) 

X(3,K) = 

X 

X coordinate in z plane (dimen- 
sionless) 

X(4,K) = 

Y 

Y coordinate in z plane (dimen- 
sionless) 

X(5,K) = 

5 

£ - coordinate in w plane (dimen- 
sionless) 

X(6,K) = 

n 

_ coordinate in w plane (dimen- 
sionless) 

X(7,K) = 


Streamwise derivative of £ 
(dimensionless) 

X(8,K) = 

n s 

Streamwise derivative of n 
(dimensionless) 
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COMMON BLOCK NAME 
(OBJECT) 


BLEED (wall bleed) 


BTHIK 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


X(9,K) ■ 

■ X(S+AS ) 

X(10,K) 

= Y(S+AS) 

X(11,K) 

= £ (S+DS ) 

X(12,K) 

= n(S-tDS) 

X(13 ,K) 

= Y 

X(14,K) 

^ss 

X(15,K) 

^sn 

X(16,K) 

" V n 

X(17,K) 

= v s 


Coordinates at station S+DS 

Coordinates at station S+DS 

Coordinates at station S+DS 

Coordinates at station S+DS 

Metric scale coefficients 
(dimensionless) 

Second derivative of £ 
(dimensionless) 

Cross derivative of £ (dimen- 
sionless) 

Normal derivative of V (dimen- 
sionless) 

Streamwise derivative of V 
(dimensionless) 


AHAS 

CDISH 

FT? 

TTP 


■ V A s 

- c 



Ratio of hole area to surface area 
Discharge coefficient 
Plenum total pressure 
Plenum total temperature 


KB LADE 


9 

Number of points 

XK(I) 

' X I 

» 

Fractional chordwise distance 

YK(I) 

“ Y I 

9 

Thickness/Chord 

I » 

1, KB LADE 
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COMMON BLOCK NAME 
(OBJECT) 


CCPLX (Parameters 
for Complex Trans- 
form) 


DESCRIPTION OF 


VARIABLE NAME 

VARIABLES 

NPTS = N 

P 

Number of singularities in 
complex transformation 

NS0URC = N s 

Number of sources in z plane 

0RDER1 = 0-,^ 

Absolute magnitude of largest 
term 

0KDER2 = 0 2 

Absolute magnitude of largest 
term 

0RDER3 = 0 3 

Absolute magnitude of largest 
term 

SLO = S LQ 

Length of streamwise coordinate 

XDS = S 

Step size for complex integra- 
tion 

XDN = S 

Step size for complex integra- 
tion 
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I 


COMMON BLOCK NAME 
(OBJECT) 

CINF (Parameters 
for Poisson Equation) 


DESCRIPTION OF 

VARIABLE NAME VARIABLES 


AMINF 


M 

GO 

Freestream Mach number 


AN(K) 


n k 

Transverse coordinate 

DEDN(K) 


(dn/dn) k 

Transverse coordinate stretching 

D2EDN(K) 


(d 2 n/dn 2 ) k 

Transverse coordinate stretching 

PINF 


ji 

Freestream static pressure 


Ji oo 

(dimensionless; 

PSIKL 


*« 

Freestream stream tunction 
(dimensionless) 

PTINF 


n« 

Freestream total pressure 
(dimensionless) 

RH0INF 


P« 

Freestream density 




(dimensionless) 

R0TINF 


P 

Freestream total density 


* o°° 

(dimensionless) 

TEMINF 



Freestream static temperature 
(dimensionless) 

TT1NF 



Freestream total temperature 
(dimensionless) 

UINF 


TT 

Freestream velocity 


u 

CD 

(dimensionless) 

UP INF 


v 

Freestream tangential velocity 
(dimensionless) 

US INF 


u e „ 

Freestream streamwise velocity 


Soo 

(dimensionless) 

UUO 

WO 


U o 

Reference velocity 
(dimensionless) 

* 

V o 

Reference metric coefficient 
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COMMON BLOCK NAME 
(OBJECT) 

VARIABLE NAME 

DESCRIPTION OF 
VARIABLES 

C0NST (Flow Constants) 

A CHI = x 


Clauser constant (0.016) 


AKAPPA = 

K 

von Karman constant (0.41) 


APLUS = 

+ 

A 

van Driest constant (26.0) 


CPR = 

C 

P r 

Specific heat at constant 
pressure (5997.0 ft^/sec^/deg R) 


CVR = 

C 

vr 

Specific heat at constant 
volume (3283.0 ft^/sec^/deg R) 


EP * 

0 e 

2-7182818 


GAMMA « 

a 

Ratio of specific heats (1.4) 


GASR = 

R - 

2 2 

Gas constant (1714.0 ft /sec / 
deg R) 


GRAVR 


Gravitational constants (32.2 
f t/sec^) 


PI = 

77 

3.1415926 


PRESR = 

P 

r 

Reference static pressure (psfa) 


PRL = 

Pr 

L 

Prandtl number laminar 0.70 


PRT = 

Pr 

T 

Prandtl number turbulent 0.72 


RH0R = 

P 

r 

Reference density (slugs/ft ) 


SNDR « 

C 

r 

Reference speed of sound (1116.0 
ft/sec) 


TEMPR = 

Tr 

Reference temperature (deg Raukine) 


TI 


(0.1745329 radians/deg) 


VISCR - 


Reference molecular viscosity 
(0.370 x 10" 6 ) 
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COMMON BLOCK NAME 
(OBJECT) 

C£)RE (Coordinate 
Functions) 


DESCRIPTION OF 


VARIABLE NAME 

VARIABLES 

Q(1,K) = R 

Radius (dimensionless) 

Q(2,K) = Z 

Axial distance (dimensionless 

Q(3,K) = 9R/9n 

Normal derivative of radius 
(dimensionless) 

Q(4,K) = 9R/9S 

Streamwise derivative of 
radius (dimensionless) 

Q(5,K) = 9 2 R/9n/9S 

Second derivative of radius 
(dimensionless) 

Q(6,K) = V 

Metric scale coefficient 
(dimensionless) 

Q ( 7 , K) « 9V/9 n 

Curvature of potential line 
(dimensionless) 

Q(8,K) = 9V/9S 

Curvature of streamline 
(dimensionless) 

Q(9,K) = 9 2 V/9n/9S 

Second derivative of metric 
scale coefficient (dimensionless) 

Q(10,K) = Y 

Physical distance across duct 
(dimensionless) 

Q(11,K) = Y/Y t 

Fractional distance across 
duct (dimensionless) 

Q(12,K) = A 

Area between adjacent stream- 
lines (dimensionless) 

Q(13,K) = G 

Gap between blade surfaces 
(dimensionless) 

Q(14,K) = 9G/9n 

Normal derivative of blade 
surface (dimensionless) 
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COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


QU5,K) = 

9G/9S 

Streamwise derivative of blade 
surface (dimensionless) 

Q(16,K) = 

9n/9n 

Transform of normal coordinate 
(dimensionless) 

Q (17 ,K) = 

9 2 n/9n 2 

Second derivative (dimensionless) 

Q(18,K) = 

n 

Normal coordinate (dimensionless) 

Q(19,K) = 

n 

Transformed normal coordinate 
(dimensionless) 

K = 1,KL 


Number of streamlines (dimensionless) 
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COMMON BLOCK NAME 
(OBJECT) 


CORE 2* (Wall Value of 
Coordinate Functions) 


* Note: Actually three 


DERIV (Force Functions) 


DESCRIPTION OF 


VARIABLE NAME 

VARIABLES 

R0(1,K) = R 

Radius (dimensionless) 

R0(2,K) = dR/dZ 

Derivative of radius (dimen- 
sionless) 

R0(3,K) = d 2 R/dZ 2 

Second derivative of radius 
(dimensionless) 

R0(4,K) = Z 

Axial distance (dimensionless) 

R0(5,K) = V 

Metric scale coefficient 
(dimensionless) 

R0 ( 6 , K) = dY/dS 

Derivative of metric scale 
coefficient (dimensionless) 

R0(7,K) = Y t 

Distance across duct (dimen- 
sionless) 

R0 ( 8 , K) = S 

Streamwise coordinate (dimen- 
sionless) 

R0 ( 9 , K) = m 

Mass flow bleed (dimensionless) 

R0(1O,K) =© w 

Wall temperature (dimensionless) 

R0S (I) = R0(I,K) 

Dummy Storage Vector 


arrays defined where 0 takes on the value H, M, T 


0 = H, Hub wall 
= M, Mean line 
= T, Tip wall 


DF (1 ,K) = 

H /XV 
L s J 

J 

K-h 

Streamwise blade force/volume 
(dimensionless) 

DF (2 ,K) = 

. H * /xv . 

J 

K-^2 

Tangential blade force/volume 
(dimensionless) 

DF(3,K) = 

[v ct 
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J 

K-l 

Total pressure loss/volume 
(dimensionless) 



COMMON BLOCK NAME 
(OBJECT) 


DRED1 (Store Flow 
Variables) 


DRED2 (Store Flow 
Variables) 


DUC0UT (Wall 
Coordinates) 


VARIABLE NAME 

DESCRIPTION OF 
VARIABLES 

df(4 ’ k) ■ [ x ];U 

Coordinate distortion 
(dimensionless) 

DF(5,K) = [H/Xvf- 1 

L s Jk - } 5 

Streamwise blade f orce/volume 
(dimensionless) 

DF(6,K) - [V»g 

Tangential blade force/volume 
(dimensionless) 

df(7,k) - [yxy]^ 

Total pressure loss/volume 
(dimensionless) 

DF(8 - K) ■ L X JS 

Coordinate distortion (dimen- 
sionless) 

BL0CK (I) 

Grid structure variables 

11 = 19 * KL + 35 


I = 1, 11 


BL0CK1 (I) 

Grid structure variables 

11 = 19 * KL + 35 


I = 1, 11 


R(l, 1, J) = R t (Z j ) 

Radius of hub (dimensionless) 

R(2,1,J) = R h (Z j ) 

Radius of tip (dimensionless) 

R(1,2,J) = S T (Z j ) 

Mass flow of tip bleed (dimension- 
less) 

R(2 ,2 , J) = ^(Zj) 

Mass flow of hub bleed (dimen- 
sionless) 

R(1,3,J) = ^(Zj) 

Wall temperature of tip 
(dimensionless) 

R(2 ,3, J) = 0 T (Z J ) 

Wall temperature of hub 
(dimensionless) 

R(1,4,J) = Zj 

Axial location of hub 

R(2,4,J) = Zj 

Axial location of tip 
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COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


DUCTIN (Input Functions) 


DUCTIX (Input Functions) 


EBLAD (Blade Row 
Parameter) 


BINI(L) 

I 

J 

K 

L 

DUCTI(I) I 


BINPUT (I , J,K) 

1,5 
1,2 
1,KLL 

5* (K-1)+I+(J-1)*5*KLL 
1,15 


RD1I(L) 

= R(1,1,L)*RADR 

RD2I(L) 

= R(1,2,L)*RADR 

ZD1I(L) 

= R(1,4,L)*RADR 

ZD2I (L) 

= R(2,4,L)*RADR 

L 

- 1, JL 

AINI (L) 

= AINDUCT (I,J,K) 


L = 1, JL 

NR0W 

LR0W 
LBLD 
NBLADE 


, see C0MM0N/SPI0/ 


, Arbitrary duct geometry 
parameters 

, Tip wall coordinates (ft) 

, Hub wall coordinates (ft) 

, Tip wall coordinates (ft) 

, Hub wall coordinates (ft) 

, See C0MM0N/SPI0X/ 


, Number of Streamwise 
stations 

Maximum number of propeller 
Rows 

Indicates propeller row 
Not used 

Number of blades in propeller 
row 
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COMMON BLOCK NAME 
(OBJECT) 

VARIABLE NAME 

DESCRIPTION OF 
VARIABLES 

FC0R (Truncation 

^ =E U0 

Error in swirl velocity 

Error Estimates) 

(dimensionless) 


EAS =E US 

Error in streamwise velocity 
(dimensionless) 


EEN = Ej 

Error in entropy (dimensionless) 


EPO - Ep 0 

Error in total pressure 
(dimensionless) 


ESI - E^ 

Error in stream function 
(dimensionless) 


ETO = E to 

Error in total temperature 
(dimensionless) 
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COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


FIVC (Inviscid Flow FIV(1,L,K) 

m 

* 


Stream function (dimensionless) 

Variables I/O) 

FIV(2,L,K) 

m 

u 

s 


Streamvlse velocity 
(dimensionless) 

FIV(3,L,K) 

m 

U <> 


Tangential velocity 
(dimensionless) 

FIV(4,L,K) 

m 

n 


Static pressure (dimensionless) 

FIV(5,L,K) 

m 

I 


Entropy (dimensionless) 

FIV(6,L,K) 

m 

e 


Static temperature (dlmensionlet 

FIV(7,L.K) 

m 

p 


Density (dimensionless) 

FIV(8,L,K) 

m 

M 


Mach number 

FIV(9,L,K) 

m 




FIV(10,L,K) 

m 




W 

@ 

j-: 

station 


L-2 


j 

station 


L-3 

Q 

j+l station 


K-1,KL 



streamlines 


FIPARM(l) 

m 

'r 


Reference density (slugs/ft^) 

FIPARM(2) 

m 

T 

r 


Reference temperature (deg R) 

FIPARM(3) 

m 

P r 


Reference pressure (psfa) 

FIPARM(4) 

m 

g 


Gravitational constant (ft/sec^ 

FIPARM(5) 

m 

u 

r 


Reference viscosity 
(slugs/ft/sec) 

FI P ARM (6) 

m 

S 


Specific heat constant pressure 

FIPARM(7) 

m 

C v 


(ft 2 /sec 2 /deg) 

Specific heat constant volume 

FIPARM(8) 

m 

R 


(f t 2 /eec 2 /deg) 

Gas constant (ft 2 /sec 2 /deg) 
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COMMON BLOCK 
(OBJECT) 


NAME 


DESCRIPTION OF 

VARIABLE 

NAME 

VARIABLES 

FIPARM(9) 

Pr I 

Prandtl number (turbulent) 

FIPARM(IO) 



FIPARM(ll) 

u r 

Reference velocity (ft/sec) 

FIPARM(12) 

N0PT7 

Number of stations stored 
in inviscid solver 

FIPARM(13) 



FIPARM(14) 

I0PT15 

First station 

FIPARM(15) 

I0PT16 

Last station 
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COMMON BLOCK NAME 
(OBJECT) 


FLAGS (Flags) 

FL0WI (Flow Input 
Functions) 


F0RS (Blade Force 
Variables) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


N0PT0 0 = 1,28 

FG(1,K) = a 
FG(2,K) = n Q 

FG(3,K) =© 0 

FG(4,K) = M 
FG(5,K) = P Q 

FG(6,K) = U 

K = 1, KL 

F0RC(1,K) = H s 

F0RC(2,K) = 
F0RC(3,K) = 

F0RC(4,K) = ^ 

F0RC(5,K) = $ B 

FORC ( 6 , K) = 


Flags to regulate calculation 
flow (see subroutines) 

Inlet swirl angle (deg) 

Inlet stagnation pressure 
(dimensionless) 

Inlet stagnation temperature 
(dimensionless) 

Inlet Mach number (dimensionless) 

Inlet stagnation density (dimen- 
sionless) 

Inlet magnitude of velocity 
(dimensionless) 

Number of streamlines 

Streamwise force/area (dimen- 
sionless) 

Swirl force/area (dimensionless) 

Streamwise force/span (dimen- 
sionless) 

Swirl force/span (dimension- 
less) 

Blade dissipation/area (dimen- 
sionless) 

Blade dissipation/span (dimen- 
sionless) 
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COMMON BLOCK NAME 

(OBJECT) VARIABLE NAME 


FORS2 (Blade Force Variables) 


FF(1,I,K) 

A 

* M 

FF(2,I,K) 

A 

« n 

FF(3,I,K 

A 

= e 

FF(4,I,K) 

A 

- e o 

FF(5,I,K) 

1) 

3> 

o 

FF(6,I,K) 

A 

= P 

FF(7 ,1 >K) 

A 

- u s 

FF(8,I,K) 

A 

FF(9,I,K) 

A 

FF(10,I,K) 

R 

FF(11,I,K) 

A 

= a 

FF(12,I,K) 

A 

= ft 

FF(13,I,K) 

A 

« I 

FF(14,I,K) 

A 

= U 


DESCRIPTION OF 
VARIABLES 


Inviscid Mach number (dimen- 
sionless) 

Inviscid static pressure 
(dimensionless) 

Inviscid static temperature 
(dimensionless) 

Inviscid total temperature 
(dimensionless) 

Inviscid total pressure 
(dimensionless) 

Inviscid density (dimen- 
sionless) 

Inviscid streamwise 
velocity (dimensionless) 

Absolute swirl velocity 
(dimensionless) 

Relative swirl velocity 
(dimensionless) 

Blade velocity (dimen- 
sionless) 

Absolute angle to axis 
(deg) 

Relative angle to axis 
(deg) 

Inviscid flow entropy 
(dimensionless) 

Magnitude of relative 
inviscid flow velocity 
(dimensionless) 
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(Continued) 


COMMON BLOCK 
(OBJECT) 

NAME 

VARIABLE NAME 

DESCRIPTION OF 
VARIABLES 

FORS2 (Blade 

Force Variables) 

FF(15,I,K) 

II 

Loss coefficient (dimen- 
sionless) 



FF(15,2,K) 

II 

> 

Blade entropy rise 
(dimensionless) 



FF(16,I,K) 

= ^ 

Stream function (dimen- 
sionless) 



FF(17,1,K) 

= C L 

Lift coefficient (dimen- 
sionless) 



FF(17,2,K) 

= C D 

Drag coefficient (dimen- 
sionless) 



I 

= 1 

Upstream of blade row 



I 

= 2 

Downstream of blade row 



K 

= 1, KL 

Number of streamlines 
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COMMON BLOCK NAME 
(OBJECT) 


FUNC (Dependent Flow 
Variables) 


INTINP (Integer Input 
Variables) 


VARIABLE NAME 
K * 1,KL 

F(1,I,K) « ifr 
F(2 ,I,K) = Ug 

F(3,I,K) = 

F(4,I,K) = n 
F(5,I,K) = I 
F(6,I,K) = 0 

F ( 7 , 1 , K) = P 
F(8,I,K) = E ns 

F(9,I,K) = 
F(10,I,K) = Q 
I =1 

I =2 

I =3 

K = 1, KL 

I0PT0 = 1, 17 

IDBG0 = 1, 23 


DESCRIPTION OF 
VARIABLES 

Number of streamlines (dimen- 
sionless) 

Stream function (dimensionless) 

Streamwise velocity (dimension- 
less) 

Swirl velocity (dimensionless) 

Static pressure (dimensionless) 

Entropy (dimensionless) 

Static temperature (dimension- 
less) 

Density (dimensionless) 

Streamwise stress (dimension- 
less) 

Swirl stress (dimensionless) 
Heat flux (dimensionless) 

Inlet conditions 
S=S 

S=S+dS 

Number of streamlines 
Input /Output options 
Debug printout options 


ISHAPE 

JL 

KDS 

KL 


Blade shape option 

Number of streamwise stations 

Number of steps per streamwise 
station 

Number of streamlines 
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COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


KLL Number of input data streamlines 

NB Number of blades 


INTPAR (Parameter 
Variables) 

These variables are treated as PARAMETER statements in UNIVAC programs 
and as integer variables in IBM or CDC programs. 


Parameters for Solution Arrays 

PARAMETER 1ST = 100 
PARAMETER IS = 100 
PARAMETER NEQ = 10 
PARAMETER NBCH - 5 
PARAMETER NBCT = 5 
PARAMETER NBH1 = NBCH+1 = 6 
PARAMETER NEQD = 2*NEQ =20 

Parameters for Slots 

PARAMETER ISL0T = 15 
PARAMETER ISL0T2 = 2*ISL0T = 30 
PARAMETER IS1 = 6*ISL0T+2 = 92 
PARAMETER IS2 = 2*IS+IS1 = 292 


Parameters for Coordinate Arrays 

PARAMETER ISM = 19*IST = 1900 
PARAMETER ISL = ISM+35 = 1935 
PARAMETER IS3 = ISM+2 = 1902 
PARAMETER IS4 = IS 3+10 =1912 
PARAMETER IS5 = IS4+10 = 1922 
PARAMETER IS 6 = IS5+10 = 1932 


Parameters for Matrix Inversion 
PARAMETER KKLP =30 

PARAMETER LNGT0 = NEQ*NEG*KKLP+NEQ*KKLP =3300 
PARAMETER LNGT1 = NEQ*NEQ* KKLP+1 = 3001 
PARAMETER LNGT2 = 2*LNGT0 = 6600 
PARAMETER LNGT3 = 3*IST*NEQ = 3000 
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COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


j Parameters for Force Variables 

| PARAMETER IFFS = 34*IST = 3400 

i 

Parameters for Smooth 


PARAMETER ISSD = 20 
PARAMETER ISSD1 = ISSD+1 =21 

PARAMERER ISSD2 = ISSD1+ISSD = 41 
PARAMETER ISSD3 = ISSD2+ISSD = 61 


REALIN (Real 
Input Variables) 


ACI 

* X 

Clauser constant 

AKI 

= K 

von Karman constant 

ALP1 

" a x 

Inlet swirl angle hub (deg to 
z axis) 

AMS1 

= m i 

Inlet Mach number (dimensionless) 

ANH 

= n H 

Power law of hub boundary layer 

ANT 

= tliji 

Power law of tip boundary layer 

API 

+ c 

II 

van Driest constant 

CPRI 

= c 

P r 

SDecific heat constant pressure 
(f t 2 /sec 2 /deg R) 

CVRI 

= c v 

r 

Specific heat constant volume 
(ft 2 /sec 2 /deg R) 

DDS 

= 

Mesh distortion parameter 

DSHI 


Displacement thickness hub (ft) 

DSTI 

* 

6 t 

Displacement thickness tip (ft) 

PRESO 

II 

O 

Inlet stagnation pressure 

PRLI 

= P RL 

Prandtl number laminar 

PRTI 

= P RT 

Prandtl number turbulent 
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COMMON BLOCK NAME 
(OBJECT) 


REALIX (Real Input 
Variables) 


DESCRIPTION OF 

VARIABLE NAME VARIABLES 


TEMPO 

II 

H 

O 

Inlet stagnation temperature 
(deg R) 

VISCRI 

= y r 

Molecular viscosity reference 
(Ib/sec ft 3 ) 

ALPHSI 

= a CH 

Blade stagger angle hub 
(deg to z axis) 

ALP SMI 

a CM 

Blade stagger angle mid (deg 
to z axis) 

ALP ST I 

= a CT 

Blade stagger angle tip (deg 
to z axis) 

C0RDHI 

" b h 

Blade chord hub (ft) 

C0RDMI 

= B 

M 

Blade chord midpoint (ft) 

C0RDTI 

H 

pq 

ii 

Blade chord tip (ft) 

PHICHI 

= *CH 

Blade camber hub (deg) 

PHICMI 

^CM 

Blade camber midpoint (deg) 

PHICTI 

^CT 

Blade camber tip (deg) 

RCLHI 

r CLH 

Hub radius of blade centerline 
(ft) 

RCLMI 

r CLM 

Midpoint radius of blade center 
line (ft) 

RCLTI 

= r CLT 

Tip radius of blade centerline 
(ft) 

THIKHI 

= t H /B H 

Blade thickness to chord ratio 


hub (dimensionless) 
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COMMON BLOCK NAME 



DESCRIPTION OF 

(OBJECT) 

VARIABLE NAME 

VARIABLES 


THIKMI 

" t M /B M 

Blade thickness to chord 
ratio midpoint (dimensionless) 


THIKTI 

” t ip / Bij 

Blade thickness to chord 
ratio tip (dimensionless) 


ZCLHI 

= Z CLH 

Blade Axial location of center- 
line (Hub) 


ZCLMI 

Z CLM 

Blade axial location of center- 
line (Midpoint) 


ZCLTI 

= Z CLT 

Blade axial location of center- 
line (Tip) 


ZCLI 

Z CL 

Blade centerline location (ft) 

SPCFD (Variables 

F(K,J) 

- i/(p v)£ 

Coefficients of Poisson equation 

in Poisson Equation) 


T 



P(K,J) 


Stream function 

SPCGD (Variables 
for Streamline 

F(K) 

- 1/(P V) K 

Coefficient of Poisson equation 

Curvature) 

G(K) 

- V/(G >K 

Coefficient for velocity 


P(K) 

K 

Density ratio (p/p r ) 


T(K) 

K 

Temperature ratio (T/T^) 


V(K) 


Metric coefficient 

SPI0 (Flow Variables) 

AVE(1,J) 

- Z 

Average axial location (dimen- 
sionless) 


AVE(2,J) 

= AR 

Area ratio (dimensionless) 


AVE(3 , J) 

- 

Mass flow (dimensionless) 


AVE(A,J) 

- 5 s 

Average streamwise velocity 
(dimensionless) 


AVE(5,J) 

= {L 

0 

Average swirl velocity 
(dimensionless) 
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COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 
AVE(6, J) = n 
AVE ( 7 , J ) = I 

AVE(8,J) = 0 

AVE(9, J) = p 
AVE(10,J) = M 

AVE(11,J) = ff 0 

AVE(12 , J) = e Q 

AVE(13 , J) = Cp 

AVE (14, J) = C P;L 

AVE(15, J) = Z 

AVE(16, J) = B 

AVE(17 , J) = M 

J - 1, JL 

B INPUT (1,J,K) = Y 

BINPUT (2 , J,K) = II 0 
BINPUT(3, J,K) = n 
BINPUT (4, J,K) * a 


DESCRIPTION OF 
VARIABLES 


Average entropy (dimensionless) 

Average entropy (dimensionless) 

Average static temperature 
(dimensionless) 

Average density (dimensionless) 

Average Mach number (dimen- 
sionless) 

Average total pressure (dimen- 
sionless) 

Average total temperature 
(dimensionless) 

Average pressure coefficient 
(dimensionless) 

Average total pressure loss 
(dimensionless) 

Diffuser effectiveness (dimen- 
sionless) 

Blockage (dimensionless) 

Area average Mach number 
(dimensionless) 

Spanwise location (dimension- 
less) 

2 

Total pressure (lb/ft abs) 

2 

Static pressure (lb/ft abs) 
Swirl angle (deg to axis) 
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VARIABLE NAME 
BINPUT(5,J,K) = 0 Q 
J = 1 
J = 2 

K - 1, 

AINPUT(1,J,K) = Y 
AINPUT(2,J,K) = n Q 
AINPUT(3,J,K) = n 
AINPUT (4 , J ,K) = a 
AINPUT (5,J,K) = 0 

o 

J = 1 
J = 2 

K - 1, 

STRMES (Poisson BP0IS = B 

Stretching Parameter) 

BP0ISI = Bj 


DESCRIPTION OF 
VARIABLES 

Total temperature (deg R) 

Inlet flow 
Exit flow 

KLL Number of spanwise stations 

Spanwise location (dimensionless) 

2 

Total pressure (lb/ft abs) 

2 

Static pressure (lb/ft abs) 

Swirl angle (deg to axis) 

Total temperature (deg R) 

Upstream of blade row (dimen- 
sionless) 

Downstream of blade row (dimen- 
sionless) 

KLL Number of spanwise stations 

Stretching parameter 
Input stretching parameter 


COMMON BLOCK NAME 
(OBJECT) 


SPI0X (Flow Variables) 
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COMMON BLOCK NAME 
(OBJECT) 


STRES (Functions of 
Orthogonal Coordinates) 


SVARB (Parameters 
and Variables) 


DESCRIPTION OF 


VARIABLE 

NAME 

VARIABLES 

G(1,K) 

= 

Hk- 

-*5 

G(2,K) 

= 

MiU 

G(3,K) 

= 

[g(XV)J 

J 

K-h 

G(4,K) 

= 

ri 3v* 

|_XV 5n. 

u 

G(5,K) 

= 

r*l 3R' 

Lxr an. 

J 

K-h 

G(6,K) 

= 

r 1 9R- 

Lxr as. 

u 

G(7,K) 

•1 

_a_ 

.XG 3n 

MU 

G(8,K) 

= 

"v_ _d_ 

_XG ^n 

(G.i avl-J 
v xy 'anjK-h 

G (9 ,K) 

= . 

'JLjl 

XG Sn 

G 1 3Rl J 

V”XR anjK-Js 

G (9+1 ,K) 

= 

G(I,K) 

@ J-l, 1=1,9 


K - 1 , KL 
ALPHS = a 

ALPLUM = 

AMACHR = M 

r 

AMACH1 = M x 
AMPLUM 

AMPLUS = M + 

APRES1 = P 

AREA1 = A l 

AREAR = A 

r 
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Stagger angle to axis (deg) 

Not used in this version 

Reference Mach number (dimen- 
sionless) 

Average inlet Mach number 
(dimensionless) 

Not used in this version 

Mass flow bleed parameter 
(dimensionless) 

Average inlet static pressure 
(dimensionless) 

o 

Inlet area (ft ) 

2 

Reference area (ft ) 


COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


ATEMP1 

* Tl 

Average inlet static pressure 
(dimensionless) 

CH0RD 

= B 

Local strut chord (dimension- 
less) 

DETA 

= An 

Step size in normal coordinate 
(dimensionless) 

DMPULS 

n 

Step size in asymptotic constant 
table (m + ) (dimensionless) 

DPLUSM 

= AP + 

Step size in asymptotic constant 
table (p + ) (dimensionless) 

DS 

= AS 

Streamwise step size between 
stations (dimensionless) 

DSH 

* 

= a h 

Displacement thickness hub 
(dimensionless) 

DSS 

= dS 

Streamwise step size (dimension- 
less) 

DST 

* 

Displacement thickness tip 
(dimensionless) 

DZ 

= AZ 

Axial step size (dimensionless) 

DYNP1 

- % 

Average inlet dynamic pressure 
(dimensionless) 

GAP 

= G 

Gap between blades (dimensionless) 

GMRl 

= (Y _1 )M r 2 


GMR2 

= YM r 2 


JLAST 


Number of streamwise stations 

JLPTS 

= 

Number of input wall points 

JSEP 

s 

Number of stations to last 
calculated point 

KMH 

= 

Hub matching point 
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COMMON BLOCK NAME 
(OBJECT) 


VARIABLE NAME 


DESCRIPTION OF 
VARIABLES 


KMT 

s 

Tip matching point 

KMO 

= 

Not used in this version 

KSEP 


Not used in this version 

LM 

= 

Size of table of constants for 
inner layer (dimensionless) 

LMM 

= 

Midpoint of table of constants 
for inner layer (dimensionless) 

PHIC 

= 't’c 

Blade camber (deg) 

PLUSM 

= 

Not used in this version 

PPLUS 

- p + 

Pressure gradient parameter 
(dimensionless) 

RADR 

= r 

r 

Reference radius of blade 
centerline (dimensionless) 

REY 

= n r 

Reynolds number P^R^U /U 
r r r r 

RH0R 

= P r 

3 

Reference density (slugs/ft ) 

SL 

= S L 

Length of duct in streamline 
coordinates (dimensionless) 

S0LD 

= a 

Solidity (dimensionless) 

THICK 

= t 

Local blade thickness (dimen- 
sionless) 

USTARH 

* 

= U H 

Friction velocity hub (dimen- 
sionless) 

U START 

* 

" U T 

Friction velocity tip (dimen- 
sionless) 

USR 

= U r 

Reference radius (dimensionless) 

WFL0 

= 

Weight flow (lb/sec) 

YPLUSM 

~ 

Matching point for table 
asymptotic constants (dimension- 
less) 


242 



COMMON BLOCK NAME 
(OBJECT) 

VARIABLE NAME 

DESCRIPTION OP 
VARIABLES 


ZLE 

= 

Axial location of blade 
trailing edge (demensionless) 


ZTE 

= 

Axial location of blade leading 
edge (dimensionless) 


Z1 

= 

Duct axial length (dimension- 
less) 

SVARBX (Blade 
Variables) 

ALPSH 

* °H 

Hub stagger angle to axis 
(deg) 


ALPSM 

“ a M 

Midpoint stagger angle to axis 
(deg) 


ALPST 

= 

Tip stagger angle to axis (deg) 


CH0RDH 

- b h 

Blade chord hub (dimensionless) 


CH0RDM 

" ®M 

Blade chord midpoint (dimension- 
less) 


CH0RDT 

= b t 

Blade chord tip (dimensionless) 


PHICH 

' 4>ch 

Blade camber - hub (deg) 


PH I CM 

= ^CM 

Blade camber - midpoint (deg) 


PHICT 

“ ^CT 

Blade camber - tip (deg) 


RCLM 

= r clm 

Midpoint radius of blade 
center (dimensionless) 


RCLT 

" r clt 

Tip radius of blade centerline 
(dimensionless) 


RCLH 

r clh 

Midpoint radius of blade center 
(dimensionless) 


THICKH 

= C H 

Blade thickness hub (dimension- 
less) 


THICKM 

- c m 

Blade thickness midpoint (dimen- 
sionless) 


THICKT 

fc T 

Blade thickness tip (dimension- 
less) 
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COMMON BLOCK NAME 
(OBJECT) 


TITLIN (Input Title) 

TURBS (Turbulent 
Viscosity and Conduc- 
tivity) 


DESCRIPTION OF 

VARIABLE NAME VARIABLES 


ZCLH 

Z CLH 

Hub axial location of blade 
center line 

ZCLM 

Z CLM 

Midpoint axial location of blade 
centerline 

ZCLT 

Z CLT 

Tip axial location of blade 
centerline 

ZCL 

= Z CL 

Axial distance to blade center- 
line (dimensionless) 

TITLE(12) 

a , J-l 
3 /M 

Any alphanumeric characters 

DHF ( 1 , K) 

= ilfe) ^ 

Derivative of viscosity 

DHF(2,K) 

j-i 

= (y T /y r ) K-i 5 

Turbulent viscosity (dimensionless) 

DPF(1,K) 

= 


DPF (2 ,K) 

1 *1—1 

= (— — ) j Turbulent conductivity (dimen- 
RE "^sionless) 


K = 1, KL 
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List of Flags N0PT0 


Flag Name 

N0PT1=O 

N0PT1=1 

N0PT2=1 

N0PT2=2 

N0PT3=1 

N0PT3=2 

N0PT4=O 

N0PT4=1 

*N0PT5=O 

N0PTS>O 

N0PT6 

N0PT7=1 

N0PT7=2 

N0PT8=O 

N0PT8=1 

N0PT9=O 

N0PT9=1 

N0PT1O= 

N0PT11-O 

N0TP11-1 

*N0PT5 is given 
In addition, a 


Purpose 

Setup initial flow 
Calculate flow at station J 
Adiabatic wall 
Wall heat transfer 

Compute and store coordinate functions 

Compute local coordinate functions 

Duct with centerbody 

Duct with no centerbody 

Continue calculating 

Stop calculating 

Not used 

Flow separated from tip wall 

Flow separated from hub wall 

Read inviscid flow variables 

Read viscous flow variables 

Full complex function calculation 

Shorten complex function calculation 

Counts number of cases calculated 

Computer graphics I/O 

Batch I/O 

a value to locate error in program 
diagnostic message is printed. 
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N0PT13 

N0PT14=O 

N0PT14=1 

N0TP15=O 

N0PT15=1 

N0PT16 

N0PT17=O 

N0PT17=1 

N0PT18=O 

N0PT18=1 

N0PT19=INFL1 

N0PT2O=O 

N0PT2O-1 

N0PT21=O 

N0PT21=1 

N0PT22=O 

N0PT22=1 

N0PT22=2 


FCPLX flag 
UN I VAC 
IBM 

Integrate along S to obtain streamline 

Integrate along n to obtain streamline 

Station counter 

No blade force calculation 

Compute blade force 

Turbulent flow 

Laminar flow 

NFLO iteration i n subroutine FL0WIN 
Optimize KDS 
Fix KDS 

Greater than critical Reynolds number 

Less than critical Reynolds number 

Stator 

Rotor 

Propeller 
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APPENDIX A 


Sample Input Setups 

The input data to three sample cases is shown in this section to demonstrate 
how the input data must be structured in order for the PANPER analysis program to 
operate correctly. To illustrate the setup of the data input for certain modes of 
operation of the PANPER program, three different setups are shown. The three 
setups are: 

(1) Isolated Propeller Configuration 

(2) Combined Propeller-Nacelle Configuration 

(3) Combined Coaxial Propeller-Nacelle Configuration 

The job control language (JCL) and the sample data input for each of these config- 
urations are shown in Tables (II), (III), and (IV), respectively. The JCL commands 
are for a UNIVAC 1110 operating system. An isolated nacelle configuration is 
essentially the same as the second configuration (Table III) with the appropriate 
changes in the mode control input and the removal of the propeller input data. 
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JCL AND INPUT DATA SETUP FOR ISOLATED 
PROPELLER CONFIGURATION 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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JCL AND INPUT DATA SETUP FOR SINGLE 
PROPELLER-NACELLE CONFIGURATION 
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JCL AND INPUT DATA FOR COAXIAL 
PROPELLER-NACELLE CONFIGURATION 
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APPENDIX B 


Example of PANPER Analysis Program Output 


The purpose of this appendix is to give an example of the printout 
that is output from the PANPER analysis program for the case. The case 
which is used to demonstrate this printout is the analysis of a nacelle- 
propeller configuration (Table III) . This type of configuration was used 
as a test case in reference 1. The entire printout from this sample case 
will not be included in this appendix but rather specific areas of the 
printout are listed in order to illustrate the output of the various tasks 
that the program performs. Printout will be shown for the following tasks: 

(1) Main Program Heading and Output of Initial Propeller Input Data 

(2) Output’ of the Options and Input Data used in Nacelle Analysis 

(3) Output of the Nacelle Geometry 

(4) Output of the Inviscid Flow Solution (At Axial Station ZH = .76489) 

(5) Output of the Lifting Line Noninduced Inflow Conditions 

(6) Output of the Options and Reference Data used in Propeller Analysis 

(7) Output of Selected Intermediate Calculation Results for the Propeller 
Analysis 

(8) Output of the Propeller Performance Results 

(9) Output of the Propeller Blade Forces 

(10) Output of the Nacelle Viscous Flow Solution (At Axial Station ZH = .76489) 
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( 2 ) Output of the Options and Input Data used in Nacelle Analysis 
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(3) Output of the Nacelle Geometry 
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UPSTREAM STATION : H SI OoWNSIPEAM STATION 

BEGIN iNtf IiC 10 T"LO* CALCULATION 
TIME 0 77.J8B40 


( 4 ) Output of the Inviscid Flow Solution (At Axial Station ZH = . 76489 ) 
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(5) Output of the Lifting Line Noninduced Inflow Conditions 


ORIGINAL F. n Gc 

OF 


RS 


POOR QUALITY 


t-oc QOOOQOOO^h 
zoo 300000QC □□ 
*-<♦♦♦♦♦♦♦♦♦♦ I | 
r-mm 3 V> lO — «rg O 3 
m —tuo 3 uv* cw mr- oo ao 
3 3 3 0 -«rN *3 
mmr-.Ki Qp-tn k>~*o r-m 
mm CscrN (\imm 


h- oo oo o oo ooo ao 
W^OO oo ooo ooo oo 
©♦♦♦♦♦♦♦♦♦♦♦♦ 
kit- MmcM o® 
Ift «A> O ® fO O*® (\*K"> c 
fsi x* -*cvjo 3f— mmin uii/i 
3 3 ONN ®® OOO O' C* 
0*0 O OO oo OOO o o 


ooo ooaoaoaaa 
a oo coo oa coo o 

♦ ♦♦♦♦♦ ♦♦♦♦♦♦ o 

3 ccr*Nr<-^^NN^i/i 

>'<j X) J* 3 -D 45 m ul <N 

O^fO -H 3 O') CC CO 

S) 45 00 

r-r-r- c— • cc to cc cc co co cc cd 

UJ •••••* 

z o 


o o oo a o oo oo o < 
o oeoooocoao oi- 

Z ♦♦♦♦♦♦♦♦♦♦♦I/) 

O^yjl^CDOHCOvOul^X 
u.— I Omp- o 3" 3 m<M f— M v or'< 
-tZ OAJO OOO 33 3 - Ui 

-J 2rOK.M3 3 3ff3 4’3 4 , a 


3 

z o a oooo oooao ao 
o rz o oooo oo ooo oo 
—• ♦♦♦♦♦+♦+♦+♦+ 

►» JJf^T 4 f-t/'N MhjOQf- 

< r-r- o -»30 ® -ho r- 

tiOOOO crOKicnm^Kioi^ 
O'-j m — i — » mmrj r\im 3 lt ' ji ,j O 

_J JL J.O,' XU^a: 5 D_ 45 aCCL 0 C 3 


U 

ooo ao ooo o z 

0 co ooo aoooooo 

1 t i ♦ 

vC 3 iT , sOr-OJ <\jm 3 cc 3 ® ►— 
45 •—•rj cm in® (NJlo® — ® m < 
_i 3 oo r-cco — *r. ^lowT — < ►— 
— 3 3 tT- — ^U'crorSi3- > r ,£ t/1 
45 f- O' «— 1 mm «— «f\j fM CMPsi (M 


333 33 
ooo ooo oa ooo 00 

z#*+*++«+*+*+ 

*0 3 <NJ ® m 3 Of- (NiKlf— ■ m ® 
®r- ® r-O 00 m,*^* (Nj CM 

0 !> o >000 ooo 00 

O O O OO Om m^^ m 



<\J fNry (NJCNJ fMCNjfSj CNfN 

ocooooaoooo 00 

X ( 1 I I I | 1 I 1 | | 1 
CKOGOf-^COsO^OvO fMM 
cc lt r— uiCMOf^r— co«"“» ®o 
3ui ®ui p- * 0 ® ® ® a 00 
r-r- r-f-r»f^r-.r*»r-r- r-r- 


iS) 

X 

o 

hjmm mm m mm mm n mm 
►-Non oa ooo ooo ao 
+ + ♦♦♦ ♦ ♦ 
O -)(NJ 1 Mf-^Ou'KN 3 

x 30 «^^iaccwK 1 cw 1 fl 

O Otv 33 JO ®co o 00 03 

o o o o o >0 o o o o > 

3 •»••••••••«« 

o 


_jmr*> Kim mmm mmoj r\#M 

uc oa 00 o a ooao 00 

Z 3 *+ ♦♦♦♦♦♦♦♦ ♦♦ 

3 3 30 41 <W) U*1 m^ Ok/| 
mr\j m — 1 cr (MO ^ao^tNO 
O O ®3 ® fN*m NSCJ^Js 
m O Ul^co 03 NQ30S 

•O N CNfN m mm aoco 


a 

z 


er 

m 

wl 

a. 

3 


260 



UNITED TECHNOLOGIES RESEARCH CENTER ; PRESCRIBED HAKE PROP-FAN PROGRAM 
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TIME TO COMPUTE PERFORMANCE 1$ 1.21 SEC 
TIME TO COMPLETE TOTAL PR 06 RAM: <• 9 . 3 *• SECONDS 



(9) Output of the Propeller Blade Forces 
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( 10 ) Output of the Nacell Viscous Flow Solution (At Axial Station ZH = . 76489 ) 
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List of Symbols 

transformation variable 

Area, a/r 2 (dimensionless) 
r 

Van Driest constant (26.0) 
Schwartz-Christoff el parameters 
Block tridiagonal matrix (dimensionless) 
Diagonal block matrix (dimensionless) 
chord, b/r r (dimensionless) 

Left diagonal block matrix (dimensionless) 
Speed of sound (ft/sec) 

Wall friction drag coefficient 

Specific heat pressure (f t 2 /sec 2 /deg R) 

Pressure drag coefficient 

Specific heat volume (ft 2 /sec 2 /deg R) 

Friction drag 

Pressure drag 

Streamwise strain (1/sec) 

Streamwise strain, r r e ns / u r (dimensionless) 
Tangential strain (1/sec) 

Tangential strain, r r e n wu r (dimensionless) 
force/span 

Complex potential (s 4- in) , body force 
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Gap between blades (ft) 

Gap between blades, g R /r r (dimensionless) 

Duct height 

Entropy 

Mass flow (slugs/sec) 

Mass flow, m/(N B r^p r U r ) (dimensionless) 

Mach number, U/C (dimensionless) 

Mass flow/ area (slugs/f t^/sec) 

Universal mass flow parameter, m w / (P W U ) 
(dimensionless) 

Mass flow/ area, m/(P r U r ) (dimensionless) 
Normal coordinate (dimensionless) 

Number of blades (dimensionless) 

Reynolds number, r r P r U r /v r (dimensionless) 
Pressure (lb/ft^) 

Total pressure 

Turbulent Prandtl number 

Heat flux, q/(p r U r C p T r ) (dimensionless) 

Heat flux 

Radial coordinate 

Radius, r/r r (dimensionless) 

Gas constant (ft^/sec^/deg R) 

Streamwise coordinate ( dimensionless) 
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Streamwise coordinate, s/(r r V r ) 
(dimensionless) 

Stanton number (dimensionless) 

Temperature (deg R) 

Total temperature 

Magnitude of velocity or velocity 
component 

Tangential velocity 

Universal velocity, U/U* (dimensionless) 

Friction velocity, /rp"/U r (dimensionless) 

Potential flow velocity (1/V metric scale 
coeff ic ient) 

Complex coordinates in duct plan 
(r + iz) 

Distance along S and n coordinates 
Universal distance 
Axial coordinate 

Axial distance, z/r r (dimensionless) 

Swirl angle to axis (deg) 

Ratio of specific heats, Cp/C^ 
(dimensionless) 

Boundary layer thickness, 6/r r 
(dimensionless) 

Displacement thickness, 6 /r r 
(dimensionless) 

Normal coordinate, n/(r r V r ) 
(dimensionless) , or 
Transformed normal coordinate 
(dimensionless) 
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ns 


n<j> 
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Angle of streamline to axis (deg) 
Temperature ratio, T/T r (dimensionless) 
Momentum thickness, 0*/T r (dimensionless) 

v/-l 

Entropy, (I-I n )/R (dimensionless) 
von Karman constant (0.41) 

Thermal conductivity (lb/sec/deg R) 
Viscosity (slugs/f t/sec) 

3.14159 

Pressure ratio, p/p r (dimensionless) 

Density (slugs/ft^) 

Density ratio, p/p r (dimensionless) 

2 

Streamwise stress, T__/(p U ) 

ns r r 

(dimensionless) 

Tangential stress, T n< j ) /(p r Uy) 
(dimensionless) 

Stress components 

Stress, t/t w (dimensionless) 

Tangential coordinate (radians) 

Blade dissipation function 

Clauser constant (0.016) (dimensionless), 
or Transformation function, dn/dn 
(dimensionless) 

Stream function (dimensionless) 
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Matrix Operators 


T Transpose 

-1 Inverse 

Superscripts 

v Iteration number 

Mean or average quantity 

A Variables for blade force calculation 

* Deviation from mean quantity 


Subscripts 

0 

1 

2 

3 

A 

E 

H 

I 


Stagnation conditions 
Inlet conditions 
Upstream of strut blade 
Downstream of strut blade 
Adiabatic 

Effective turbulent 
Hub conditions 

Incompressible conditions, singularity 
index 


n 


in the direction of the normal 


r Reference conditions, based on stan- 

dard sea level atmosphere conditions 
for all thermodynamic quantities. The 
reference radius r r is the inlet outer 
radius, and the velocity is the mean 
inlet velocity 


274 




00 


in the streamwise direction 
Tip conditions 
Wall conditions 

Free Stream or edge of boundary layer 
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Figure 1. Flow Diagram of the Combined Propeller-Nacelle Analysis 
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Figure 2. Cylindrical and Cartesian Coordinate Systems for Propeller Geometry 
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Figure 3. Computer Wake Representation for One Blade of a Hovering Rotor, Classical 
and Generalized Distorted Wake Models 
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Figure 5. Generalized Wake Geometry Equations Containing Input Wake Constants 
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I0PT3 « 3 STRAIGHT WALL ANNULAR DIFFUSER 



I0PT3 - 5 CURVED WALL DIFFUSER NO. 1 


Figure 6. Preprogrammed Duct Wail Contours 
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Figure 7. Inlet FLow Distribution 
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Figure 8. Addition of Straight Annular Channel Inlet 
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Figure 12. Flow Diagram of Subroutine GCWAKE 
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Figure 13. Flow Diagram of Subroutine PROP 


290 


85 - 4 - 90 - 






Figure 14. Flow Diagram of Subroutine SOLVEL 
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Figure 15. Flow Diagram of Subroutine SOLVEN 
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Figure 16. Flow Chart for ALTMN 
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Figure 17. Iteration Procedure 
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Figure 18 . Flow Chart for Subroutine COORST 
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Figure 19. Calculating Ducts with Slots 
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Figure 20. Flow Chart for Subroutine Flowin 
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Figure 21. Flow Chart for GBIade 
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Figure 23. Locating Lifting Line in Steamline Coordinates 
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